Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,127 T→G 5.8% Q221P (CAG→CCG)  PP_3439 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,1270TG 5.8% 224.6 / 4.7 79Q221P (CAG→CCG) PP_3439AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base G (3/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACT  >  NC_002947/3895979‑3896275
                                                                                                                                                    |                                                                                                                                                    
ttGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAAAGCCCCGGGGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTg                                                                                                                                                     >  1:936189/1‑150 (MQ=255)
  gCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgc                                                                                                                                                   <  1:658269/150‑1 (MQ=255)
  gCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgc                                                                                                                                                   <  1:598866/150‑1 (MQ=255)
    aCGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGgcgc                                                                                                                                                 >  2:1231285/1‑150 (MQ=255)
         tCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTg                                                                                                                                            >  1:1032457/1‑150 (MQ=255)
               aGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCg                                                                                                                                      >  2:1386577/1‑150 (MQ=255)
                ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     <  2:1212125/150‑1 (MQ=255)
                ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     <  2:521560/150‑1 (MQ=255)
                       ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                              >  1:1071022/1‑150 (MQ=255)
                          ggCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCa                                                                                                                           >  2:1210327/1‑150 (MQ=255)
                           gCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCAc                                                                                                                          >  1:1292393/1‑150 (MQ=255)
                               gCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACgcgc                                                                                                                      >  2:353748/1‑150 (MQ=255)
                                   cGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGt                                                                                                                  <  2:1450938/150‑1 (MQ=255)
                                       aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGAGTCATCTGCAGTTTCAGCAAGGGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAg                                                                                                              >  1:192549/1‑150 (MQ=255)
                                           cTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                          <  1:353748/150‑1 (MQ=255)
                                           cTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                          <  1:740036/150‑1 (MQ=255)
                                            tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGt                                                                                                         >  1:518428/1‑150 (MQ=255)
                                            tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGt                                                                                                         >  2:276742/1‑150 (MQ=255)
                                              tAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAg                                                                                                       >  2:1274795/1‑150 (MQ=255)
                                                   ccccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     >  2:1272338/1‑115 (MQ=255)
                                                   ccccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     <  1:1272338/115‑1 (MQ=255)
                                                    cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACt                                                                                                                                             >  1:1387407/1‑106 (MQ=255)
                                                    cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACt                                                                                                                                             <  2:1387407/106‑1 (MQ=255)
                                                    cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgc                                                                                                 <  1:1025792/150‑1 (MQ=255)
                                                       ggTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                              <  2:1292393/150‑1 (MQ=255)
                                                       ggTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                              <  2:1071022/150‑1 (MQ=255)
                                                       ggTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                              <  2:1323810/150‑1 (MQ=255)
                                                             gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGc                                                                                                                                   >  2:1121931/1‑107 (MQ=255)
                                                             gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGc                                                                                                                                   <  1:1121931/107‑1 (MQ=255)
                                                              ccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGc                                                                                       <  1:1386577/150‑1 (MQ=255)
                                                                  cAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                   >  1:396837/1‑150 (MQ=255)
                                                                       gCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGtt                                                                              <  1:1231285/150‑1 (MQ=255)
                                                                        cTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTg                                                                             <  1:1274795/150‑1 (MQ=255)
                                                                         tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAg                                                                                                           >  1:1618946/1‑119 (MQ=255)
                                                                         tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAg                                                                                                           <  2:1618946/119‑1 (MQ=255)
                                                                         tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                            <  1:1622671/150‑1 (MQ=255)
                                                                         tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                            <  1:669160/150‑1 (MQ=255)
                                                                           tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     <  1:1283845/91‑1 (MQ=255)
                                                                           tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                     >  2:1283845/1‑91 (MQ=255)
                                                                           tttGAATGCCCGCTGCAGGTGTGCACGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGATAGCCATGCGGCGTCAGGCCGAAGCGTTGGCg                                                                          >  1:1205651/1‑150 (MQ=255)
                                                                              gAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                       <  1:1210327/150‑1 (MQ=255)
                                                                               aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                           >  2:336167/1‑145 (MQ=255)
                                                                               aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                           <  1:336167/145‑1 (MQ=255)
                                                                               aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAg                                                                      >  1:1420068/1‑150 (MQ=255)
                                                                                aTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                    <  1:1155691/103‑1 (MQ=255)
                                                                                aTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                    >  2:1155691/1‑103 (MQ=255)
                                                                                 tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGc                                                                    <  2:1032457/150‑1 (MQ=255)
                                                                                             gtgtGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGGGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCCGCGTCAGGCCGAAGCGTTGGCGGACGGCACGAATCACGTa                                                        >  1:839882/1‑150 (MQ=255)
                                                                                              tgtgCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAg                                                       <  1:1311705/150‑1 (MQ=255)
                                                                                                 gCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                    >  2:1161255/1‑150 (MQ=255)
                                                                                                  ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                   >  2:1362454/1‑150 (MQ=255)
                                                                                                   cTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCg                                                  >  2:856838/1‑150 (MQ=255)
                                                                                                    tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                            <  2:174074/91‑1 (MQ=255)
                                                                                                    tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                            >  1:174074/1‑91 (MQ=255)
                                                                                                    tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                 <  1:1362454/150‑1 (MQ=255)
                                                                                                      aTCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCgagc                                               >  2:763811/1‑149 (MQ=255)
                                                                                                         aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAg                                            <  1:431322/150‑1 (MQ=255)
                                                                                                           cGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGc                                          <  1:763811/150‑1 (MQ=255)
                                                                                                            gAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                       >  2:941144/1‑120 (MQ=255)
                                                                                                            gAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                       <  1:941144/120‑1 (MQ=255)
                                                                                                            gAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGcc                                         <  2:936189/150‑1 (MQ=255)
                                                                                                                ccGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCaa                                     <  1:1161255/150‑1 (MQ=255)
                                                                                                                        ctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                             >  1:939059/1‑150 (MQ=255)
                                                                                                                        ctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                             >  2:568068/1‑150 (MQ=255)
                                                                                                                           tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                          >  2:1373258/1‑150 (MQ=255)
                                                                                                                                     gCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAgg                <  1:349298/150‑1 (MQ=255)
                                                                                                                                        aTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTg             >  1:1398942/1‑150 (MQ=255)
                                                                                                                                         ttGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                   <  2:879740/111‑1 (MQ=255)
                                                                                                                                         ttGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                   >  1:879740/1‑111 (MQ=255)
                                                                                                                                             gCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCgg        <  2:192549/150‑1 (MQ=255)
                                                                                                                                              tggCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCggg       >  1:233473/2‑150 (MQ=255)
                                                                                                                                                gCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTc     >  2:649314/1‑150 (MQ=255)
                                                                                                                                                 cccTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                    <  1:1602732/102‑1 (MQ=255)
                                                                                                                                                 cccTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                    >  2:1602732/1‑102 (MQ=255)
                                                                                                                                                   cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                    >  1:1451829/1‑100 (MQ=255)
                                                                                                                                                   cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                    <  2:1451829/100‑1 (MQ=255)
                                                                                                                                                   cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAgg                >  1:1280123/1‑136 (MQ=255)
                                                                                                                                                   cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAgg                <  2:1280123/136‑1 (MQ=255)
                                                                                                                                                   cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACt  >  1:991309/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACT  >  NC_002947/3895979‑3896275

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