Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,911,699 T→A 9.0% intergenic (+48/+80) nirD → / ← PP_1707 nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,911,6990TA 9.0% 612.2 / 2.5 221intergenic (+48/+80)nirD/PP_1707nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base T (72/115);  new base A (24/2);  total (102/119)
Fisher's exact test for biased strand distribution p-value = 1.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

TTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCC  >  NC_002947/1911550‑1911848
                                                                                                                                                     |                                                                                                                                                     
ttCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCt                                                                                                                                                       >  2:8734295/1‑150 (MQ=255)
 tCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTa                                                                                                                                                      <  1:4635041/150‑1 (MQ=255)
     ccTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGa                                                                                                                                                  >  2:4593303/1‑150 (MQ=255)
       tgctgAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGGCACCCCGCACCAACggcaa                                                                                                                                                >  1:8224096/1‑147 (MQ=255)
        gctgAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATc                                                                                                                                               <  2:7367428/150‑1 (MQ=255)
         ctgAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCt                                                                                                                                              <  1:8468353/150‑1 (MQ=255)
         ctgAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCt                                                                                                                                              <  1:2645309/150‑1 (MQ=255)
           gAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGc                                                                                                                                            <  2:1667915/150‑1 (MQ=255)
            aGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCa                                                                                                                                           >  2:7097404/1‑150 (MQ=255)
             gTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCAc                                                                                                                                          <  2:13123/150‑1 (MQ=255)
              tGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCAcc                                                                                                                                         >  2:4891931/1‑150 (MQ=255)
                   aGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCtt                                                                                                                                    <  2:5362363/150‑1 (MQ=255)
                   aGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCtt                                                                                                                                    <  1:2879396/150‑1 (MQ=255)
                   aGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCtt                                                                                                                                    <  2:7300519/150‑1 (MQ=255)
                    gtgtCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTg                                                                                                                                   >  2:4867331/1‑150 (MQ=255)
                    gtgtCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCACCTACGGAATCTGCACCGCCGTg                                                                                                                                   >  2:471042/1‑150 (MQ=255)
                       tCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTaa                                                                                                                                <  2:7154569/150‑1 (MQ=255)
                         tGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATa                                                                                                                              <  1:7325016/150‑1 (MQ=255)
                         gggaggaTGCCGGGCAGAGGGTGCGGGTGTGGCCGGTGCGGTTGCAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCGACGGAATCTGCACCGCCTTGTAATa                                                                                                                              <  2:7762916/149‑1 (MQ=255)
                          ggaggaTGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGACTCTGCACCGCCTTGTAATAc                                                                                                                             >  2:3473162/1‑150 (MQ=255)
                          ggaggaTGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCCCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATAc                                                                                                                             <  1:4394218/150‑1 (MQ=255)
                              gaTGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCAACCGAATATGCACCGCCTTGTAATACGgcg                                                                                                                         >  2:4672396/1‑150 (MQ=255)
                              gaTGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTAGCGGGCAACCCCCCCCCAACCGAAAAAGCACCCCCATGTAAAAAGgcg                                                                                                                         >  2:784558/1‑150 (MQ=255)
                                tGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGcc                                                                                                                       <  2:1937709/150‑1 (MQ=255)
                                 gCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGGCATCCCCCGCCAACCGAAAATTCACCCCCGGGTAAAAAGGCGCCg                                                                                                                      >  1:5952195/1‑150 (MQ=255)
                                  ccGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGa                                                                                                                     <  1:5632950/150‑1 (MQ=255)
                                   cGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGAc                                                                                                                    <  2:4793965/150‑1 (MQ=255)
                                   cGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGAc                                                                                                                    <  2:6498114/150‑1 (MQ=255)
                                    gggCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGAcc                                                                                                                   <  2:3515787/150‑1 (MQ=255)
                                      gCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCTACGGAAAATGCACCGCCTTGTAATACGGCGCCGACCCt                                                                                                                 >  2:7596165/1‑150 (MQ=255)
                                        agagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCtgt                                                                                                               <  1:7107042/150‑1 (MQ=255)
                                        agagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCAACGGCATCTGCACCGCCGTGTAATAAGGCGCCGccccagg                                                                                                               >  2:3991858/1‑147 (MQ=39)
                                         gagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTACACCGCCTTGTAATACGGCGCCGACCCtgtg                                                                                                              >  2:8034800/1‑150 (MQ=255)
                                         gagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCACCTACGGAAACAACACCGCCTGGTAATACGGCGCCGACCCtgtg                                                                                                              >  2:3598688/1‑150 (MQ=39)
                                          agGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCGTGTAATACGGCGCCGACCCTGTgg                                                                                                             >  2:8559093/1‑150 (MQ=255)
                                          agGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCTACCGAAACTGCACCGCCTTGTAATACCGCGCCGACCCAGTgg                                                                                                             >  2:4512215/1‑150 (MQ=255)
                                            gTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGa                                                                                                           <  1:1117081/150‑1 (MQ=255)
                                            gTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCAGTGGGa                                                                                                           >  2:4204423/1‑150 (MQ=255)
                                             ttGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAg                                                                                                          >  1:8807360/1‑150 (MQ=255)
                                               gCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCg                                                                                                        >  2:4225872/1‑150 (MQ=255)
                                                cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                       <  2:1285228/150‑1 (MQ=255)
                                                cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                       <  2:8455390/150‑1 (MQ=255)
                                                 gggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCggg                                                                                                      <  1:5013700/150‑1 (MQ=255)
                                                 gggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCggg                                                                                                      <  2:4379997/150‑1 (MQ=255)
                                                  ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                     <  2:867100/150‑1 (MQ=255)
                                                  ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                     <  1:3950076/150‑1 (MQ=255)
                                                  ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                     <  2:1519543/150‑1 (MQ=255)
                                                  ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                     <  1:2014465/150‑1 (MQ=255)
                                                  ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                     <  1:916110/150‑1 (MQ=255)
                                                   gtgtGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCt                                                                                                    <  2:6670946/150‑1 (MQ=255)
                                                    tgtgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCtt                                                                                                   <  2:6897153/150‑1 (MQ=255)
                                                    tgtgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCtt                                                                                                   <  1:6071451/150‑1 (MQ=255)
                                                       ggCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGcc                                                                                                <  1:3769857/150‑1 (MQ=255)
                                                        gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGccc                                                                                               <  2:3721266/150‑1 (MQ=255)
                                                           ggTGCGTTTGACGGGGGTGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCTGTCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAACACGGCGCCGACCGTGTGGGAGCGGGCTTGCCcgcg                                                                                            <  2:2567549/150‑1 (MQ=38)
                                                            gTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGa                                                                                           <  2:5581244/150‑1 (MQ=255)
                                                            gTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCAAGCCcgctc                                                                                           <  1:1916421/150‑3 (MQ=255)
                                                             tGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCAAGCCcgctcc                                                                                          <  2:4262075/150‑4 (MQ=255)
                                                               cGGTTGACGGGGGCGGCGGTGGAGCTGGCGGTGGCCTGTTACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCACGCCCGCTCCTACGTAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                        <  2:3910587/150‑1 (MQ=39)
                                                               cGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                        <  2:3404573/150‑1 (MQ=255)
                                                                  ttGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCATCGCGGGCCAGCCCCCGCCAAACGAATAAGCCACCCCCTGTAATAAGGCCCCCCCCCCGTGGGGGCGGGCGTTGCCGCCCAGAGGc                                                                                     >  1:1336375/1‑150 (MQ=255)
                                                                            aGGCGGTGGAGCGGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCACGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGc                                                                           <  2:5772626/150‑1 (MQ=255)
                                                                               cGGTGGAGGTGGCGGTGGCCTGATACCCCTTTGCGGCCTGGGCCGGCCTCTTCGGGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGa                                                                        <  2:6106998/150‑1 (MQ=255)
                                                                                  tGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGGGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                     <  2:1408972/150‑1 (MQ=255)
                                                                                        gggCGGTGGCCTGATACCCCTTTGCGGCCTGTGCCGGCCTCGTCGCGGGCACGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaac                                                               <  2:6294383/149‑1 (MQ=255)
                                                                                          gCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCACCGACGGAAACAGCACCGCCGTGTAAAAAGGCGCCCCCCCAGTGGGAGCCGGCGTGCCCCCCCAGAGGCCAGCACAACCCGACGAGATaaaac                                                             >  2:3089139/1‑147 (MQ=255)
                                                                                                  cTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGGAAGACCCCGCCAAACGAAAATGCACCGCCGTGTAATACGGCGCCCCCCCCGTGGGGGCGGGGTTGGCCCCGGAGAGGCCCGCCCAAGCCCCCTAGAGAAACAAacccagg                                                     >  2:52948/1‑147 (MQ=255)
                                                                                                  cTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCAGCACCGCCCTGTAATACCGCCCCGACCCGGTGGGAGCGGGCATTCCCCCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCGgc                                                     >  2:8280816/1‑150 (MQ=38)
                                                                                                  cTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCCGGCAAGCCCCCCCCAACCGAATCAGCCCCCCCCTGTAAAACGGCGCCCCCCCAGGGGGGGCGGGCGTGCCCGCGCAGAGGCCACCCCAAGCCCCCGAGAGAACCCACCCCCgc                                                     >  2:3001208/1‑150 (MQ=255)
                                                                                                     aTACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcggc                                                  <  2:2783325/150‑1 (MQ=255)
                                                                                                       cccctTTTGCTGCCTGTGCGGGCCTCTTCGGGGGCAAGCCCGCTCCACCGGAACCTGCACCCCTTTTTAAACCGCCCCCGCCCGTGTGGGCGCGGCCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGt                                                <  1:7801638/149‑1 (MQ=255)
                                                                                                       aCCCGTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGTAATCTGCACCGCCTTTTAATAGGGCGCCGACCGTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGt                                                <  2:2941234/150‑1 (MQ=39)
                                                                                                       aCCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGt                                                >  2:7911534/1‑150 (MQ=255)
                                                                                                        ccccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCGACGGAATCTGCACCCCCGTGTAATACGGCGCCGCCCCGGTGGGAGCGGGCTTGCCCGCGCACAGGCCAGCACAAGCCGACTCTATAACCAACCCCAgcggcggg                                               >  2:6951527/1‑148 (MQ=255)
                                                                                                            tttGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCGGCACCGCCTTGTACTACGGCGCCGCCCCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGcc                                           <  1:6492091/150‑1 (MQ=255)
                                                                                                              tgctgcCTGTGCCGCCCTATTCGCGGTCACGCCCGCTCAAACGGAATCTGCACCGCTTTGTAATACGGCGCCGACCCTTTGGTAGGGGGCTTGCCGGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                         <  2:3781869/150‑1 (MQ=255)
                                                                                                                   ccTGTGCCGGCCTCTTCGCGGCCAAGCCCGCTCCTACGGAATCTGCACCGCCTTTTAATACGGCGCCGCCCCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCaaa                                    <  2:3232969/150‑1 (MQ=255)
                                                                                                                         ccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATc                              <  2:1794602/150‑1 (MQ=255)
                                                                                                                          gggCCTGTCCGCGGGCACGCCCGCTCATACGGAATCTGCACCGCCTTTTAATCCGGCGCCGACCCTGGGGGCGCGGGCTTCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATcc                             <  1:6326942/149‑1 (MQ=255)
                                                                                                                           ggCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATCCGGCGCCGCCCCGGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATcca                            <  1:3600386/150‑1 (MQ=255)
                                                                                                                            gCCTCTTCGCGGGCAACCCCCCTCCTACGGAATATGCACCGCCTTGTAATACGGCGCCGCCCCTGTGGGAGCGGGCGTGCCCCCCGAGAGGCCACCACAACCCGACTAGATAACCAAACCCGGCGGCGGGCGCCCGGCCAAcacaaacac                           >  2:8324972/1‑145 (MQ=255)
                                                                                                                             ccTCTTCGCGGGCAAGCCCCCTCCAACCGAAACTGCCCCGCCGTGTAATACGGCGCCGCCCCAGTGGGAGCGGGCGTGCCCCCGAAGAGGCCCGCACAAGCCCCCTAGAAAACCAAACCCAGCGGCGGGCGCCCCGCAAACACAAacacc                          >  2:1074754/1‑150 (MQ=255)
                                                                                                                              ctctTCGCGGGCACGCCCGCTACTACGTAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGTGTGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccacca                         <  2:7619603/150‑1 (MQ=255)
                                                                                                                              ctctTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccacca                         <  2:4080821/150‑1 (MQ=255)
                                                                                                                                cttggggggCTAGCCCGTTCCTACGGTATCTGCACCGCCTTTTATTCCGGCGCCGACCCTGTGGGCGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAgg                       <  1:6570876/144‑1 (MQ=255)
                                                                                                                                   cgcgGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTg                    >  2:7136316/1‑150 (MQ=255)
                                                                                                                                   cgcgGGCAAGCCCGCTCCCACGGACTCGGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTg                    <  2:830018/150‑1 (MQ=255)
                                                                                                                                          aaGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGAc                                                                       >  1:5072274/1‑92 (MQ=255)
                                                                                                                                          aaGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacca             >  2:5491518/1‑150 (MQ=255)
                                                                                                                                            gCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacc           <  2:5830223/150‑1 (MQ=255)
                                                                                                                                              ccGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt         >  1:6725419/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt         >  1:248240/1‑150 (MQ=255)
                                                                                                                                              ccGCCCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt         >  1:3341540/1‑150 (MQ=255)
                                                                                                                                                gCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtct       >  2:2274129/1‑150 (MQ=255)
                                                                                                                                                gCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtct       >  1:4261942/1‑150 (MQ=255)
                                                                                                                                                  tCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg     >  1:475976/1‑150 (MQ=255)
                                                                                                                                                    cTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc   >  1:8114314/1‑150 (MQ=255)
                                                                                                                                                    cTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc   >  1:7995862/1‑150 (MQ=255)
                                                                                                                                                     tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc  >  2:5622728/1‑150 (MQ=255)
                                                                                                                                                     tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc  >  1:6449924/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCAGCCTGCTGAGTGGGGAGTGTCTGGAGGATGCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCC  >  NC_002947/1911550‑1911848

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: