Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,906,818 A→G 5.6% intergenic (+116/+12) 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,906,8180AG 5.6% 294.7 / 3.7 99intergenic (+116/+12)nirD/PP_1707nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (49/44);  new base G (3/3);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTC  >  NC_002947/1906669‑1906964
                                                                                                                                                     |                                                                                                                                                  
gAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCACTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                                    <  2:298850/150‑1 (MQ=255)
   tgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCCCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGca                                                                                                                                                 >  2:1922259‑M1/53‑150 (MQ=255)
      ggAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                              >  1:408827/1‑150 (MQ=255)
      ggAGGCGGAGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTAGCGGGCAAGCCCGCTCCTCCGAAATCTGCACCACCTTCTACTACAGCGCCAACCCTGTGGGAGCGGGCTTGCCCCCGAACAGGCCAGCCCaa                                                                                                                                              >  2:281962/1‑150 (MQ=37)
             gTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                              >  1:292244/1‑143 (MQ=255)
             gTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                              <  2:292244/143‑1 (MQ=255)
              tGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                                                      <  1:1628688/150‑1 (MQ=255)
                    tGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaac                                                                                                                                <  1:1445767/150‑1 (MQ=255)
                    ggtcgttGGGCTGAGACAGCTTGGCTGGCTGTGCCAGCCTCTCCGCCGCGAAGCGCGCATCTAGGGGGCCTGCCGCGCCAGGCAATACGCCGCCGACCCGGTGGCACCGCGCTAGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaac                                                                                                                                <  2:684910/144‑1 (MQ=255)
                          tGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCCGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAAccc                                                                                                                          <  2:1252635/150‑1 (MQ=255)
                              cTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgc                                                                                                                      >  2:1336991/1‑150 (MQ=255)
                                 aTACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcggc                                                                                                                   <  1:1813426/150‑1 (MQ=255)
                                     cccTTTGCTGCCTGTGCCGGCCTCTTCGCTGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTgc                                                                                                               >  2:86842/1‑150 (MQ=255)
                                     cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGTGCGGGCTTGCCCGCGAAGAGGCCAGCGCCAGCCGACTAGATAACCAACCCCTGCGGCGTgc                                                                                                               >  1:1308879/1‑150 (MQ=255)
                                     cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaa                                                                                                                                 >  1:699080/1‑132 (MQ=255)
                                     cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaa                                                                                                                                 <  2:699080/132‑1 (MQ=255)
                                     cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCATGCGGCGTgc                                                                                                               >  1:1050822/1‑150 (MQ=255)
                                      cccttGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGcc                                                                                                                                           >  2:310323/4‑121 (MQ=255)
                                      cccttGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGcc                                                                                                                                           <  1:310323/118‑1 (MQ=255)
                                      ccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTgcg                                                                                                              >  2:556381/1‑150 (MQ=255)
                                        tttGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACGCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGcc                                                                                                            <  1:589839/150‑1 (MQ=255)
                                             tgcCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCGTTGTAATACGGCGCCGACCCTGTGGGAGCGGAGTTGCCGGCGAAGAGGCCAGCACAAGCCGACTAGATAACCACCCCCTGCGGCGTGCGCCCCGCa                                                                                                       >  1:496302/1‑150 (MQ=39)
                                                   tgctggcctcttcgcgggcaagcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaaggggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAAcaca                                                                                                 >  2:1922259‑M2/99‑150 (MQ=255)
                                                    gCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                              <  1:733501/104‑1 (MQ=255)
                                                    gCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                              >  2:733501/1‑104 (MQ=255)
                                                    gCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACAt                                                                                                <  2:408827/150‑1 (MQ=255)
                                                             TTCGCGGGCAATCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGc                                                                                       <  2:1637864‑M1/150‑1 (MQ=255)
                                                                 CGGCCAAGCCCGCTCCTACCGAATCTGCTCCGCCTTGTAATACGTCGCCGACCCTGTGGTAGCGGGCGTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGc                                                                                   >  2:496302‑M1/1‑150 (MQ=255)
                                                                      aagcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATa                                                                                                                                  <  1:352957‑M2/19‑1 (MQ=255)
                                                                      aagcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATa                                                                                                                                  >  2:352957‑M2/80‑98 (MQ=255)
                                                                      AAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                                                    <  1:635309‑M1/128‑1 (MQ=255)
                                                                      AAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                                                    >  2:635309‑M1/1‑128 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacc                                                                            >  2:1356967‑M2/78‑150 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgaagattccgtaggagcgggcttgcccgcgaagaggccGGCACAGGCCGACTAGATACCCAACCCCTGCGGCGTGCGCCCCGCCCACCCATCCACCAGGatggcaacaccc                                                                            >  1:1852866‑M2/78‑146 (MQ=255)
                                                                         cccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                                                          >  1:316056‑M2/77‑119 (MQ=255)
                                                                         cccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                                                          <  2:316056‑M2/43‑1 (MQ=255)
                                                                         CCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacca                                                                           >  1:1147083‑M1/1‑150 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccac                                                                                >  2:875838‑M2/76‑144 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccac                                                                                <  1:875838‑M2/69‑1 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                          >  1:1042272‑M2/76‑150 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGc                                                                                   <  1:652122‑M2/66‑1 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGc                                                                                   >  2:652122‑M2/76‑141 (MQ=255)
                                                                               CCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGc                                                                     >  1:556381‑M1/1‑150 (MQ=255)
                                                                                gcagggAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                                                    >  1:1523766/5‑150 (MQ=255)
                                                                                CTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGCCTTGCCTGCGCAGAGGCCAGCCCAAGCCGCCTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAGCACATCCACCAGGCTGGCGACCACCATCtcggc                                                                    <  1:280152‑M1/150‑3 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACACGCCGACTAGATaaccaac                                                                                                                            >  1:1689788/1‑93 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaaccaac                                                                                                                            <  2:1689788/93‑1 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGAGGGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCACACACAACCACCACGCTGGCCACCACCATCTCGccc                                                                   >  1:1903423/1‑150 (MQ=255)
                                                                                  aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcggc                                                                                                                   <  2:1734529/101‑1 (MQ=255)
                                                                                  aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcggc                                                                                                                   >  1:1734529/1‑101 (MQ=255)
                                                                                  aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCAAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCa                                                                  >  1:325778/1‑150 (MQ=255)
                                                                                   cagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCAt                                                                 <  2:1042272‑M2/84‑1 (MQ=255)
                                                                                         cTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgt                                                           >  1:848440/1‑150 (MQ=255)
                                                                                           gCACCGCCTTGTAGTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATACACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         <  1:883046/150‑1 (MQ=255)
                                                                                           gCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         >  2:1221549/1‑150 (MQ=255)
                                                                                            cACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcg                                                        >  1:1040772/1‑150 (MQ=255)
                                                                                             cccGCCTTGTAATACGGCGCCGACGCCGTCGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcgt                                                       <  2:1903423/149‑1 (MQ=255)
                                                                                              ccGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGcgtgcgtg                                                      >  1:246067/1‑150 (MQ=255)
                                                                                               cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                                                     >  1:1193471/1‑150 (MQ=255)
                                                                                                gCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                                                      >  2:1312848/1‑68 (MQ=255)
                                                                                                gCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                                                      <  1:1312848/68‑1 (MQ=255)
                                                                                                 ccTTGTAATACGGCGCCGACCCTGTTGGAGCGAGCTTGCCTGCGACGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTACGCCCGGCAAACACATCTACCAGGCTGGCCACCACCATCTCACCTATGCGCGTGCGTGTCt                                                   >  2:31780/1‑150 (MQ=255)
                                                                                                 ccTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCt                                                   <  1:1336991/150‑1 (MQ=255)
                                                                                                  cTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGGCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTg                                                  >  1:647174/1‑150 (MQ=255)
                                                                                                    tGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaaccaa                                                                                                                             >  2:1381906/1‑73 (MQ=255)
                                                                                                    tGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaaccaa                                                                                                                             <  1:1381906/73‑1 (MQ=255)
                                                                                                       aaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAg                                             <  1:1922259‑M2/104‑1 (MQ=255)
                                                                                                       aaTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAg                                             >  2:1278054/1‑150 (MQ=255)
                                                                                                         tACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTg                                           <  1:744033/150‑1 (MQ=255)
                                                                                                         gaCGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTg                                           <  2:1308879/149‑1 (MQ=255)
                                                                                                              cgcCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACt                                      >  1:864435/1‑150 (MQ=255)
                                                                                                                ccGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTgg                                    <  2:820614/150‑1 (MQ=255)
                                                                                                                     ccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         >  1:158948/1‑124 (MQ=255)
                                                                                                                     ccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         <  2:158948/124‑1 (MQ=255)
                                                                                                                      cTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGAt                              <  2:848440/150‑1 (MQ=255)
                                                                                                                              gcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         <  1:546647‑M2/92‑1 (MQ=255)
                                                                                                                              gcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                         >  2:546647‑M2/24‑115 (MQ=255)
                                                                                                                               cgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCa                     <  1:86842‑M2/128‑1 (MQ=255)
                                                                                                                                gggctttcccgcgaagaggccGGCACACCCCGACTAGATAGCCAACCCCAGCGGCGTGCGCCCTGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCGCATGCGCGTGCGTGTCGGTACAGTGGTACTGGCTGTTTGCGGTTGCAg                    <  2:1852866‑M2/129‑1 (MQ=255)
                                                                                                                                gggcttgcccgcgaagaggccAGCACAGGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGc                                                                                       >  1:989428‑M2/22‑83 (MQ=255)
                                                                                                                                gggcttgcccgcgaagaggccAGCACAGGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGc                                                                                       <  2:989428‑M2/62‑1 (MQ=255)
                                                                                                                                 ggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccac                                                                                            >  2:1660929‑M2/21‑77 (MQ=255)
                                                                                                                                 ggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccac                                                                                            <  1:1660929‑M2/57‑1 (MQ=255)
                                                                                                                                  gcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGa                  <  1:44132‑M2/131‑1 (MQ=255)
                                                                                                                                    ttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAgg                   >  1:1060827‑M2/18‑147 (MQ=255)
                                                                                                                                    ttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAgg                   <  2:1060827‑M2/130‑1 (MQ=255)
                                                                                                                                    ttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACt                >  1:756891‑M2/18‑150 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtctc                                                                       <  2:1220394‑M2/78‑1 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtctc                                                                       >  1:1220394‑M2/14‑91 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAg                    >  1:1366453‑M2/14‑142 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAg                    <  2:1366453‑M2/129‑1 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACAt            >  2:1903685‑M2/14‑150 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACAt            <  2:1050822‑M2/137‑1 (MQ=255)
                                                                                                                                         cgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATc           <  1:1903685‑M2/138‑1 (MQ=255)
                                                                                                                                           cgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCt                                                   >  2:990891‑M2/11‑108 (MQ=255)
                                                                                                                                           cgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCt                                                   <  1:990891‑M2/98‑1 (MQ=255)
                                                                                                                                            gaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTc        <  2:246067‑M2/141‑1 (MQ=255)
                                                                                                                                              agaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTcgc      <  2:980796‑M2/143‑1 (MQ=255)
                                                                                                                                               gaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTcgcg     <  1:1637864‑M2/144‑1 (MQ=255)
                                                                                                                                                  gccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTc  >  1:971239‑M2/4‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTC  >  NC_002947/1906669‑1906964

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: