Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,019,460 G→A V359I (GTC→ATC)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,4600GA96.8% 171.6 / NA 65V359I (GTC→ATC) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/32);  minor base C (2/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGG  >  NC_000913/2019318‑2019604
                                                                                                                                              |                                                                                                                                                
gggCATTCAGTTCGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGCCGATACGCATCCGGTTCCCAGctcta                                                                                                                                            >  2:237227/2‑148 (MQ=255)
  gCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTa                                                                                                                                          <  1:552494/149‑1 (MQ=255)
   cATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTAc                                                                                                                                         >  1:1067297/1‑149 (MQ=255)
             gggCAAAGCAACATCAGTTGCGAAAGCTTTTGTGGTCAGCAGCAGTCCGCTTCCGAGATACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGTTACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                               <  1:764034/149‑1 (MQ=255)
              ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCCACGCACAGCACACCATGAACCTCTGGCGGGGGAAGACGACGATACGCGTCCGGTGCCCAGCGCTGTACAAGGGCGTGTa                                                                                                                              >  1:394075/1‑149 (MQ=255)
              ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTa                                                                                                                              >  1:398948/1‑149 (MQ=255)
              ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCAAAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTa                                                                                                                              >  2:567419/1‑149 (MQ=255)
               gCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTaa                                                                                                                             >  1:39822/1‑149 (MQ=255)
                    ccAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAACCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGc                                                                                                                        <  2:27931/148‑1 (MQ=255)
                     cAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCa                                                                                                                       >  2:477308/1‑149 (MQ=255)
                          tCAGTGGCGAAAGCTTTAGTGGGCAGCCGCAGGACGCTTCCGATCAACTGCATAGTCAAAGCAGTGGACAGCATGAAGCTCTGGCGGGGGTAGACGTCGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                                                                                                  <  2:314604/149‑1 (MQ=255)
                             ggggCGAAAGCTTTAGTGCTCAGCAGCAGGCCACTTCCCAGAAACACCAAAGCCAACGCCCAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGATCCCATCTCTTTACAAGGGCGAGTAACAGGCAACAgcggc                                                                                                               <  2:1237769/147‑1 (MQ=255)
                                gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                            <  2:809044/149‑1 (MQ=255)
                                gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTAACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                            >  2:264614/1‑149 (MQ=255)
                                 cGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTg                                                                                                           <  1:779757/149‑1 (MQ=255)
                                  gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                                                                                          <  2:379197/149‑1 (MQ=255)
                                  gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCAACTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtggt                                                                                                          >  1:667836/1‑146 (MQ=255)
                                  gAAAGCTTTAGTGGTCAGCAGCAGGCCGCATCCCCGCAACAGCAAAGCCAACGAACAGCAAACCATGAACCTCTGGCGGGCGAAGACGACGATACGCTTCCGGATCCCATCTCTTTACAAGGGAGTGTAACAGGCAACAGCGGCGTTGa                                                                                                          >  1:1276910/1‑149 (MQ=255)
                                      gcgttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                                                                                      <  1:991068/146‑1 (MQ=255)
                                            gTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGACCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTTCAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                                                                                                >  1:217547/1‑149 (MQ=255)
                                            gTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCCTCTCTGTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAGTTTCGcc                                                                                                >  1:1183863/1‑149 (MQ=255)
                                                tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                            <  1:1180340/149‑1 (MQ=255)
                                                 cagcagcagGCCGCTTCCCAGCAACCGCAAAGCCACCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACg                                                                                           <  1:986833/149‑1 (MQ=255)
                                                  agcatcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGCCGATACGCTTCCGGTTCCCTTCTCTATACAAGGGCGTGGAACAGGCAACAGCGGCGTTGCTATTTTCTCCGAACGt                                                                                          >  1:152867/1‑149 (MQ=255)
                                                  agcagcagGCCGCTTCCCAGCAACAGCAAATCCTACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGCTACGATTCCTGTTCCCTTGTCTTTAAAAGGGCGTGTAACAGGCCACAGCGGCGAAGATATTTTCGCCAAAcct                                                                                          >  2:769865/1‑147 (MQ=255)
                                                     agcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGCCGCTACGCTTCCGGTTCCCCTCTCTTTACAAGTGCGTGTAACAGGCAACAGCGGGGTGGAGATTTTCGCCTAACGTCag                                                                                       >  1:256040/1‑149 (MQ=255)
                                                        agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg                                                                                    <  2:283185/149‑1 (MQ=255)
                                                           ccGCATCCCACCAACACCAAAGCCAACCCACAGCAAACCAAGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAACGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAg                                                                                 <  2:753889/149‑1 (MQ=255)
                                                            cGCTTCCCAGCAACAGCAACGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGc                                                                                <  2:736635/149‑1 (MQ=255)
                                                              cTTCCCTGCAACAGCAAAGCCAACGCACAGCAAACCATGATCCTCTGGCTGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCTACAGCGGCGTTGATATTTTCGCCTAACGTCTGAGGTAGCAc                                                                              <  2:1072027/149‑1 (MQ=255)
                                                              cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGGTGATAATTTCGCCTAACGTCAGAGGTAGCAc                                                                              <  1:1012444/149‑1 (MQ=255)
                                                              cTTCCCAGCAACAGCAAAGCCAACCCACAGCAAACCATGTACCTCTGGCGGGGGAAGACGACGATCCGCTTCCGGTTCCCATCTCTTTACAAGGGCGGGGACCAGGCAACAGCGGCGGTGAGATTTTCGGCTAACGTCAGAGGTAGGAc                                                                              >  2:1123034/1‑149 (MQ=255)
                                                               ttCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAcc                                                                             >  2:1059997/1‑149 (MQ=255)
                                                                  ccAGCAACAGCAAAGCCAACGCACAGCAACCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCAATCTCTTTACAAGGGCGTGTAACAGGCAACAGCCGCGTGGATATTTTCGCATAACGTCAGAGGTATCACCTTa                                                                          >  2:819215/1‑149 (MQ=255)
                                                                   cagcaagcgcaaAGCCAACGCACAGCAAACCATCAACCTCTGGCGGGGGAAGTCGACGATACGCTTCCGGTTCCCATCTCTTTACTAGGGCGTGTATCAGGCAACAGCGGCGTTGATATTTTCGCGTAACTTCTGAGGTAGCACCGTaa                                                                         <  2:1276910/149‑1 (MQ=255)
                                                                   cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGaaa                                                                         >  2:20815/1‑149 (MQ=255)
                                                                      caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGCCGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACCGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                      >  2:109421/1‑149 (MQ=255)
                                                                      caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                      <  2:59578/149‑1 (MQ=255)
                                                                          cgcaaAGCCAACGCACCGCACCCCCTGAACCTCTGGCGGGGGAAGACGACGAGACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGt                                                                  <  2:386468/148‑1 (MQ=255)
                                                                              aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTTTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGACGTAGCACCGTAATCCGCGTCatt                                                              >  1:526042/1‑149 (MQ=255)
                                                                                   aaCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTcccc                                                         >  1:264087/1‑149 (MQ=255)
                                                                                     cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGc                                                       <  2:264087/149‑1 (MQ=255)
                                                                                        acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCttt                                                    >  1:286773/1‑149 (MQ=255)
                                                                                        acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCttt                                                    >  2:412374/1‑149 (MQ=255)
                                                                                           gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                                 <  2:32645/149‑1 (MQ=255)
                                                                                            cAAACCATGAACCTCTGGCGGGGGAACACGACGATACGCTTCCGGTTCCCATCTCTATACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCAAACGTCAGAGCTAGCAACGTAATCCGCCTCTTTTACCCGCTttgttg                                                >  2:442361/1‑149 (MQ=255)
                                                                                            cAAACCATGAACCTCTGGCCGGGGCAGACGACGATACGCTTCCGGTTCCCATCTCATTACAAGGGCGTGGAACAGGCAACCGCGGCGGTCAAAATTTAGCCTAACCTCACAGGTCGCCCCGTAATCCGCTTCTATTCCCCGCTttgttg                                                >  2:1242347/1‑149 (MQ=255)
                                                                                               aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTgcgc                                             >  2:218796/1‑149 (MQ=255)
                                                                                               aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTgcgc                                             <  1:60324/149‑1 (MQ=255)
                                                                                               aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGAAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTgcgc                                             >  1:397492/1‑149 (MQ=255)
                                                                                                  aTGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCa                                          <  1:76307/149‑1 (MQ=255)
                                                                                                   ttaaCCTCTGGCGGGGGAAGACGACGATACGCTTCAGGTTACCATCTCTTTACAAGGGCGTGAATCAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCaa                                         <  1:309548/147‑1 (MQ=255)
                                                                                                     aaCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTCTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGa                                       >  2:1095996/1‑149 (MQ=255)
                                                                                                         tctGGCGTGGTACGACGACTATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAAAGGCAACAGCGGCGTTGATATTTTCGCCTATCGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGAcaca                                   <  1:769865/149‑1 (MQ=255)
                                                                                                          ctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAg                                  <  2:491010/149‑1 (MQ=255)
                                                                                                                   gAAGACGACGATACTCTTCCGGTTCCCATCTCATTACAAGGGCGTTTAACAGTCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAg                         <  2:685183/149‑1 (MQ=255)
                                                                                                                    aaGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGc                        >  1:688004/1‑149 (MQ=255)
                                                                                                                     aGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGcc                       <  2:643256/149‑1 (MQ=255)
                                                                                                                         gatgaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATa                   >  1:1127143/4‑149 (MQ=255)
                                                                                                                         gacgaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATa                   >  1:515225/1‑149 (MQ=255)
                                                                                                                            gaTACGCTTCCGGTTACCATCTCTTTACTAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTGCCCGCTTTGTTGCGCTCTAGACTCAGGATAATTAGCCGATAAGc                <  1:196955/149‑1 (MQ=255)
                                                                                                                               aCGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGt             <  1:567419/149‑1 (MQ=255)
                                                                                                                                    tCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGa        <  1:264614/149‑1 (MQ=255)
                                                                                                                                       ggTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGacac     <  1:590753/149‑1 (MQ=255)
                                                                                                                                          tCCCATCTCTGTACAGGGGCGTGTAGCAGGCCGCAGCGGCGTTGAGATTTTCGCCTAACGTCTGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGAAGGTTTATTAGCAGATAAGCAGTAGCGACACAgg  >  2:829120/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGG  >  NC_000913/2019318‑2019604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑