Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,911,767 A→G 7.9% 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,911,7670AG 7.9% 578.6 / 26.6 193intergenic (+116/+12)nirD/PP_1707nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (86/87);  new base G (9/6);  total (100/93)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 196 total aligned reads displayed.

GAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTCGA  >  NC_002947/1911618‑1911915
                                                                                                                                                     |                                                                                                                                                    
   gggggCGGGGGTGGAGCGGGGGGGGGCCTGATACCCCTTTGCGGCCTGGGCCGGCCTCTTCGCGGGCACGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCAGca                                                                                                                                                   <  2:1685014/150‑1 (MQ=37)
      ggAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                                                                <  2:3214155/150‑1 (MQ=255)
         ggCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGCACAAGcc                                                                                                                                             >  2:1752958/1‑150 (MQ=255)
           ccgggcgcacgccgccggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGa                                                                                                                                           >  2:303246‑M1/45‑150 (MQ=255)
           cGGTGGAGCTGGGGGTGGCCTGATACCCCTTGGCTGCCTGGGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGGGGGAGCGGGATTGCCCGCGAAGAGGCCAGCACAAGCCGa                                                                                                                                           <  2:4855943/150‑1 (MQ=38)
              tGGGGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCGCCTACGGCATATGCACCGCCGTGTAAAAACGCGCCCACCCCGTGGGAGCGGGCGTGCCCCCGAAGAGGCCCGCACAAGCCCAcaa                                                                                                                                        >  2:5175034/1‑148 (MQ=255)
              ggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                                                        >  2:4604013‑M1/42‑150 (MQ=255)
                    tGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaac                                                                                                                                  <  2:3748358/150‑1 (MQ=255)
                             ccTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCGCCGACGGAATCTGCACCGCCGTGGAATACGGCGCCGACCCAGGGGGAGCGGGCTTGCCCGCGAAGAGGCCCCCACAAGCCGACAAGATAACCAAccccag                                                                                                                         >  2:3597007/1‑148 (MQ=38)
                               tGATACCCCTTTGCTGCCTGTGCCGGCCTCGTCGCGGGCCATCCCCCTCCTACGGGATCTTCACCGCCGTGTAAAACGGCGCCCGCCCGGTGGGAGCGGGCGTTCCCGCGCAGAGGCCAGCACAAACCGACTAGCTAACCAACCCCAgcg                                                                                                                       >  2:2932877/1‑150 (MQ=255)
                                       cTTTGCTGCCTGTGCCGGCCTCTTAGCGGGCCACCCCCCGCCGACGGGAACTGCACCCCCGGTGAAAACCGCGCCCCCCCCGGGGGGGCGGGCGTGGCCGCCGCGGGGGCCGCCCAAGCCCCCGAGAAAAACAAACCCCCCGGCGTgcgc                                                                                                               >  2:1472293/1‑150 (MQ=255)
                                        ttggcggcCTGGGCCGGCCTCTCCGCGGGCAAGCCCGCGCCGACGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGcc                                                                                                              <  2:1107861/144‑1 (MQ=38)
                                          tgctgcCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                                                            <  2:4403169/150‑1 (MQ=255)
                                          tgctgcCTGGGCGGGCCCCTCCGCGGGCAACCCCGCTCCTCCGGAATCGGCACCGCCTTGTAAGCGGGCGCCGACCCTGTGGGCGCGGGCTGCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                                                            <  1:867669/150‑1 (MQ=255)
                                               ccTGGGCCGGCCTCTTCGGGGGCAAGCCCGCACCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCaaa                                                                                                       <  1:1761245/150‑1 (MQ=255)
                                                cTGTGCCGGCCTCTTCGCGGGCCAGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAAGACGGCGCCGACCCTGGGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                                                      <  2:2850262/150‑1 (MQ=39)
                                                                  ggccacccccGCCCCCCCGGAACGTGCACCGCTTTTAACTCCGGCGCCGCCCGGGGGGGGGCGGGCGCCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                                                    <  1:2364123/143‑1 (MQ=255)
                                                                  GGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                                                    <  2:4356581‑M1/150‑1 (MQ=255)
                                                                   GGCAAGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                                                   >  2:3040426‑M1/1‑150 (MQ=255)
                                                                   GGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                                                   <  2:3217883‑M1/150‑1 (MQ=255)
                                                                     CAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacc                                                                                 <  1:2505067‑M1/150‑1 (MQ=255)
                                                                      AAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacca                                                                                <  2:5433155‑M1/150‑1 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAGGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacc                                                                              >  1:1225504‑M2/78‑150 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAGGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacc                                                                              >  2:2755105‑M2/78‑150 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                                                      >  1:2598955‑M2/78‑126 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                                                      <  2:2598955‑M2/49‑1 (MQ=255)
                                                                        gcccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacc                                                                              >  2:1012390‑M2/78‑150 (MQ=255)
                                                                         cccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACGTCCACCAGGCTGGccaccacca                                                                             >  2:1986646‑M2/77‑150 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                            >  1:1614001‑M2/76‑150 (MQ=255)
                                                                          ccgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                            >  2:2404394‑M2/76‑150 (MQ=255)
                                                                          CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                            <  2:3889762‑M1/150‑1 (MQ=255)
                                                                          CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                            <  2:2382646‑M1/150‑1 (MQ=255)
                                                                          CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                                                            <  2:2382317‑M1/150‑1 (MQ=255)
                                                                           cgctcccacagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtc                                                                           >  2:2245221‑M2/75‑150 (MQ=255)
                                                                           CGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtc                                                                           >  1:3217852‑M1/1‑150 (MQ=255)
                                                                              TCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGTGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg                                                                        <  1:5245773‑M1/150‑1 (MQ=255)
                                                                              TCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg                                                                        <  1:3257540‑M1/150‑1 (MQ=255)
                                                                              TCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg                                                                        <  1:4604572‑M1/150‑1 (MQ=255)
                                                                                CTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                                                      <  1:2614018‑M1/150‑1 (MQ=255)
                                                                                CTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                                                      <  1:4507805‑M1/150‑1 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCaaa                                                                                                       <  2:4223689/116‑1 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                                                     >  2:2057779/1‑150 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                                                     >  1:4861119/1‑150 (MQ=255)
                                                                                 tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                                                     >  1:381160/1‑150 (MQ=255)
                                                                                  acagggtcggcgccgtattacaaggcggtgcagattccgtaggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCa                                                                    <  1:4282665‑M2/83‑1 (MQ=255)
                                                                                  aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCa                                                                    >  2:2436780/1‑150 (MQ=255)
                                                                                    ggAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATg                                                                  >  2:3949055/1‑150 (MQ=255)
                                                                                        tCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGcgcg                                                              >  1:1978964/1‑150 (MQ=255)
                                                                                           gccccGCCTTTTAATCCGGCGCCCCCCCGGGGGGCGGGGGCGTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                           <  1:2500911/147‑1 (MQ=255)
                                                                                           gCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcg                                                                                                                       >  1:320255/1‑90 (MQ=255)
                                                                                             aCCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcgt                                                         <  1:5599550/150‑1 (MQ=255)
                                                                                              ccGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGcgtgcgtg                                                        >  2:4130034/1‑150 (MQ=255)
                                                                                               cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                                                       <  2:4507816/150‑1 (MQ=255)
                                                                                               cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                                                       >  2:2162170/1‑150 (MQ=255)
                                                                                                  cTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTg                                                    >  2:4027393/1‑150 (MQ=255)
                                                                                                   ttGTAATCCGGCGCCGCCCCGGGGGGAGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                                                   <  2:1233486/150‑1 (MQ=255)
                                                                                                   ttGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                                                   >  1:5121828/1‑150 (MQ=255)
                                                                                                   ttGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                                                   >  2:4961366/1‑150 (MQ=255)
                                                                                                   ttGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTga                                                   >  2:800324/1‑149 (MQ=255)
                                                                                                       aaTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAg                                               >  2:851365/1‑150 (MQ=255)
                                                                                                       aaTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAg                                               >  1:278610/1‑150 (MQ=255)
                                                                                                          aCGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTgg                                            >  1:1220994/1‑150 (MQ=255)
                                                                                                           cGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGc                                           <  2:5126409/150‑1 (MQ=255)
                                                                                                            ggCGCCCCCCCGGTGGGCGCGGGCTGCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCa                                          <  1:3918595/150‑1 (MQ=255)
                                                                                                             gcgcCGACCCTGTGGGTGGGGGGTTGCCCGCGAAGTGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACTCATCCACCAGGCTGGCCTCCTGGTTCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCAc                                         <  2:1717910/150‑1 (MQ=255)
                                                                                                                 cGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgc                                     >  2:3588150/1‑150 (MQ=255)
                                                                                                                 cGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgc                                     <  2:1600529/150‑1 (MQ=255)
                                                                                                                  gACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgcg                                    >  1:467440/1‑150 (MQ=255)
                                                                                                                     ccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGa                                 <  1:3588150/150‑1 (MQ=255)
                                                                                                                        gtgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGc                              <  2:3789685/150‑1 (MQ=255)
                                                                                                                         tgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccac                                                                                              >  1:990745/1‑85 (MQ=255)
                                                                                                                         tcttcgcgggcaagcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCg                             >  1:4931002‑M2/29‑150 (MQ=255)
                                                                                                                           ttcgcgggcaagcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGt                           >  2:2465195‑M2/27‑150 (MQ=255)
                                                                                                                           ggcgcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGt                           <  1:3156083‑M2/124‑1 (MQ=255)
                                                                                                                           ggagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGt                           >  2:3567219‑M2/27‑150 (MQ=255)
                                                                                                                            gagcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGtt                          <  2:3382856‑M2/125‑1 (MQ=255)
                                                                                                                             agcgggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTg                         >  1:1709017‑M2/25‑150 (MQ=255)
                                                                                                                               cgggctggcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCa                       <  1:3624876‑M2/128‑1 (MQ=255)
                                                                                                                                gggcttgcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAg                      <  2:3861238‑M2/129‑1 (MQ=255)
                                                                                                                                   ctggcccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGAc                   <  1:2537498‑M2/132‑1 (MQ=255)
                                                                                                                                      ccccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTAc                <  1:5175034‑M2/135‑1 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGAt                                <  2:1217801‑M2/119‑1 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGAt                                >  1:1217331‑M2/14‑132 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACAt              <  2:3322473‑M2/137‑1 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACAt              >  1:5323337‑M2/14‑150 (MQ=255)
                                                                                                                                        ccgcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACAt              >  2:56961‑M2/14‑150 (MQ=255)
                                                                                                                                          gcgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCt            >  1:3065270‑M2/12‑150 (MQ=255)
                                                                                                                                           cgaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                                                           <  2:4626130‑M2/92‑1 (MQ=255)
                                                                                                                                            gaagaggccGGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTc          <  1:1012390‑M2/141‑1 (MQ=255)
                                                                                                                                            gaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTc          <  2:5076414‑M2/141‑1 (MQ=255)
                                                                                                                                            gaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTc          <  2:939234‑M2/141‑1 (MQ=255)
                                                                                                                                            gaagaggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTc          >  2:257406‑M2/10‑150 (MQ=255)
                                                                                                                                               gaggccGGCACAGGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTcgcg       <  1:2755105‑M2/144‑1 (MQ=255)
                                                                                                                                                aggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCgcgc      <  1:2404394‑M2/145‑1 (MQ=255)
                                                                                                                                                 ggccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCt     <  1:1141391‑M2/146‑1 (MQ=255)
                                                                                                                                                  gccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTc    >  2:1289770‑M2/4‑150 (MQ=255)
                                                                                                                                                   ccAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTCg   >  1:4135280‑M2/3‑150 (MQ=255)
                                                                                                                                                    cAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTCGa  >  2:1159354‑M2/2‑150 (MQ=255)
                                                                                                                                                    cAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTCGa  <  1:2932877‑M2/149‑1 (MQ=255)
                                                                                                                                                     aGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGc                                                                                     <  2:3124374/66‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCGCGATGCGGTTGCAGGACTACATCTTCGCGCTCGA  >  NC_002947/1911618‑1911915

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: