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.

ATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGC  >  NC_000913/1977081‑1977379
                                                                                                                                                      |                                                                                                                                                     
atcaAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑Acgtcgt                                                                                                                                                       <  1:1322181/150‑1 (MQ=255)
   aaaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTAc                                                                                                                                                    >  1:830791/1‑150 (MQ=255)
    aaaaaCCGCCGAAAACTGTACCGAGAACAACCAGCTAACCTTATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGGCGTTAcc                                                                                                                                                   <  2:30047/147‑1 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTAcc                                                                                                                                                   >  1:689475/1‑150 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTAcc                                                                                                                                                   >  1:1319396/1‑150 (MQ=255)
    aaTAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTAcc                                                                                                                                                   >  2:258699/1‑150 (MQ=255)
      tAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACTACATCATCCTTCCACAGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgc                                                                                                                                                 <  1:101618/150‑1 (MQ=255)
      tAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgc                                                                                                                                                 >  2:1435430/1‑150 (MQ=255)
         ccgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCCTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCgctgc                                                                                                                                              <  1:1286243/150‑1 (MQ=255)
         ccgccgAAAACTGGAACGAGATCAACAAGGTACCCTAACAAGATCCGCACGACATATTGCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACTCTCGTTTGTGGCTATATAAAGCAGCGGTAACGACG‑TACCgctgc                                                                                                                                              <  2:848084/150‑1 (MQ=255)
          cgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgct                                                                                                                                             <  1:1153312/150‑1 (MQ=255)
          cgccgAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCAGTAACGACG‑TACCgctgct                                                                                                                                             <  1:165531/150‑1 (MQ=255)
            ccggacaCTGTACCGCGACCACCCAGGTCCCCTAATAAGCTAAGCACGACCTCCTCCTACCACTGTTGACCATGACCGGCTGGGCAGTCGGCAGGCGTTAACGCGCGATTGGGGCAAAAAACAGCAGCGGT‑ACGGCGTTACCGCTGCTgg                                                                                                                                           <  1:269624/144‑1 (MQ=255)
            ccgAAAACTGTACCGATAACAACCAGGTAACCTAATAAGATAAGCCCGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgcttt                                                                                                                                           <  1:516029/150‑3 (MQ=255)
              gAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATTACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                                                                         >  1:352307/1‑150 (MQ=255)
              gAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGaa                                                                                                                                         >  1:169962/1‑150 (MQ=255)
               aaaaCTGTACCGAGAACAACCAGGGAACCTAATAAGATAAGCACGACATCATCCTTCCAATGTTGACCATGACAGGATGTTCAGTGGGCAGGCGTTAACGCGCGATTTGGGCAAAAAAAAGCAGCGGTAACGACG‑TACCgctgctttttt                                                                                                                                        <  1:1254202/150‑6 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                >  2:88341/1‑150 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                <  2:909358/150‑1 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                >  2:1365548/1‑150 (MQ=255)
                       aCCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                                                                                <  1:184893/150‑1 (MQ=255)
                        ccGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCt                                                                                                                               >  2:970855/1‑150 (MQ=255)
                         cGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                                                                                                                              >  2:790222/1‑150 (MQ=255)
                         cGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                                                                                                                              >  2:195963/1‑150 (MQ=255)
                          gagaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTcc                                                                                                                             >  1:1536236/1‑149 (MQ=255)
                          gagaACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCa                                                                                                                             >  2:1489971/1‑150 (MQ=255)
                                aacCATGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCTATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                                                                                       <  1:997929/150‑1 (MQ=255)
                                aacCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                                                                                       >  1:694832/1‑150 (MQ=255)
                                 acCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc                                                                                                                      <  1:448261/150‑1 (MQ=255)
                                  cctggTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     <  1:1624971/147‑1 (MQ=255)
                                  ccAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     >  1:1527901/1‑150 (MQ=255)
                                  ccAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCa                                                                                                                     >  1:1358467/1‑150 (MQ=255)
                                     ggTAACCTAATAAGATAAGCACGACATAATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                                                                                  >  1:1301113/1‑150 (MQ=255)
                                      ttAACCTATTAAGATAATCACTACATCATCCTTCCACTGTTGACCATGACAGTATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGATTGTTGCGCCTCACCGTATCAGTTa                                                                                                                 <  2:1508318/149‑1 (MQ=255)
                                       aatccTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                                                                                <  2:1525010/147‑1 (MQ=255)
                                        aaCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                                                                               <  1:936302/150‑1 (MQ=255)
                                        aaCCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaaa                                                                                                               >  2:1628069/1‑150 (MQ=255)
                                           cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                                                                            >  1:176122/1‑150 (MQ=255)
                                             aataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                                                                          <  2:1275792/150‑1 (MQ=255)
                                              ataagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                                                                         >  1:1271052/1‑150 (MQ=255)
                                               taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                                                                        <  1:549989/150‑1 (MQ=255)
                                               taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                                                                        >  1:576997/1‑150 (MQ=255)
                                                    taagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGAActc                                                                                                   >  2:86925/1‑150 (MQ=255)
                                                      agCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                                                                                 <  2:1319396/150‑1 (MQ=255)
                                                       gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                                                                                >  2:1461996/1‑150 (MQ=255)
                                                             catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                                                                          <  2:1706905/150‑1 (MQ=255)
                                                                  tcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGTGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagc                                                                                     <  1:1365548/150‑1 (MQ=255)
                                                                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagca                                                                                    <  1:978558/150‑1 (MQ=255)
                                                                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACAGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACCTCCTGTACTCTCTGTTCATAcagca                                                                                    >  1:1122410/1‑150 (MQ=255)
                                                                      tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                                                                                 >  2:292389/1‑150 (MQ=255)
                                                                      tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                                                                                 >  2:17602/1‑150 (MQ=255)
                                                                         aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                                                                              <  2:694832/150‑1 (MQ=255)
                                                                           tGTTGACCATGACAGGATTGTCAGTCGTCAGGCGTTATCGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                                                                            <  2:923840/150‑1 (MQ=255)
                                                                           tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                                                                            <  1:195963/150‑1 (MQ=255)
                                                                                  cATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat                                                                     >  2:117773/1‑150 (MQ=255)
                                                                                   aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc                                                                    <  1:1322753/150‑1 (MQ=255)
                                                                                    tGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg                                                                   <  2:1261514/150‑1 (MQ=255)
                                                                                     gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg                                                                  <  2:830791/150‑1 (MQ=255)
                                                                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg                                                              >  2:324104/1‑150 (MQ=255)
                                                                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg                                                              >  1:1718459/1‑150 (MQ=255)
                                                                                           aTGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAt                                                            >  2:1335031/1‑150 (MQ=255)
                                                                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  1:790222/150‑1 (MQ=255)
                                                                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  2:1353346/150‑1 (MQ=255)
                                                                                               tCAGTCGGCAGGCGTTAAAGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCGGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTg                                                        <  1:1346128/150‑1 (MQ=255)
                                                                                                 aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTTAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTggg                                                      >  1:472110/1‑150 (MQ=255)
                                                                                                   tcgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaa                                                    <  2:382601/150‑1 (MQ=255)
                                                                                                    cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGaaa                                                   >  1:1198073/1‑150 (MQ=255)
                                                                                                       cAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGtt                                                <  1:1435430/150‑1 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTAc                                             <  1:705918/149‑1 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTAc                                             >  2:705918/1‑149 (MQ=255)
                                                                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACg                                            <  1:970855/150‑1 (MQ=255)
                                                                                                              tAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCt                                         >  1:597351/1‑150 (MQ=255)
                                                                                                                   cgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTa                                    >  1:503671/1‑150 (MQ=255)
                                                                                                                   cgcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTa                                    >  1:824568/1‑150 (MQ=255)
                                                                                                                    gcgATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTAc                                   <  2:1094970/150‑1 (MQ=255)
                                                                                                                    gcgATTGGGGCAAAAAAAAGCAGCGGG‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTAc                                   <  1:679007/150‑1 (MQ=255)
                                                                                                                          ggggCAAAAAAAAGCAGCGGT‑ACTTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGTGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCg                             <  2:471588/150‑1 (MQ=255)
                                                                                                                           gggAAAAAAAAAGCATCTGT‑ACTTCGTTACCGCTGGTTGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCATTAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCgg                            <  2:724111/150‑1 (MQ=255)
                                                                                                                            ggCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCggg                           >  1:1353093/1‑150 (MQ=255)
                                                                                                                              cAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAAAATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                         >  2:1203889/1‑150 (MQ=255)
                                                                                                                              aaaaaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAt                         <  1:258699/149‑1 (MQ=255)
                                                                                                                                 aaaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGc                      <  2:955986/150‑1 (MQ=255)
                                                                                                                                  aaaaaGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCg                     <  1:88341/150‑1 (MQ=255)
                                                                                                                                    ccagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGt                   >  2:420889/3‑150 (MQ=255)
                                                                                                                                     cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt                  >  2:124977/2‑150 (MQ=255)
                                                                                                                                     cagcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGtt                  >  2:1269324/2‑150 (MQ=255)
                                                                                                                                     aaGCAGCGGT‑ACGTCGGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGAGGGCGGtt                  <  2:1576677/150‑1 (MQ=255)
                                                                                                                                      agcagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTg                 >  1:819012/1‑150 (MQ=255)
                                                                                                                                      agcagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTg                 <  1:1628069/150‑1 (MQ=255)
                                                                                                                                        cagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTgcc               >  2:975069/1‑148 (MQ=255)
                                                                                                                                        cagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  1:614046/1‑150 (MQ=255)
                                                                                                                                        cagcGGTAACGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  2:1485404/1‑150 (MQ=255)
                                                                                                                                        cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGAGAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  1:359208/1‑150 (MQ=255)
                                                                                                                                        cagcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGAATCAGTTAAACAGCCTGTACTCTCTGTTCCTCCAGAAGTTGTGGGATAATATCGGCAGGATTCTGGCAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGAc               >  2:1548853/1‑150 (MQ=255)
                                                                                                                                         agcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTGAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACa              <  2:1358467/150‑1 (MQ=255)
                                                                                                                                         agcGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACa              <  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)
                                                                                                                                                      |                                                                                                                                                     
ATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑ACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGGGAAAGTTTACGTCTTTTTACTGCCCGGGATGGCGGTTGACATAAGCTGCAGGC  >  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: ‑