Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,542,945 A→T 5.4% intergenic (+16/+17) PP_1348 → / ← pduO MutT/nudix family protein/thiamine‑phosphate pyrophosphorylase/ATP:cob(I)alamin adenosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,542,9450AT 5.4% 954.4 / 29.1 277intergenic (+16/+17)PP_1348/pduOMutT/nudix family protein/thiamine‑phosphate pyrophosphorylase/ATP:cob(I)alamin adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (141/121);  new base T (11/4);  total (152/125)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 311 total aligned reads displayed.

CGCAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGTGGC  >  NC_002947/1542796‑1543094
                                                                                                                                                     |                                                                                                                                                     
cgcAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCaaa                                                                                                                                                       <  2:2283913/150‑1 (MQ=255)
 gcAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAag                                                                                                                                                      <  2:1303108/150‑1 (MQ=255)
 gcAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAag                                                                                                                                                      >  2:4107758/1‑150 (MQ=255)
   aGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                                    <  2:748686/150‑1 (MQ=255)
   aGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                                    <  1:1001618/150‑1 (MQ=255)
   aGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                                    <  1:3935304/150‑1 (MQ=255)
       gTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGccc                                                                                                                                                >  1:4497745/1‑150 (MQ=255)
           aCTGCCGGGTTCAACTAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                            <  2:1272773/150‑1 (MQ=255)
           aCTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                            <  2:5078007/150‑1 (MQ=255)
           aCTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                            <  2:220518/150‑1 (MQ=255)
             tGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTa                                                                                                                                          <  1:4691295/150‑1 (MQ=255)
                     tCAACAAGCCTGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCTGGCGTTCTTGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGc                                                                                                                                  <  2:6747751/150‑1 (MQ=255)
                     tCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCACGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGc                                                                                                                                  <  2:6041789/150‑1 (MQ=255)
                       aacaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCt                                                                                                                                >  1:1747839/1‑150 (MQ=255)
                       aacaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCt                                                                                                                                >  2:3640506/1‑150 (MQ=255)
                       aacaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCt                                                                                                                                >  2:5314604/1‑150 (MQ=255)
                        acaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                            <  1:1396511/137‑1 (MQ=255)
                        acaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                            >  2:1396511/1‑137 (MQ=255)
                        acaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCtt                                                                                                                               >  2:4596190/1‑150 (MQ=255)
                          aaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCa                                                                                                                             >  2:4376846/1‑150 (MQ=255)
                          aaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCa                                                                                                                             <  1:3482114/150‑1 (MQ=255)
                           aGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAg                                                                                                                            <  2:3264601/150‑1 (MQ=255)
                                 tttACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTgg                                                                                                                      <  2:3667793/150‑1 (MQ=255)
                                 tttACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTgg                                                                                                                      <  1:6215140/150‑1 (MQ=255)
                                  ttACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTggg                                                                                                                     <  2:4303773/150‑1 (MQ=255)
                                   tACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTgggg                                                                                                                    <  1:2473226/150‑1 (MQ=255)
                                    aCCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTggggg                                                                                                                   >  1:6672395/1‑150 (MQ=255)
                                     ccTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggc                                                                                                                  >  1:1619408/1‑150 (MQ=255)
                                      ctgctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggcg                                                                                                                 <  1:743658/150‑1 (MQ=255)
                                        gctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggcggg                                                                                                               >  1:2614598/1‑150 (MQ=255)
                                        gctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggcggg                                                                                                               >  2:4182395/1‑150 (MQ=255)
                                          ggggtgggGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCt                                                                                                             <  1:4045681/149‑1 (MQ=255)
                                                 gggaGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCATTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCacaa                                                                                                      <  2:4805956/147‑1 (MQ=255)
                                                  gTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCacaac                                                                                                     >  2:6501194/1‑150 (MQ=255)
                                                      ggCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACtt                                                                                                 >  2:5813403/1‑150 (MQ=255)
                                                      ggCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACtt                                                                                                 >  2:2649571/1‑150 (MQ=255)
                                                         cGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTCGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTg                                                                                              <  1:6576271/150‑1 (MQ=255)
                                                          gAGTGAGCGCGTGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGc                                                                                             <  1:2750822/150‑1 (MQ=255)
                                                             tGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACCGGGTGTAGCCGGGATCCGGGCGTTCGGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAAt                                                                                          <  1:3289255/150‑1 (MQ=255)
                                                             tGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAAt                                                                                          >  1:4847506/1‑150 (MQ=255)
                                                             tGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAAt                                                                                          <  1:6329946/150‑1 (MQ=255)
                                                                  gcgAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGccct                                                                                     <  1:2720941/150‑3 (MQ=255)
                                                                  gcgAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACGTCTGCAATGCCtt                                                                                     >  2:6098004/1‑150 (MQ=255)
                                                                    gAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTgc                                                                                   <  2:3158390/150‑1 (MQ=255)
                                                                             tGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCg                                                                          >  1:910319/1‑150 (MQ=255)
                                                                                     cGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGt                                                                 <  1:1260312/150‑1 (MQ=255)
                                                                                      ggggCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGt                                                                 >  1:4032088/1‑150 (MQ=255)
                                                                                           acaGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGtgctgc                                                            >  1:4303402/1‑150 (MQ=255)
                                                                                           acaGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGtgctgc                                                            <  2:3330524/150‑1 (MQ=255)
                                                                                           acaGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGtgctgc                                                            >  1:6390643/1‑150 (MQ=255)
                                                                                             aGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCca                                                          >  2:2309517/1‑150 (MQ=255)
                                                                                                 tgtAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACaaa                                                      <  2:1012789/150‑1 (MQ=255)
                                                                                                     gCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGc                                                  >  1:1480103/1‑150 (MQ=255)
                                                                                                     gCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGc                                                  >  1:4724561/1‑150 (MQ=255)
                                                                                                       cGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAg                                                <  2:2587710/150‑1 (MQ=255)
                                                                                                          gATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTc                                             >  1:4957118/1‑150 (MQ=255)
                                                                                                              cGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGaa                                         >  2:5403768/1‑150 (MQ=255)
                                                                                                              cGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGaa                                         >  2:2968234/1‑150 (MQ=255)
                                                                                                               gggCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGaaa                                        <  2:685613/150‑1 (MQ=255)
                                                                                                               gggCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGaaa                                        >  1:3708092/1‑150 (MQ=255)
                                                                                                               gggCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGaaa                                        >  1:6027392/1‑150 (MQ=255)
                                                                                                                       tGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTg                                >  2:5657178/1‑150 (MQ=255)
                                                                                                                       tGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTg                                <  2:1588543/150‑1 (MQ=255)
                                                                                                                       tGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTg                                >  1:5000629/1‑150 (MQ=255)
                                                                                                                         gCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAACGCCGGTTGat                              >  2:346057/1‑150 (MQ=255)
                                                                                                                          cccGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGata                             >  2:5229628/1‑150 (MQ=255)
                                                                                                                           ccGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGatat                            >  1:2280341/1‑150 (MQ=255)
                                                                                                                           ccGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGatat                            >  1:5906337/1‑150 (MQ=255)
                                                                                                                            cGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGAtata                           >  1:1631883/1‑150 (MQ=255)
                                                                                                                             gAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGACCGCCGGTTGAtatag                          >  2:4512549/1‑149 (MQ=255)
                                                                                                                               ggTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGc                        <  1:805484/150‑1 (MQ=255)
                                                                                                                               ggTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGc                        <  1:1874105/150‑1 (MQ=255)
                                                                                                                               ggTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGc                        <  1:529403/150‑1 (MQ=255)
                                                                                                                                 tttGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCa                      <  1:4376846/150‑1 (MQ=255)
                                                                                                                                   tGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGa                    <  1:3640506/150‑1 (MQ=255)
                                                                                                                                      taaagggggCTGCTTTGCAGCCCCTTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCa                 >  2:2747722/6‑150 (MQ=255)
                                                                                                                                      taaagggggCTGCTTTGCAGCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCaa                >  1:2080181/6‑150 (MQ=255)
                                                                                                                                      gTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCa                 >  1:5353002/1‑150 (MQ=255)
                                                                                                                                         agggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGa                                          <  1:459220/120‑1 (MQ=255)
                                                                                                                                         aaGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACg              >  2:2093319/1‑150 (MQ=255)
                                                                                                                                         aaGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACg              >  1:5509397/1‑150 (MQ=255)
                                                                                                                                         aaGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACg              >  1:5247878/1‑150 (MQ=255)
                                                                                                                                         aaGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACg              <  2:5930079/150‑1 (MQ=255)
                                                                                                                                         aaGGGGCTGCAAAGCAGCCCCCTTTAGTACCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACg              >  1:3612008/1‑150 (MQ=255)
                                                                                                                                          gggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAGAGCCGGTTGATATACGCCAGACCa                 >  1:2839908/2‑146 (MQ=255)
                                                                                                                                          gggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAGAGCCGGTTGATATACGCCAGACCa                 <  2:2839908/145‑1 (MQ=255)
                                                                                                                                           ggggCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGcc            >  1:2956487/1‑150 (MQ=255)
                                                                                                                                            gggCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCt           >  1:6116686/1‑150 (MQ=255)
                                                                                                                                              gCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTc         <  2:3549472/150‑1 (MQ=255)
                                                                                                                                               cTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTcc        <  2:1747839/150‑1 (MQ=255)
                                                                                                                                                tGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCgat                                                                        <  2:1476398/85‑1 (MQ=255)
                                                                                                                                                tGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCa       >  2:2453176/1‑150 (MQ=255)
                                                                                                                                                 gCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAg      >  1:2689051/1‑150 (MQ=255)
                                                                                                                                                 gCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAg      >  1:5672992/1‑150 (MQ=255)
                                                                                                                                                  cAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCTCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGt     >  1:4195314/1‑150 (MQ=255)
                                                                                                                                                  cAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGt     >  2:3925175/1‑150 (MQ=255)
                                                                                                                                                  cAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGt     <  1:3245473/150‑1 (MQ=255)
                                                                                                                                                   aaaGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGTg    >  2:3754271/1‑150 (MQ=255)
                                                                                                                                                     agcagcCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCTGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGTGGc  <  1:4107758/150‑1 (MQ=255)
                                                                                                                                                     agcagcCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGTGGc  >  2:4976221/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTTCCAGTGGC  >  NC_002947/1542796‑1543094

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