Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 968,878 G→T 8.5% intergenic (‑308/+23) PP_0831 ← / ← PP_t68 hypothetical protein/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947968,8780GT 8.5% 199.1 / 5.9 67intergenic (‑308/+23)PP_0831/PP_t68hypothetical protein/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base T (2/4);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TAACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATC  >  NC_002947/968729‑969025
                                                                                                                                                     |                                                                                                                                                    
tAACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAaccgaccg                                                                                                                                                      >  1:1168672/1‑150 (MQ=255)
            ggaaaTATATCTGACGAATGAGAGCTTAAGTGCTAGCGCAAGGGGTGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGCGAGTCAGGGGCGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGttt                                                                                                                                         <  2:130979/147‑1 (MQ=255)
            cGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGttt                                                                                                                                         >  2:1622847/1‑150 (MQ=255)
            cGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGAAGGCACAAAAAAACCGACCTTACGGTCGGttt                                                                                                                                         <  1:1106359/150‑1 (MQ=255)
             ggAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGtttt                                                                                                                                        >  1:2130095/1‑150 (MQ=255)
               aaTATATCTGACTAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCAAATAGAATGCATGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCTTACGGTCGGtttttt                                                                                                                                      <  2:961020/150‑1 (MQ=255)
                atatatCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTc                                                                                                                                     >  2:14790/1‑150 (MQ=255)
                 tatatCTGACAAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCAAATAAAATGCAGTTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACTGGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTcc                                                                                                                                    <  2:2140997/150‑1 (MQ=255)
                         aggAATGAAAGTTTAAGTGGTAGCGCACGGGGGGGGCAATAGAATGGAGGTATGGGGGGTGGCGGGGGAGTAAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTgg                                                                                                                            <  2:1168672/148‑1 (MQ=255)
                         aCGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGAttt                                                                                                                              >  2:1482227/1‑148 (MQ=255)
                         aCGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGAttt                                                                                                                              <  1:1482227/148‑1 (MQ=255)
                           gAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCGGGATATGGTg                                                                                                                          >  2:976745/1‑150 (MQ=255)
                           gAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTg                                                                                                                          >  2:1576164/1‑150 (MQ=255)
                              tGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCCAAGGCACCAAAAGACCAACCGGAAGGGTGGTGTGGTGCGGATTTGGTGccc                                                                                                                       >  2:22565/1‑150 (MQ=255)
                                  aGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAg                                                                                                                   <  1:14790/150‑1 (MQ=255)
                                  aGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAg                                                                                                                   <  1:1576164/150‑1 (MQ=255)
                                   gCTTAAGTGCTAGCGCAGGGGGCGGCAAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAgg                                                                                                                  <  2:1754992/150‑1 (MQ=255)
                                    cTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGAGGAAAAGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAggg                                                                                                                 <  1:22565/150‑1 (MQ=255)
                                     ttAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGga                                                                                                                <  2:2130095/150‑1 (MQ=255)
                                            cTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGa                                                                                                         >  1:952059/1‑150 (MQ=255)
                                             tAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACGTTACGGTCGGtttttt                                                                                                                                      <  1:1532468/120‑1 (MQ=255)
                                             tAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACGTTACGGTCGGtttttt                                                                                                                                      >  2:1532468/1‑120 (MQ=255)
                                               gcgcAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACt                                                                                                      <  1:1312149/150‑1 (MQ=255)
                                                      gggCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTc                                                                                               >  1:1029328/1‑150 (MQ=255)
                                                          ggcTAATGGAACGCATGTGTGGAGGGCGGCGCGCGAGTCAGGGGCGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCttt                                                                                           <  1:1121061/150‑1 (MQ=255)
                                                           gcCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTc                                                                                          >  2:2251374/1‑150 (MQ=255)
                                                             cAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCg                                                                                                                                   >  2:1940006/1‑107 (MQ=255)
                                                             cAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCg                                                                                                                                   <  1:1940006/107‑1 (MQ=255)
                                                               aTAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaa                                                                                      >  2:212798/1‑150 (MQ=255)
                                                                 aGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTTCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGaaatc                                                                                    <  2:944725/150‑3 (MQ=255)
                                                                   aaTGCAGGTATGGGGGGAGGCGCGGGAGACAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACgg                                                                                  >  2:74123/1‑150 (MQ=255)
                                                                     ggCAGGTAGGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAt                                                                                <  2:2340440/149‑1 (MQ=255)
                                                                       cAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAttt                                                                              >  1:1849642/1‑150 (MQ=255)
                                                                        aGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACGAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGAtttt                                                                             >  1:1750984/1‑150 (MQ=255)
                                                                        aGGTATGGGGGGAGGCGCGGGAGTCAGGAGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGCACTCCCACTCCTTTCGAAAACGGAtttt                                                                             >  2:309927/1‑150 (MQ=255)
                                                                          gTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                                    <  1:1973800/141‑1 (MQ=255)
                                                                          gTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAAc                                                                                    >  2:1973800/1‑141 (MQ=255)
                                                                          gTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGa                                                                           >  1:1058283/1‑150 (MQ=255)
                                                                          gTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGa                                                                           >  1:1337163/1‑150 (MQ=255)
                                                                           tATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACGAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCCAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                          <  2:1750984/150‑1 (MQ=255)
                                                                                            gAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCaa                                                         >  1:1695016/1‑150 (MQ=255)
                                                                                                 aGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGGAAAACAACCCAGTGAGCTAGGTCGGTTTTTTCCGGGTGGGGTGCACCGAGGGAGAATCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCg                                                    >  1:1252423/1‑150 (MQ=255)
                                                                                                  ggggggCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGc                                                   >  1:1333412/1‑150 (MQ=255)
                                                                                                            gTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCg                                         >  1:1884952/1‑150 (MQ=255)
                                                                                                             ttaaCGACGTGCGTGAATCGCAGGCACAAAAAAGCCGACCTTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGTGTCTCGAATTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGt                                        <  1:976745/148‑1 (MQ=255)
                                                                                                             tCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGt                                        <  2:952059/150‑1 (MQ=255)
                                                                                                               aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                          >  2:422516/1‑114 (MQ=255)
                                                                                                               aaCGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                          <  1:422516/114‑1 (MQ=255)
                                                                                                                     gtgcgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCa                                <  1:2251374/150‑1 (MQ=255)
                                                                                                                        cgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGGCTCTAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGta                             >  1:1976098/1‑150 (MQ=255)
                                                                                                                        cgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                          >  2:447372/1‑105 (MQ=255)
                                                                                                                        cgtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGaa                                                                          <  1:447372/105‑1 (MQ=255)
                                                                                                                         gtgAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGtat                            >  1:1474377/1‑150 (MQ=255)
                                                                                                                         gtgAATCGCAGGCACAAAAAAACCGACCATAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGtat                            <  1:538605/150‑1 (MQ=255)
                                                                                                                           gAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATaa                          >  2:693141/1‑150 (MQ=255)
                                                                                                                            aaTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGCAACCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATaaa                         <  1:54047/150‑1 (MQ=255)
                                                                                                                            aaTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATaaa                         <  2:1029328/150‑1 (MQ=255)
                                                                                                                            aaTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATaaa                         <  2:1474377/150‑1 (MQ=255)
                                                                                                                             aTCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGTGCGTGGCGCGCAGTATAAag                        <  1:74123/150‑1 (MQ=255)
                                                                                                                               cGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAagag                      <  2:1337163/150‑1 (MQ=255)
                                                                                                                                    gCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc                  <  2:1820169/149‑1 (MQ=255)
                                                                                                                                    gCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGc                  >  1:1820169/1‑149 (MQ=255)
                                                                                                                                        aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTg             >  1:961024/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGGCGCGCGCGCTATAAAGAGGTGCCGCTg             >  2:1331252/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCGGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggt           <  2:1338701/150‑1 (MQ=255)
                                                                                                                                             aaCCGACCTTACGGTCGGTTTTTTTCCGGGTTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcg         >  1:678003/1‑150 (MQ=255)
                                                                                                                                             aaCCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcgg        <  2:1849642/150‑1 (MQ=255)
                                                                                                                                             aaCCGACCGTAAGGTCGGTTTTTTCCGGATGTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTggtcgg        <  1:693141/150‑1 (MQ=255)
                                                                                                                                                   ccgTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATc  <  2:1252423/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAACAATACTTTCGGAATATATCTGACGAATGACAGCTTAAGTGCTAGCGCAAGGGGCGGCCAATAGAATGCAGGTATGGGGGGAGGCGCGGGAGTCAGGGGGGCGGCGTCAACGACGTGCGTGAATCGCAGGCACAAAAAAACCGACCGTAAGGTCGGTTTTTTCCGGATTTGGTGCCCCGAGGGAGACTCGAACTCCCACTCCTTTCGAAAACGGATTTTGAATCCGCCGCGTCTACCAATTCCGCCATCAGGGCGTGGCGCGCAGTATAAAGAGGTGCCGCTGGTCGGTCAATC  >  NC_002947/968729‑969025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: