Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,828 T→G 5.4% intergenic (+68/+17) zapB → / ← rraA cell division factor ZapB/ribonuclease E inhibitor protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,8280TG 5.4% 287.2 / 7.9 93intergenic (+68/+17)zapB/rraAcell division factor ZapB/ribonuclease E inhibitor protein A
Reads supporting (aligned to +/- strand):  ref base T (41/47);  new base G (2/3);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

TGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC  >  NC_000913/4118679‑4118973
                                                                                                                                                     |                                                                                                                                                 
tGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCt                                                                                                                                                   >  1:850500/1‑150 (MQ=255)
     acaacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCgg                                                                                                                                              <  2:739912/150‑1 (MQ=255)
      caacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                             >  1:493102/1‑150 (MQ=255)
      caacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGt                                                                                                                                             >  1:1566367/1‑150 (MQ=255)
       aacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTg                                                                                                                                            <  2:663512/150‑1 (MQ=255)
       aacCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTg                                                                                                                                            <  2:1477683/150‑1 (MQ=255)
         ccATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGcc                                                                                                                                          <  1:1182903/150‑1 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGATCGTCTGCAGGGCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttt                                                                                                                                     <  1:145323/150‑3 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCGGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                     >  1:344029/1‑150 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                     <  2:826858/150‑1 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                     >  1:1643282/1‑150 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                     <  2:1019456/150‑1 (MQ=255)
              tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                     >  2:1371839/1‑150 (MQ=255)
               gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCtt                                                                                                                                    >  1:1302748/1‑150 (MQ=255)
               gAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCtttttt                                                                                                                                    <  1:144332/150‑1 (MQ=255)
                  aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATc                                                                                                                                 <  1:1359951/150‑1 (MQ=255)
                  aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttctg                                                                                                                                 <  1:1143111/150‑4 (MQ=255)
                   gAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCa                                                                                                                                <  2:1391294/150‑1 (MQ=255)
                    aaCAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAt                                                                                                                               <  2:169396/150‑1 (MQ=255)
                        gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                           <  2:932068/150‑1 (MQ=255)
                        gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                           <  2:943085/150‑1 (MQ=255)
                        gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGATTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCaa                                                                                                                           <  2:557242/150‑1 (MQ=255)
                          agAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAata                                                                                                                         <  2:656210/150‑1 (MQ=255)
                               ggCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                    <  1:970137/150‑1 (MQ=255)
                               ggCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                    <  2:494275/150‑1 (MQ=255)
                                  tGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCgg                                                                                                                 >  2:1215531/1‑150 (MQ=255)
                                  tGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCgg                                                                                                                 <  1:352147/150‑1 (MQ=255)
                                     cAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtc                                                                                                              >  1:973164/1‑150 (MQ=255)
                                         aaCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttct                                                                                                          <  2:353078/150‑1 (MQ=255)
                                          aCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTg                                                                                                         <  2:1449520/150‑1 (MQ=255)
                                                tcAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCTCGTAATTTCTTTACTGTGTATGACATAAGGTAACGCAGAAAAAATTCACCTTGCTGTGCCTTTCTTATCATTCAATATCCAGCTGATCTTCTGAAAgta                                                                                                   <  1:916419/149‑3 (MQ=255)
                                                gCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaa                                                                                                   <  1:878721/150‑1 (MQ=255)
                                                 cAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaat                                                                                                  >  2:565066/1‑150 (MQ=255)
                                                  aGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaata                                                                                                 <  1:1146082/150‑1 (MQ=255)
                                                      ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                             <  2:102077/150‑1 (MQ=255)
                                                       ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATccc                                                                                            <  1:1084088/150‑1 (MQ=255)
                                                        tgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCg                                                                                           >  2:1659024/1‑150 (MQ=255)
                                                          ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGt                                                                                         <  2:1552554/150‑1 (MQ=255)
                                                           tgGGTAGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACATTTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGTTATTCTTAAATAATAATCCCTGTa                                                                                        <  2:234720/150‑1 (MQ=255)
                                                              gTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTg                                                                                     >  1:895084/1‑150 (MQ=255)
                                                                   attgAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAACGCACCTTGCTGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                <  2:1618417/147‑1 (MQ=255)
                                                                   aTGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                <  1:1290499/150‑1 (MQ=255)
                                                                   aTGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                >  2:158952/1‑150 (MQ=255)
                                                                    tGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCGTTCTTATCATTCAATATCCAGCGGATCTTCTGATAGAATAATCCCGGTATTGTCGGCa                                                                               >  1:1050489/1‑150 (MQ=255)
                                                                      gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                    <  1:987320/111‑1 (MQ=255)
                                                                      gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAg                                                                                                                    >  2:987320/1‑111 (MQ=255)
                                                                       aaGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATaa                                                                            <  2:1643282/150‑1 (MQ=255)
                                                                           ggTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGAt                                                                        >  2:233424/1‑150 (MQ=255)
                                                                             tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                      >  1:1361333/1‑150 (MQ=255)
                                                                              cTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGt                                                                     >  1:314856/1‑150 (MQ=255)
                                                                                   tACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                             <  2:1146559/121‑1 (MQ=255)
                                                                                   tACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                             >  1:1146559/1‑121 (MQ=255)
                                                                                    aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCgg                                                               <  1:1659024/150‑1 (MQ=255)
                                                                                      gCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGgag                                                             <  2:263991/150‑1 (MQ=255)
                                                                                       cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                               <  2:528278/131‑1 (MQ=255)
                                                                                       cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                               >  1:528278/1‑131 (MQ=255)
                                                                                            tCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                       >  1:1405624/1‑150 (MQ=255)
                                                                                            tCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                       >  1:829038/1‑150 (MQ=255)
                                                                                               tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCTGTATTGTCGTCATAAAGATTGTCGCCGGAGAAGAAGGTg                                                    <  2:973164/150‑1 (MQ=255)
                                                                                                ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                   <  2:895084/150‑1 (MQ=255)
                                                                                                ccGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                   >  1:1568404/1‑150 (MQ=255)
                                                                                                   cgTAATTTCTTTACTGTGGAGGCCAACAGGTCACGCAGAAATATGGCACCTTGCGGTGCCTTTCTTAGCATTCAATATCCAGCGGATCTTCTGAAAGAATAATGCCGGTAGTGTCGGCATAAAGATGGGCGCCGGAGAAGAAGGGTacac                                                <  2:975888/150‑1 (MQ=255)
                                                                                                         ttCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaa                                          >  2:416339/1‑150 (MQ=255)
                                                                                                          tCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                         >  2:198839/1‑150 (MQ=255)
                                                                                                             ttACTGTGTAGGCCAACTGGTATCGCAGAAAAAATGCACCTTGCGTTGCCTTTCTTATCATTCAATATCGATCTGATCTTCTGAAAGAATAATCCCTGTATTGTCTGCATAAAGATTGTCGCCGGAGAAGAAGGTGTCACCGCCAAAAtt                                      <  2:1302748/150‑1 (MQ=255)
                                                                                                                cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCAGGAGAAGAAGGTGACACCGCCAAAATTGAc                                   >  2:1042783/1‑150 (MQ=255)
                                                                                                                 tgtgTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                  <  1:565066/150‑1 (MQ=255)
                                                                                                                  gtgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                 >  2:1605808/1‑150 (MQ=255)
                                                                                                                  gtgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                 >  1:373563/1‑150 (MQ=255)
                                                                                                                  gtgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                 >  2:1492275/1‑150 (MQ=255)
                                                                                                                   tgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                <  1:3252/102‑1 (MQ=255)
                                                                                                                   tgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGc                                                                                >  2:3252/1‑102 (MQ=255)
                                                                                                                   tgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                >  1:292834/1‑150 (MQ=255)
                                                                                                                   tgtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                >  1:298602/1‑150 (MQ=255)
                                                                                                                    gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCgg                               <  2:1267526/150‑1 (MQ=255)
                                                                                                                    gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCgg                               <  2:1566367/150‑1 (MQ=255)
                                                                                                                     tAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGa                              >  2:1191822/1‑150 (MQ=255)
                                                                                                                         ccAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                          >  1:18848/1‑150 (MQ=255)
                                                                                                                          cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                         >  2:1344107/1‑150 (MQ=255)
                                                                                                                          cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                         >  2:1158253/1‑150 (MQ=255)
                                                                                                                            aCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCt                       >  1:1584445/1‑150 (MQ=255)
                                                                                                                             cAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                      <  1:1057236/150‑1 (MQ=255)
                                                                                                                             cAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                      <  1:1062885/150‑1 (MQ=255)
                                                                                                                               ggTAACGCGGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                    <  2:344029/150‑1 (MQ=255)
                                                                                                                                 tAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGc                  <  1:734937/150‑1 (MQ=255)
                                                                                                                                  aaCGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGcc                 <  2:493102/150‑1 (MQ=255)
                                                                                                                                    cGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa               >  2:594999/1‑150 (MQ=255)
                                                                                                                                      cAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg             >  2:732167/1‑150 (MQ=255)
                                                                                                                                        taagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc           >  2:1282281/5‑150 (MQ=255)
                                                                                                                                          agaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt         >  2:1234733/3‑150 (MQ=255)
                                                                                                                                          aaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATTTCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt         <  2:640593/150‑1 (MQ=255)
                                                                                                                                          aaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCtt         >  1:702721/1‑150 (MQ=255)
                                                                                                                                             aaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGc      >  2:1332693/1‑150 (MQ=255)
                                                                                                                                                 cACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc  >  1:531351/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC  >  NC_000913/4118679‑4118973

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