Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,270 T→A 5.3% intergenic (‑70/+99) msrA ← / ← PP_0337 peptide methionine sulfoxide reductase MsrA/sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947403,2700TA 5.3% 829.7 / 21.4 251intergenic (‑70/+99)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (96/140);  new base A (2/12);  total (99/152)
Fisher's exact test for biased strand distribution p-value = 5.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GGCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCG  >  NC_002947/403121‑403419
                                                                                                                                                       |                                                                                                                                                     
 gCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGAACATGGCCGTTCCTCTGTGTTTATTCCGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTa                                                                                                                                                      >  2:27350/1‑150 (MQ=255)
   gtgtCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACCGGTCCa                                                                                                                                                    >  1:2936962/1‑150 (MQ=255)
   gtgtCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGAACa                                                                                                                                                    <  2:3260819/150‑1 (MQ=255)
        gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                               <  2:4502986/150‑1 (MQ=255)
        gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                               <  2:5711883/150‑1 (MQ=255)
         cgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAg                                                                                                                                              >  1:6715403/1‑150 (MQ=255)
             caggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                          <  2:3235693/150‑1 (MQ=255)
             caggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                          <  1:3056519/150‑1 (MQ=255)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTc                                                                                                                                         >  1:7288419/1‑150 (MQ=255)
              aggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCACGCCCGCTCCCACAGGAACACcacaaccttt                                                                                                                                         <  2:7254555/150‑7 (MQ=39)
                gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaa                                                                                                                                       <  1:3000387/150‑1 (MQ=255)
                gcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCaa                                                                                                                                       <  1:6745260/150‑1 (MQ=255)
                  aagccaagttctaagcaacaccctgctctgcttcacctgaagcgggaaacgccgaccccatccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggcctcttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCaa                                                                                                                                       >  2:3820561‑M1/118‑150 (MQ=255)
                    ggCCTGTTCCGCAGTGGGCTTGACGTTTTTGTTCACCAGGATTTCCGATCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGt                                                                                                                                   <  1:838318/150‑1 (MQ=255)
                    ggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGt                                                                                                                                   >  2:7224729/1‑150 (MQ=255)
                    ggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGt                                                                                                                                   <  2:3336851/150‑1 (MQ=255)
                      cgtGTTCCGCAGTGGGCATGACGTTTTTGTTCTCCAGGATTTCCGAACGCAGGACCATGGCCGGTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTACAGCACAGGTTTCATAGGTtg                                                                                                                                 <  2:4674144/148‑1 (MQ=255)
                      ccTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtg                                                                                                                                 <  1:1035287/150‑1 (MQ=255)
                       cTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                <  2:7149991/150‑1 (MQ=255)
                       cTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                >  1:4446455/1‑150 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  1:4280194/150‑1 (MQ=255)
                        tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               <  2:7085498/150‑1 (MQ=255)
                                agggggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCt                                                                                                                       <  1:3930726/147‑1 (MQ=255)
                                  tGGGCATGTCGTTGTTGTTCACCAGGATTTGCGAACGCTGGAGCATGGCCGTTGCTCTGTGTTTATTCAGTGGTGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACATCTCTGGTTTCTTATGTTGTGGTGTTCCtgt                                                                                                                     <  2:5775260/150‑1 (MQ=255)
                                  tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                     <  1:6500747/150‑1 (MQ=255)
                                  tGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                     <  1:6035081/150‑1 (MQ=255)
                                   gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  2:486412/150‑1 (MQ=255)
                                   gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  1:5577425/150‑1 (MQ=255)
                                   acaccctgctctgcttcacctgacgcgggaaacgccgaccccagccggcgtttcccgcttctttttgtctgtgtgaaatgcgccggcctcttcgcgggctTGCCCGCTCCCACAGGTACAGCACAGTTCTCAaggtctgtggtgatactg                                                                                                                      >  1:1367723‑M1/101‑144 (MQ=255)
                                    ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                               >  1:4689510/1‑138 (MQ=255)
                                    ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                   <  1:7187472/150‑1 (MQ=255)
                                    ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                   >  1:4752530/1‑150 (MQ=255)
                                    ggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTgg                                                                                                                   >  2:1475307/1‑150 (MQ=255)
                                      cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGa                                                                                                                 <  2:769293/150‑1 (MQ=255)
                                      cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                 <  1:4779764/150‑1 (MQ=21)
                                        tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGc                                                                                                               <  2:5106890/150‑1 (MQ=18)
                                          aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                             <  2:3574011/150‑1 (MQ=255)
                                          aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                             <  2:1410082/150‑1 (MQ=255)
                                          aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                             <  2:3329152/150‑1 (MQ=255)
                                            gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCgg                                                                                                             <  1:2894739/150‑1 (MQ=17)
                                              ttttGTTCACCAGGATTTCCGACCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                         <  1:1940790/150‑1 (MQ=255)
                                              ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                         <  1:4865301/150‑1 (MQ=255)
                                              ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                         <  2:1145649/150‑1 (MQ=255)
                                              ttttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGGTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCCCAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                         <  2:4140532/150‑1 (MQ=255)
                                                ttGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGGTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGAGCGGGCATGc                                                                                                       <  2:5946638/150‑1 (MQ=18)
                                                ttGATCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGc                                                                                                       >  2:3207410/1‑150 (MQ=255)
                                                cttcacctgaagcgggaaacgccgaccccatccgccgtttcccgcttctatttgtctgtgtgaaatgcgccggccccttcgcgggctTGCCCGCTCCCACAGGTACAGCACACTTCTCAAGGTCTGTGGTGATACTGTGGGAGCGggcaa                                                                                                         >  2:5564302‑M1/88‑148 (MQ=255)
                                                 tGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGcc                                                                                                      <  2:4903580/150‑1 (MQ=255)
                                                  gTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGccc                                                                                                     >  2:5238100/1‑150 (MQ=255)
                                                       ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                               >  1:5064103/1‑103 (MQ=255)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                         >  2:1530814/1‑124 (MQ=255)
                                                        cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                         <  1:1530814/124‑1 (MQ=255)
                                                                ccGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGa                                                                                       >  2:4757397/1‑150 (MQ=255)
                                                                  gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCCCAGGTACAGCCCAGGGTTCGAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGact                                                                                     >  1:3446591/1‑148 (MQ=255)
                                                                   aaCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGTTCGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGTTTGTGGTTTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                    <  2:873366/150‑1 (MQ=255)
                                                                      gCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                                 >  1:4819714/1‑150 (MQ=255)
                                                                      gCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                                 >  2:3884950/1‑150 (MQ=255)
                                                                          gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                             >  2:937938/1‑150 (MQ=255)
                                                                          gACCATGACCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                             >  2:1669952/1‑150 (MQ=255)
                                                                            ccATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACACGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCCAATGGCTTCtata                                                                           >  1:5639012/1‑150 (MQ=255)
                                                                              aTGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGt                                                                         <  2:6715403/150‑1 (MQ=255)
                                                                                       ccTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGa                                                                <  2:4331223/150‑1 (MQ=255)
                                                                                        ctctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACCGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAt                                                               <  2:1605540/150‑1 (MQ=255)
                                                                                         tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                              >  2:375948/1‑150 (MQ=255)
                                                                                         tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGGATCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGCAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                              >  2:4029103/1‑150 (MQ=255)
                                                                                           tgtgTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTAc                                                            >  2:1303143/1‑150 (MQ=255)
                                                                                            gtgtTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACg                                                           <  2:228188/150‑1 (MQ=255)
                                                                                            gtgtTTATACAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATCCCCGCTCCCACAGGAACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGAGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACg                                                           >  2:67671/1‑150 (MQ=255)
                                                                                              tttttttCTGTGGGGTTGGGGTTGGGGCG‑‑TTCGGGGGCAGGCCCGCTCCCACTGGTACAGCACTGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCa                                                         <  1:7224729/149‑1 (MQ=255)
                                                                                              gtTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    <  2:2519424/91‑1 (MQ=255)
                                                                                              gtTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                    >  1:2519424/1‑91 (MQ=255)
                                                                                               tttATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATTCCCCCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCCAATGGCTTCTATAGTGCCCGAGGATTACGCAg                                                        >  1:2612799/1‑150 (MQ=255)
                                                                                                    gctggggggTCGGGGTTTGGGCC‑‑TTCGGGGGCTTGCCCGCGCCCACAGGTACCGCCCCGGTTTCCCTGGTTGGGGTGTTCCTGTGGGCGGGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTc                                                   <  1:3186380/145‑1 (MQ=255)
                                                                                                      aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACCCCACAGGTTTCAAAGGTTTTGGTGTTCCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCCGAATTGCTTCTATAGTGCCCGGGGAATACGCACCTGTCAg                                                 >  1:2100539/1‑150 (MQ=38)
                                                                                                      aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                 <  2:6380961/150‑1 (MQ=255)
                                                                                                            gTAGTGGTTTGTGCC‑‑TTCGCGGGCATGCCAGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGAGACCGAGGATTACGCAGCTGTCAGGCCAcc                                           >  2:6146853/1‑150 (MQ=38)
                                                                                                                gggTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGc                                       >  2:3206355/1‑150 (MQ=255)
                                                                                                                     tGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcgc                                  <  2:2896971/150‑1 (MQ=255)
                                                                                                                          c‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAATGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGAtc                             >  2:3963756/1‑149 (MQ=255)
                                                                                                                          c‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                             <  2:2990138/150‑1 (MQ=255)
                                                                                                                             ggagcgGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAg                            >  2:4989582/4‑150 (MQ=255)
                                                                                                                              tCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAgc                           <  2:5598519/150‑1 (MQ=255)
                                                                                                                                      aTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATg                                                                                    <  2:6929672/85‑1 (MQ=255)
                                                                                                                                        GCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCg                 <  1:5231699‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGt                                                                                                         <  1:705968‑M1/60‑1 (MQ=255)
                                                                                                                                          CCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt               <  1:6786053‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCt               <  1:6883187‑M1/150‑1 (MQ=255)
                                                                                                                                             CTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCgg            <  2:6982865‑M1/150‑1 (MQ=255)
                                                                                                                                                CCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa         <  1:6328334‑M1/150‑1 (MQ=255)
                                                                                                                                                CCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa         <  1:625023‑M1/150‑1 (MQ=255)
                                                                                                                                                CCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa         <  2:6591414‑M1/150‑1 (MQ=255)
                                                                                                                                                 CACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAg        <  1:1853340‑M1/150‑1 (MQ=255)
                                                                                                                                                  cCAGGTACAGCACTGGTTTCAATGGTTGTGGGGTTCCTGTGGGCGGGTGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGGGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  1:7255620‑M1/2‑150 (MQ=255)
                                                                                                                                                  acgggTACAGCACAGGTTTCAAAGGTTGTGGTCTTCCGGTGGGAGCGGGCGTGCCCGCGCAAGGGCCCGAATGGGGTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       >  2:5711201/4‑150 (MQ=255)
                                                                                                                                                  ACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGc                                               >  1:1273775‑M1/1‑110 (MQ=255)
                                                                                                                                                  ACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGc                                               <  2:1273775‑M1/110‑1 (MQ=255)
                                                                                                                                                  ACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       <  2:4275740‑M1/150‑1 (MQ=255)
                                                                                                                                                  ACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt       <  2:3054435‑M1/150‑1 (MQ=255)
                                                                                                                                                      gTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGc                                                                                           >  1:2059078/1‑62 (MQ=255)
                                                                                                                                                      gTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgc   >  1:3634867/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGCGTGTCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCG  >  NC_002947/403121‑403419

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: