Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,230 Δ1 bp 12.7% intergenic (‑91/+36) motA ← / ← flhC proton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,2300T.12.7% 288.0 / 36.9 102intergenic (‑91/+36)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base T (44/43);  new base . (7/6);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 107 total aligned reads displayed.

ATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGC  >  NC_000913/1977081‑1977379
                                                                                                                                                      |                                                                                                                                                     
atcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgtcgt                                                                                                                                                       <  1:1322181/150‑1 (MQ=255)
   aaaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAc                                                                                                                                                    >  1:830791/1‑150 (MQ=255)
    aaaaaCCGCCGAAAACTGTACCGAGAACAACCAGCTAACCTTATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGGCGTTAcc                                                                                                                                                   <  2:30047/147‑1 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAcc                                                                                                                                                   >  1:689475/1‑150 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAcc                                                                                                                                                   >  1:1319396/1‑150 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTAcc                                                                                                                                                   >  2:258699/1‑150 (MQ=255)
      tAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACTACATCATCCTTCCACAGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgc                                                                                                                                                 <  1:101618/150‑1 (MQ=255)
      tAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgc                                                                                                                                                 >  2:1435430/1‑150 (MQ=255)
         ccgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCCTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctgc                                                                                                                                              <  1:1286243/150‑1 (MQ=255)
         ccgccgAAAACTGGAACGAGATCAACAAGGTACCCTAACAAGATCCGCACGACATATTGCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACTCTCGTTTGTGGCTATATAAAGCAGCGGTAACGACGTACCgctgc                                                                                                                                              <  2:848084/150‑1 (MQ=255)
          cgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgct                                                                                                                                             <  1:1153312/150‑1 (MQ=255)
          cgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCAGTAACGACGTACCgctgct                                                                                                                                             <  1:165531/150‑1 (MQ=255)
            ccggacaCTGTACCGCGACCACCCAGGTCCCCTAATAAGCTAAGCACGACCTCCTCCTACCACTGTTGACCATGACCGGCTGGGCAGTCGGCAGGCGTTAACGCGCGATTGGGGCAAAAAACAGCAGCGGTACGGCGTTACCGCTGCTgg                                                                                                                                           <  1:269624/144‑1 (MQ=255)
            ccgAAAACTGTACCGATAACAACCAGGTAACCTAATAAGATAAGCCCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgcttt                                                                                                                                           <  1:516029/150‑3 (MQ=255)
              gAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATTACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGaa                                                                                                                                         >  1:352307/1‑150 (MQ=255)
              gAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGaa                                                                                                                                         >  1:169962/1‑150 (MQ=255)
               aaaaCTGTACCGAGAACAACCAGGGAACCTAATAAGATAAGCACGACATCATCCTTCCAATGTTGACCATGACAGGATGTTCAGTGGGCAGGCGTTAACGCGCGATTTGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgctttttt                                                                                                                                        <  1:1254202/150‑6 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                >  2:88341/1‑150 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                <  2:909358/150‑1 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                >  2:1365548/1‑150 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                <  1:184893/150‑1 (MQ=255)
                        ccGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCt                                                                                                                               >  2:970855/1‑150 (MQ=255)
                         cGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                                                                                                                              >  2:790222/1‑150 (MQ=255)
                         cGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                                                                                                                              >  2:195963/1‑150 (MQ=255)
                          gagaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTcc                                                                                                                             >  1:1536236/1‑149 (MQ=255)
                          gagaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                                                                             >  2:1489971/1‑150 (MQ=255)
                                aacCATGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCTATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                                                                                       <  1:997929/150‑1 (MQ=255)
                                aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                                                                                       >  1:694832/1‑150 (MQ=255)
                                 acCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc                                                                                                                      <  1:448261/150‑1 (MQ=255)
                                  cctggTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     <  1:1624971/147‑1 (MQ=255)
                                  ccAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     >  1:1527901/1‑150 (MQ=255)
                                  ccAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     >  1:1358467/1‑150 (MQ=255)
                                     ggTAACCTAATAAGATAAGCACGACATAATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                                                                                  >  1:1301113/1‑150 (MQ=255)
                                      ttAACCTATTAAGATAATCACTACATCATCCTTCCACTGTTGACCATGACAGTATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGATTGTTGCGCCTCACCGTATCAGTTa                                                                                                                 <  2:1508318/149‑1 (MQ=255)
                                       aatccTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                                                                                <  2:1525010/147‑1 (MQ=255)
                                        aaCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                                                                               <  1:936302/150‑1 (MQ=255)
                                        aaCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                                                                               >  2:1628069/1‑150 (MQ=255)
                                           cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                                                                            >  1:176122/1‑150 (MQ=255)
                                             aataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                                                                          <  2:1275792/150‑1 (MQ=255)
                                              ataagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                                                                         >  1:1271052/1‑150 (MQ=255)
                                               taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                                                                        <  1:549989/150‑1 (MQ=255)
                                               taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                                                                        >  1:576997/1‑150 (MQ=255)
                                                    taagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGAActc                                                                                                   >  2:86925/1‑150 (MQ=255)
                                                      agCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                                                                                 <  2:1319396/150‑1 (MQ=255)
                                                       gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                                                                                >  2:1461996/1‑150 (MQ=255)
                                                             catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                                                                          <  2:1706905/150‑1 (MQ=255)
                                                                  tcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGTGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagc                                                                                     <  1:1365548/150‑1 (MQ=255)
                                                                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                                                                                    <  1:978558/150‑1 (MQ=255)
                                                                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACAGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACCTCCTGTACTCTCTGTTCATAcagca                                                                                    >  1:1122410/1‑150 (MQ=255)
                                                                      tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                                                                                 >  2:292389/1‑150 (MQ=255)
                                                                      tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                                                                                 >  2:17602/1‑150 (MQ=255)
                                                                         aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                                                                              <  2:694832/150‑1 (MQ=255)
                                                                           tGTTGACCATGACAGGATTGTCAGTCGTCAGGCGTTATCGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                                                                            <  2:923840/150‑1 (MQ=255)
                                                                           tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                                                                            <  1:195963/150‑1 (MQ=255)
                                                                                  cATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat                                                                     >  2:117773/1‑150 (MQ=255)
                                                                                   aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc                                                                    <  1:1322753/150‑1 (MQ=255)
                                                                                    tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg                                                                   <  2:1261514/150‑1 (MQ=255)
                                                                                     gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg                                                                  <  2:830791/150‑1 (MQ=255)
                                                                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg                                                              >  2:324104/1‑150 (MQ=255)
                                                                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg                                                              >  1:1718459/1‑150 (MQ=255)
                                                                                           aTGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAt                                                            >  2:1335031/1‑150 (MQ=255)
                                                                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  1:790222/150‑1 (MQ=255)
                                                                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  2:1353346/150‑1 (MQ=255)
                                                                                               tCAGTCGGCAGGCGTTAAAGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCGGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  1:1346128/150‑1 (MQ=255)
                                                                                                 aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTTAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTggg                                                      >  1:472110/1‑150 (MQ=255)
                                                                                                   tcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaa                                                    <  2:382601/150‑1 (MQ=255)
                                                                                                    cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaaa                                                   >  1:1198073/1‑150 (MQ=255)
                                                                                                       cAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGtt                                                <  1:1435430/150‑1 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTAc                                             <  1:705918/149‑1 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTAc                                             >  2:705918/1‑149 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACg                                            <  1:970855/150‑1 (MQ=255)
                                                                                                              tAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCt                                         >  1:597351/1‑150 (MQ=255)
                                                                                                                   cgcgATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTa                                    >  1:503671/1‑150 (MQ=255)
                                                                                                                   cgcgATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTa                                    >  1:824568/1‑150 (MQ=255)
                                                                                                                    gcgATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTAc                                   <  2:1094970/150‑1 (MQ=255)
                                                                                                                    gcgATTGGGGCAAAAAAAAGCAGCGGGACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTAc                                   <  1:679007/150‑1 (MQ=255)
                                                                                                                          ggggCAAAAAAAAGCAGCGGTACTTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGTGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCg                             <  2:471588/150‑1 (MQ=255)
                                                                                                                           gggAAAAAAAAAGCATCTGTACTTCGTTACCGCTGGTTGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCATTAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCgg                            <  2:724111/150‑1 (MQ=255)
                                                                                                                            ggCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCggg                           >  1:1353093/1‑150 (MQ=255)
                                                                                                                              cAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAAAATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                         >  2:1203889/1‑150 (MQ=255)
                                                                                                                              aaaaaaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                         <  1:258699/149‑1 (MQ=255)
                                                                                                                                 aaaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGc                      <  2:955986/150‑1 (MQ=255)
                                                                                                                                  aaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCg                     <  1:88341/150‑1 (MQ=255)
                                                                                                                                    ccagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGt                   >  2:420889/3‑150 (MQ=255)
                                                                                                                                     cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt                  >  2:124977/2‑150 (MQ=255)
                                                                                                                                     cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt                  >  2:1269324/2‑150 (MQ=255)
                                                                                                                                     aaGCAGCGGTACGTCGGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAGGGCGGtt                  <  2:1576677/150‑1 (MQ=255)
                                                                                                                                      agcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTg                 >  1:819012/1‑150 (MQ=255)
                                                                                                                                      agcagcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTg                 <  1:1628069/150‑1 (MQ=255)
                                                                                                                                        cagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTgcc               >  2:975069/1‑148 (MQ=255)
                                                                                                                                        cagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  1:614046/1‑150 (MQ=255)
                                                                                                                                        cagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  2:1485404/1‑150 (MQ=255)
                                                                                                                                        cagcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGAGAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  1:359208/1‑150 (MQ=255)
                                                                                                                                        cagcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGAATCAGTTAAACAGCCTGTACTCTCTGTTCCTCCAGAAGTTGTGGGATAATATCGGCAGGATTCTGGCAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  2:1548853/1‑150 (MQ=255)
                                                                                                                                         agcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTGAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACa              <  2:1358467/150‑1 (MQ=255)
                                                                                                                                         agcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACa              <  2:689475/150‑1 (MQ=255)
                                                                                                                                                aCGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCt        <  1:1461996/150‑1 (MQ=255)
                                                                                                                                                   tcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCa     <  1:98867/150‑1 (MQ=255)
                                                                                                                                                    cgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACCGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACCTCTTTTTACTGCCCCGGATGGCGGTTGACATAAGCTGCAg    >  1:322456/1‑150 (MQ=255)
                                                                                                                                                     gtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAgg   >  2:411473/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGC  >  NC_000913/1977081‑1977379

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: