Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 403,274 G→C 5.1% intergenic (‑74/+95) 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,2740GC 5.1% 782.4 / 15.9 241intergenic (‑74/+95)msrA/PP_0337peptide methionine sulfoxide reductase MsrA/sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (100/128);  new base C (2/11);  total (102/139)
Fisher's exact test for biased strand distribution p-value = 4.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

TCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCG  >  NC_002947/403127‑403423
                                                                                                                                                     |                                                                                                                                                     
  gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                                   <  1:534935/150‑1 (MQ=255)
  gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                                   <  1:6488376/150‑1 (MQ=255)
  gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATCCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcaca                                                                                                                                                   >  2:6762119/1‑150 (MQ=255)
  gcgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCGGTGTTTATTCAGTGGGGGAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGcccc                                                                                                                                                   >  2:7642369/1‑147 (MQ=255)
   cgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAg                                                                                                                                                  <  1:1201082/150‑1 (MQ=255)
   cgGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAg                                                                                                                                                  >  2:1154108/1‑150 (MQ=255)
    ggCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAgg                                                                                                                                                 <  2:3859194/150‑1 (MQ=255)
      ccAGGCAGGGCCTGGTCCGCTGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGtt                                                                                                                                               <  2:5515724/150‑1 (MQ=255)
       caggcaggGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGttt                                                                                                                                              <  2:7882950/150‑1 (MQ=255)
        aggcaggGCCTGTTCCGCAGTGGGCATGACTTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACACcacaaccttt                                                                                                                                             <  1:731007/150‑7 (MQ=255)
              ggCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCGGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGt                                                                                                                                       >  2:5244002/1‑150 (MQ=255)
               gccgggtCCGCAGTGGGCATGACGTTTTTGTTCCCCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGtt                                                                                                                                      <  1:8044963/144‑1 (MQ=255)
               gCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGtt                                                                                                                                      <  2:3243424/150‑1 (MQ=255)
                ccTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGACCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtg                                                                                                                                     <  2:2042569/150‑1 (MQ=255)
                 cTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                    >  1:3390977/1‑150 (MQ=255)
                 cTGTTCCGCAGTGGGAATGACGATTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgt                                                                                                                                    <  1:5660940/150‑1 (MQ=255)
                  tgtgCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCTTGGCCGTTCCTCTGTGTTTAGTGAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCTCAGGTTCAGCACAGGTTTCCACGGTtggg                                                                                                                                   <  2:3010718/150‑3 (MQ=255)
                  tGTTCCGCTGTGGGCATGACGTTTTTGTTCACCAGGATTTCCTAACGCAGGACCATGTCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                   <  1:4301865/150‑1 (MQ=255)
                  tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                   <  2:173622/150‑1 (MQ=255)
                  tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                   <  2:5782777/150‑1 (MQ=255)
                  tGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTtgtg                                                                                                                                   <  1:2358699/150‑1 (MQ=255)
                   gTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtgg                                                                                                                                  <  2:321759/150‑1 (MQ=255)
                    ttCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                                 <  2:6905787/150‑1 (MQ=255)
                    ttCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTgtggt                                                                                                                                 <  2:8468353/150‑1 (MQ=255)
                        gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                             <  2:6349570/150‑1 (MQ=255)
                        gCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTc                                                                                                                             <  1:330998/150‑1 (MQ=255)
                         cAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTcc                                                                                                                            >  2:4985387/1‑150 (MQ=255)
                            ggggCATGACGTTTTTGGTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgt                                                                                                                         <  2:8669179/149‑1 (MQ=255)
                             gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                        <  2:2969494/150‑1 (MQ=255)
                             gggCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCtgtg                                                                                                                        >  1:4974356/1‑150 (MQ=255)
                                cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGa                                                                                                                     <  1:3979179/150‑1 (MQ=255)
                                cATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                     <  2:6588871/150‑1 (MQ=21)
                                  tGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCCCTTTCGCGGGCACGCCCGCTCCCACAGGAACACCACAACCTTTGAAACCTGTGCTGTACCTGTGGGa                                                                                                                     <  1:602284/150‑1 (MQ=17)
                                   gACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCg                                                                                                                  >  2:3756993/1‑150 (MQ=255)
                                    aCGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCgg                                                                                                                 <  2:4147960/150‑1 (MQ=255)
                                      gTTTTTGTTCACCAGGATTTCCGAACGCAGGGCCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                               <  2:5692498/150‑1 (MQ=255)
                                      gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                               <  1:2989794/150‑1 (MQ=255)
                                      gTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGc                                                                                                               <  1:849556/150‑1 (MQ=255)
                                            gTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGccc                                                                                                         <  2:6578044/150‑1 (MQ=255)
                                            gTTCACCAGGATTCCCGAACGCCGGACCTTGGCCGTTCCTCTGTGTTTATTCCGTGGGTTTGGGGTTTGGGCC‑‑TTCGCGGGCAGGCCCGCTCCCACAGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGccc                                                                                                         <  1:945667/150‑1 (MQ=37)
                                               cACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCcgcg                                                                                                      <  2:2734420/150‑1 (MQ=255)
                                                 ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                    >  1:4449862/1‑150 (MQ=255)
                                                 ccAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCCCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                    >  1:7810292/1‑150 (MQ=255)
                                                  cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                   >  2:339367/1‑150 (MQ=255)
                                                  cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                   >  1:7176397/1‑150 (MQ=255)
                                                  cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                   >  1:1956274/1‑150 (MQ=255)
                                                  cAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                   >  2:4284852/1‑150 (MQ=255)
                                                  cAGGATTTCCGAACGCAGGACCATGGCAGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTTTTCCTGTGGGAGCGGGCGTGCCCGCGaca                                                                                                   >  1:2942555/1‑148 (MQ=255)
                                                          ccGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTTGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGa                                                                                           <  2:1787669/150‑1 (MQ=255)
                                                           cGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCCCAGGGCCCGCCCAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGaa                                                                                          >  1:190731/1‑150 (MQ=38)
                                                            gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                         <  2:8180579/150‑1 (MQ=255)
                                                            gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                         >  1:2192621/1‑150 (MQ=255)
                                                            gAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAt                                                                                         <  2:5915564/150‑1 (MQ=255)
                                                              aCGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCCCAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATgg                                                                                       >  2:2412765/1‑150 (MQ=255)
                                                                gCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                                     <  2:8983539/150‑1 (MQ=255)
                                                                gCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCGCCCACAGGTACAGCACAGGTTTCAAAGGTTGGGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCt                                                                                     >  1:817620/1‑150 (MQ=255)
                                                                  aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATTCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAAAGGCTTc                                                                                   >  1:1098482/1‑150 (MQ=255)
                                                                  aGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTc                                                                                   <  2:7497586/150‑1 (MQ=255)
                                                                    gACCTTGGCCGTTCCTCTGTGTTTATTCAGTGGGGTCGGGGTGTGGGCC‑‑TTCGCGGGCTTGCCCGCTCCCCCAGGTACAGCACAGGTTTCAAAGGTTGTGGGGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                                 <  1:9124215/150‑1 (MQ=37)
                                                                    gACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCta                                                                                 <  2:9157915/150‑1 (MQ=255)
                                                                      ccATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCtata                                                                               <  2:5487537/150‑1 (MQ=255)
                                                                        aTGGCCGTTCCTGTGTGTTTATTCTGTGGGGTCGGGGTTTGGGCC‑‑TTCGCGGGCAGGCCCGCTCCCACTGGTACAGCACAGGTTTCAACGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGt                                                                             <  1:3680073/150‑1 (MQ=38)
                                                                          ggCCGTTCCTCTGTGTTTAGTCCGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGGACAGCACAGGTTTCAACGGGTGGGGTGTTCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGc                                                                           <  1:5130373/150‑1 (MQ=38)
                                                                              gTTCCTCTGTTTTTATACAGTGGGGTAGGGGATTGGCCC‑‑TTCGCGGGCATGCCCGCACCCACAGGTACACCACAGCTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGa                                                                       <  1:8358867/150‑1 (MQ=38)
                                                                               ttCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAg                                                                      >  2:4585392/1‑150 (MQ=255)
                                                                                   tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGAtt                                                                  >  1:4636834/1‑150 (MQ=255)
                                                                                   tctGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCAAGCCCGCTCCCACACGTACAGCACAGGTGTCAAAGGTTGTGGTGTTCCTGTGGGACCGGGCGTGCCCGCCAAAGGGCCCGAATGGCTTCAATAGTGCCCCAGGaaa                                                                  >  1:4500233/1‑148 (MQ=38)
                                                                                       tgtTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑GTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGc                                                              <  2:6504615/150‑1 (MQ=255)
                                                                                        gtTTATTCCGTGGGGTCGGGGTTGGGGCC‑‑TTGGCGGGCAGGCCCGCTCCCACTGTTACAGCACTGTTTTCAATGGTTGTGGTGTTCGTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCa                                                             <  1:6301691/150‑1 (MQ=255)
                                                                                            aTTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTg                                                         >  1:3675787/1‑150 (MQ=255)
                                                                                             ttCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCAAGCCCCCGCCCACACGTACAGCACAGGTTTCAAAGGTTGTGGTGTTACTGTGGGAGCGGGCGTGCCCCCCAAACGGCCCCAATGGCTTCTATAGTGCCCCAGGAGTACGCAGCAGt                                                        >  1:902222/1‑150 (MQ=255)
                                                                                              tCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCAGTc                                                       >  1:8053800/1‑150 (MQ=255)
                                                                                               cAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCa                                                      <  2:4386410/150‑1 (MQ=255)
                                                                                                cggggggttggggttTGGGCC‑‑GTCGGGGGCTGGCCCGCGCCCACAGGTACCGCACTGGTTTAAATGGTTGTGGGGTGCCTGTGGTCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                     <  1:529882/147‑1 (MQ=255)
                                                                                                aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                     <  2:1629310/150‑1 (MQ=255)
                                                                                                aGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGGACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAg                                                     >  1:6345546/1‑150 (MQ=255)
                                                                                                aGTGGGGTAGGGGTTTGGGCC‑‑TTCCCGGGCAAGCCCCCGCCCACACGTACACCACAGGTTTAAAAGGTTGTGGTGTTACAGTGGGAGCGGGCGTGCCCCCCACAGGGCCCCAAAGGCGTCAAAAGTGCCCGCGGAATAACCAGCAGTCac                                                     >  1:1565498/1‑149 (MQ=255)
                                                                                                      gttggggTTTGGGCC‑‑TTCGCGGGCATGCCGGCTCCCACTGGTACAGCACAGGTTTCAATGGTTGTGGTGTTCGTGTGGGTGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCAcc                                               <  1:725651/147‑1 (MQ=39)
                                                                                                               ggggCC‑‑TTCGCGGGCAGGCCCGCTCCCACAGGAACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcgc                                      <  1:8708904/149‑1 (MQ=255)
                                                                                                                gggCC‑‑TGCGCGGGCAGGCCCGCTCCCACAGGTACCGCACAGGTTTCACAGGTTGTGGTGTTCCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCAcgcg                                     <  1:6161531/150‑1 (MQ=255)
                                                                                                                  gCC‑‑TGCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTGGTGGTGTTCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCaatgagc                                   <  2:1823940/150‑7 (MQ=38)
                                                                                                                    c‑‑TTCGGGGGAAGGCCCGCTCCCACCGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                                 <  1:6900877/150‑1 (MQ=255)
                                                                                                                    c‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATa                                 <  2:7345625/150‑1 (MQ=255)
                                                                                                                                aTGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGc                       <  2:2663409/150‑1 (MQ=255)
                                                                                                                                     cGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTc                  <  2:5144093/150‑1 (MQ=255)
                                                                                                                                      gCTCCCCCCGGTACAGCCCCGGTTCCCAGGGTGGTGGGGTTCGGGGGGGCGGGGGGGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCg                 <  1:7716779/150‑1 (MQ=255)
                                                                                                                                        tCCCACAGGTACAGCACCGGGTTCAAAGGGTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc               >  2:2217808/1‑150 (MQ=255)
                                                                                                                                        tCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGc               >  2:2572521/1‑150 (MQ=255)
                                                                                                                                          ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa             >  1:6741085/1‑150 (MQ=255)
                                                                                                                                          ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa             >  2:8399346/1‑150 (MQ=255)
                                                                                                                                          ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa             <  2:8317616/150‑1 (MQ=255)
                                                                                                                                          ccACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCa             >  2:7466025/1‑150 (MQ=255)
                                                                                                                                            acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCGGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt           >  2:6991675/1‑150 (MQ=255)
                                                                                                                                            acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt           >  1:1333977/1‑150 (MQ=255)
                                                                                                                                            acaGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGt           >  2:6991253/1‑150 (MQ=255)
                                                                                                                                             caGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGAt                                  <  2:6287389/126‑1 (MQ=255)
                                                                                                                                              aGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTc         <  2:7721763/150‑1 (MQ=255)
                                                                                                                                                 tACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTcgcg      >  2:8808021/1‑150 (MQ=255)
                                                                                                                                                    aGCACAGGTTTCAAAGGTTGTGGTGTTCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGccc   <  2:1268359/150‑1 (MQ=255)
                                                                                                                                                     gCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCg  >  1:7973041/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCGGCCAGGCAGGGCCTGTTCCGCAGTGGGCATGACGTTTTTGTTCACCAGGATTTCCGAACGCAGGACCATGGCCGTTCCTCTGTGTTTATTCAGTGGGGTAGGGGTTTGGGCC‑‑TTCGCGGGCATGCCCGCTCCCACAGGTACAGCACAGGTTTCAAAGGTTGTGGTGTTCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCCGAATGGCTTCTATAGTGCCCGAGGATTACGCAGCTGTCAGGCCACCGGGCCACGCGGATAGCGCTTCAGCCGCTCGGCCAGTTCGCGCCCG  >  NC_002947/403127‑403423

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: