Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,745,790 T→G 5.6% intergenic (+32/+720) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,745,7900TG 5.6% 651.1 / 13.0 224intergenic (+32/+720)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (87/118);  new base G (8/8);  total (95/129)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

GGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  NC_002947/4745642‑4745938
                                                                                                                                                    |                                                                                                                                                    
ggCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGGCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTc                                                                                                                                                     >  1:184809/1‑150 (MQ=255)
   tGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCggg                                                                                                                                                  <  2:5512958/150‑1 (MQ=255)
   tGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCggg                                                                                                                                                  >  1:2511624/1‑150 (MQ=255)
     gATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCa                                                                                                                                                <  1:5524140/150‑1 (MQ=255)
       tctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGc                                                                                                                                              <  1:5388562/150‑1 (MQ=255)
       tctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGc                                                                                                                                              <  2:1234390/150‑1 (MQ=255)
         tcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                            >  2:1359447/1‑150 (MQ=255)
         tcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                            >  2:1532667/1‑150 (MQ=255)
            cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                         <  2:4316423/150‑1 (MQ=255)
            cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                         <  2:2892069/150‑1 (MQ=255)
             aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                        <  1:924165/150‑1 (MQ=255)
             aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                        <  2:46790/150‑1 (MQ=255)
               tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTa                                                                                                                                      >  2:487626/1‑150 (MQ=255)
               tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTa                                                                                                                                      <  2:1654565/150‑1 (MQ=255)
                ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACATCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAt                                                                                                                                     <  2:83334/150‑1 (MQ=255)
                ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGTAGCGCGACATCGGTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCATGAAGGCTGCCCTCGGGCAGCCtttttttt                                                                                                                                     <  1:1639072/150‑3 (MQ=255)
                ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAt                                                                                                                                     >  2:2330559/1‑150 (MQ=255)
                ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAt                                                                                                                                     >  2:6607175/1‑150 (MQ=255)
                ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTAtattat                                                                                                                                     <  2:3442330/150‑1 (MQ=255)
                 gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAtt                                                                                                                                    >  2:1325631/1‑150 (MQ=255)
                 gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAtt                                                                                                                                    >  1:1351104/1‑150 (MQ=255)
                 gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAtt                                                                                                                                    <  1:5274387/150‑1 (MQ=255)
                  aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   >  2:2941280/1‑150 (MQ=255)
                  aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   >  2:5452756/1‑150 (MQ=255)
                  aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   <  1:6087895/150‑1 (MQ=255)
                     gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   <  2:149884/147‑1 (MQ=255)
                     gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTc                                                                                                                                <  1:6715909/150‑1 (MQ=255)
                        aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCCGCCTTTTTTATTTGTCCtg                                                                                                                             >  1:3623817/1‑150 (MQ=255)
                        aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCtg                                                                                                                             <  1:4854767/150‑1 (MQ=255)
                        aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCtg                                                                                                                             <  1:3398228/150‑1 (MQ=255)
                            tctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGcc                                                                                                                         >  2:659936/1‑150 (MQ=255)
                            tctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGcc                                                                                                                         >  1:2384741/1‑150 (MQ=255)
                            tctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGAAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTTTCCTGTGcc                                                                                                                         >  1:2384786/1‑150 (MQ=255)
                             ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                        <  2:2244670/150‑1 (MQ=255)
                              tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                        <  2:6037167/149‑1 (MQ=255)
                               ccAGCCCGTACATGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGc                                                                                                                      <  1:1633993/150‑1 (MQ=255)
                                cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGcc                                                                                                                     >  1:3165111/1‑150 (MQ=255)
                                 aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCct                                                                                                                    <  1:6777164/150‑1 (MQ=255)
                                 aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCct                                                                                                                    <  1:5988902/150‑1 (MQ=255)
                                     cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGAAAGCTGCAGTAGCCGAACGCAAGAAGCCAGCCCTCGGGCAGCCTTATTTAAATGTCCTGTGCCGGCCTCGTc                                                                                                                >  1:517419/1‑150 (MQ=255)
                                       tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTcgc                                                                                                              >  1:437070/1‑150 (MQ=255)
                                         cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCgg                                                                                                            >  1:4985241/1‑150 (MQ=255)
                                          aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc                                                                                                           <  2:5053363/150‑1 (MQ=255)
                                             aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                        <  2:776233/150‑1 (MQ=255)
                                             aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                        <  2:1958091/150‑1 (MQ=255)
                                               cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                      <  1:6165410/150‑1 (MQ=255)
                                                ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                     <  2:6690473/150‑1 (MQ=255)
                                                  ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCGCGTGGCGGCCGTGcccc                                                                                                   >  1:1175698/1‑149 (MQ=255)
                                                   cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                  <  1:1948961/150‑1 (MQ=255)
                                                   cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                  <  1:628906/150‑1 (MQ=255)
                                                     tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTGGCGGGCGTGGCcgcgc                                                                                                >  1:1743045/1‑149 (MQ=255)
                                                      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                               <  2:1900276/150‑1 (MQ=255)
                                                      ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                               <  2:6837154/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  1:4616270/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  1:2554374/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  2:4422092/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  2:3770362/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  2:100010/150‑1 (MQ=255)
                                                        gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                            <  1:6397113/150‑1 (MQ=39)
                                                         cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                            <  2:2213819/150‑1 (MQ=255)
                                                           ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                          <  1:4943665/150‑1 (MQ=255)
                                                             aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCg                                                                                        <  2:4414165/150‑1 (MQ=39)
                                                             aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                              <  2:4191903/128‑1 (MQ=255)
                                                             aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                              >  1:4191903/1‑128 (MQ=255)
                                                                 gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCGAAGAGGCCGGcac                                                                                    >  1:4453756/1‑150 (MQ=255)
                                                                  tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCGGCCCGGCCCCCGCCGACGCCGGccca                                                                                   >  1:6014375/1‑150 (MQ=255)
                                                                        ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                            <  1:2165451/135‑1 (MQ=255)
                                                                        ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTCCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGAcacc                                                                             >  1:3966696/1‑150 (MQ=255)
                                                                                 aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCGGGGCCGCGAACAGGCCGGCCCCCACACCCCaaaaaca                                                                    >  1:1779614/1‑147 (MQ=255)
                                                                                   gcgcTACATCGTGGCCCTGAAAGACCGCATGTAAGCTGCCGTAGCCTATCGCACGTTGGCTGCCCTCGGGCAGACTTTTTTATTTGTCGTGTGCCGGCCTCTTCGGGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                  <  2:2155048/150‑1 (MQ=255)
                                                                                   gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCTGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                  <  2:4540811/150‑1 (MQ=255)
                                                                                      cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCTGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGGGCCGGCCTCTTCGCGGCATTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                               <  2:3969715/150‑1 (MQ=255)
                                                                                      cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGTCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCCGCACAGAAACCACAAAAACACAgat                                                               >  1:4889896/1‑148 (MQ=255)
                                                                                                 cccTGAACGACCGCACGTACGCTGCCGTCGCCGACCGCAAGATGGCTGCCGGCGGGCTGCCTTTTTTTTTTGTCCTGGGCCGGCCTATTCGCGGCGTGGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                    <  2:4617981/150‑1 (MQ=255)
                                                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                   >  2:3203645/1‑150 (MQ=255)
                                                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCAGTTTTTCTccgcccg                                                   >  2:345533/1‑150 (MQ=255)
                                                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                   >  2:5894067/1‑150 (MQ=255)
                                                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCCGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTGGCGGGCTTGCCCGCGTAGGGGGCGGCCCAGGCCCCACAAAACCTCAGTTTTTATccgcccg                                                   >  1:3509162/1‑150 (MQ=255)
                                                                                                      aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCGAAGTGGCCGGCACCTACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCa                                               >  2:530933/1‑150 (MQ=255)
                                                                                                        aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGATTTTCTCCGCCCGCTCAcc                                             >  2:5735187/1‑149 (MQ=255)
                                                                                                             gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGcc                                        >  2:2563344/1‑150 (MQ=255)
                                                                                                                  tacgctgctGTAGCCGAACCCAAGAAGGCTGCCCTCGGCCATCATTTTTTTTTTGTCCGGGGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                   <  2:2507193/147‑1 (MQ=255)
                                                                                                                    aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGCCGTGCCCGCGAAGAGGCCGGCACAGGCACCACAAAAACTCACTTGTTTTACGCCCGCTCACGCAGCCCTTTCAg                                 >  1:1405612/1‑150 (MQ=255)
                                                                                                                     gCTGCAGTAGCCGAACGCAATAAGGCTGCCCTGGGGCTGCATTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                <  2:3657606/150‑1 (MQ=255)
                                                                                                                        gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCTGTCAGGGTa                             >  2:2181742/1‑150 (MQ=255)
                                                                                                                               ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  2:4389921/79‑1 (MQ=25)
                                                                                                                               ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCGCACCCAGCCCATTCAGGGGATTGAACg                      >  1:1060234/1‑150 (MQ=255)
                                                                                                                                cGAACGCACGACGGCGGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGGGGCGTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgg                     <  2:4823688/150‑1 (MQ=255)
                                                                                                                                cGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgg                     >  2:1079576/1‑150 (MQ=255)
                                                                                                                                 gaccGCACGACGGCGGCCCGCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCGCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                    <  2:4080773/147‑1 (MQ=255)
                                                                                                                                 gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTTTGGGCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAATAGGCCGGCACAGACACCACAACACCTCCGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                    <  2:5867431/150‑1 (MQ=255)
                                                                                                                                 gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTCGCGGCCGTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCACAGTTTTTATCCGCCCCCGCACGCAGCCCGTTAAGGGTGGTGAACGgc                    >  1:1553800/1‑150 (MQ=255)
                                                                                                                                 gAACGCAAGAAGGCGGCCCTCGGGCAGCCTTTTTTATTGGGCCTGGGCCGGCCTCTGCGCGGCCGGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                    <  2:6432279/150‑1 (MQ=255)
                                                                                                                                     gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCCCAGCCCATTCAGGGTATTGTACGgcgcct                >  1:2715854/1‑148 (MQ=255)
                                                                                                                                     gCAAGAAGGCTGCCCTCGGGCAGCCTTTTGTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCag                >  2:6836943/1‑149 (MQ=255)
                                                                                                                                       aaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc              >  2:2773498/4‑150 (MQ=255)
                                                                                                                                        aaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGGCGTGCCCGCGGAGTGGCCCGCACCGGCACCCCAAAACCGCCGTTTTTATCCGCCCCCTCACGCAGCCCGTTGCGGGTGTTGTAAGGCGCATcca             >  1:2088395/3‑150 (MQ=255)
                                                                                                                                            ggggccccTCGGCCCGACTTTTTTATTTTGCCGGGCCCGGCCTCTCGGGGGCCTGCCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac         <  1:6606809/145‑1 (MQ=255)
                                                                                                                                                  ccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt   >  2:3838505/1‑150 (MQ=255)
                                                                                                                                                   cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCGAAGAGGCCGGCACAGAAACCACAAAACCTCAGTTTTTATACGCCCGCTCACGCACCCCATTCACGGTATTTAACGGCGCATCCCCCACGAACttg  >  2:2361960/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  NC_002947/4745642‑4745938

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