Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,010,005 G→C 15.4% intergenic (+38/‑138) PP_5259 → / → ydcJ LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,010,0050GC15.4% 831.7 / 127.5 298intergenic (+38/‑138)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (103/149);  new base C (16/30);  total (119/179)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 315 total aligned reads displayed.

CATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAAC  >  NC_002947/6009857‑6010154
                                                                                                                                                    |                                                                                                                                                     
cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                      >  1:6010218/1‑150 (MQ=255)
 aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                                                     >  1:4939977/1‑150 (MQ=255)
 aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgcg                                                                                                                                                     >  1:5322328/1‑150 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGt                                                                                                                                                    <  2:3424397/150‑1 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCGACGGGGACCGCGt                                                                                                                                                    >  1:6319782/1‑150 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGt                                                                                                                                                    <  2:5972547/150‑1 (MQ=255)
  tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGt                                                                                                                                                    <  1:4193085/150‑1 (MQ=255)
   gCAATGCGCAGTGGCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                   <  2:6614480/150‑1 (MQ=255)
   gCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                   <  2:2593722/150‑1 (MQ=255)
   gCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                   <  2:1839435/150‑1 (MQ=255)
   gCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTa                                                                                                                                                   <  1:3754925/150‑1 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAg                                                                                                                                                  <  2:1186803/150‑1 (MQ=255)
    cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAg                                                                                                                                                  <  1:2640941/150‑1 (MQ=38)
        gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                              <  2:1607232/150‑1 (MQ=255)
        gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                              <  2:1606972/150‑1 (MQ=255)
        gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                              <  1:176302/150‑1 (MQ=255)
        gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                              <  2:3986975/150‑1 (MQ=255)
        gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAGgagc                                                                                                                                              <  1:7276894/150‑1 (MQ=255)
         cgcAGTGCCGGAGCGCATCACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGagcg                                                                                                                                             <  1:1851520/150‑1 (MQ=37)
          gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCCCTCCTACGGGGACCGCGTAGGCGCgg                                                                                                                                            >  1:7067384/1‑150 (MQ=255)
           cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                           <  1:1367178/150‑1 (MQ=255)
               gCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                       >  2:199707/1‑150 (MQ=255)
               gCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                       <  1:1258600/150‑1 (MQ=255)
               gCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                       <  1:2805583/150‑1 (MQ=255)
                 cGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCGTTAGc                                                                                                                                     >  2:3193461/1‑150 (MQ=255)
                  ggAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTACGACGGGGACCGCGTAGGCGCGGCTTTAGcc                                                                                                                                    >  1:6381016/1‑150 (MQ=255)
                   gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                                   >  2:6852741/1‑150 (MQ=255)
                   gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                                   <  2:7000250/150‑1 (MQ=255)
                    aGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                                  <  2:4406416/150‑1 (MQ=255)
                       gcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaa                                                                                                                               >  1:6101292/1‑150 (MQ=255)
                            aCCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGc                                                                                                                          >  2:259354/1‑150 (MQ=255)
                             ccTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGcc                                                                                                                         <  2:2342980/150‑1 (MQ=255)
                                  cGTATGCCGAGCATGCGGCAGAGGTGCCATAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGa                                                                                                                    <  2:4651238/150‑1 (MQ=255)
                                    tttGCCGTGCATGCGGCAGAGGTGCCAACGGTGCGGGGGCTGGTGGATTGGTTTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTcaaac                                                                                                                  <  1:7300501/148‑5 (MQ=14)
                                       gCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCg                                                                                                               >  1:1664311/1‑150 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                             <  1:5504002/150‑1 (MQ=255)
                                         cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGt                                                                                                                                                    >  1:913155/1‑111 (MQ=255)
                                          gAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCGAAATCCGCACCGACGGGGCCCGCGGTGGCGCGGCGTTAGCCCGGACAGGGGCGGCGGGGGCGaga                                                                                                            >  1:2473601/1‑148 (MQ=255)
                                               tGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGCCGATc                                                                                                            <  1:3468633/145‑1 (MQ=255)
                                                 cGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                                                     <  1:6534051/150‑1 (MQ=255)
                                                    cAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                                  <  1:2585467/150‑1 (MQ=255)
                                                       agGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCTGGTTTGACGGGCTCTTCGCGGCTACAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCAGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGa                                                                                               <  1:4801266/150‑1 (MQ=255)
                                                              aaaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAATCCGCTACAAACGGGACCCCGTAGGCGCGGCTTTAGACCGGAAACGGCCGTCGAGGTCGATCAAAAAAGGAATGacaaacac                                                                                        >  1:6916948/1‑148 (MQ=255)
                                                                aGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTGGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcactcac                                                                                      <  2:2149681/150‑1 (MQ=255)
                                                                    gCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtttt                                                                                  <  2:913348/150‑1 (MQ=255)
                                                                       ggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttttt                                                                               <  2:317178/150‑1 (MQ=255)
                                                                          gcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGc                                                                            <  2:7173162/150‑1 (MQ=255)
                                                                            tggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                          <  1:7426846/150‑1 (MQ=255)
                                                                            gggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTACAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                          <  2:4553085/149‑1 (MQ=255)
                                                                             ggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                         <  2:6276176/150‑1 (MQ=255)
                                                                              gtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                        <  1:1009077/150‑1 (MQ=255)
                                                                              gtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTATTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGTGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                        <  1:6650165/150‑1 (MQ=255)
                                                                              gtggAAGGGATGGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                        <  2:4928400/150‑1 (MQ=255)
                                                                               tggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgtg                                                                       <  2:3713596/150‑1 (MQ=255)
                                                                                ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCGCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                      <  2:5892389/150‑1 (MQ=255)
                                                                                 gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                                                     >  2:6336378/1‑150 (MQ=255)
                                                                                 gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACACACTTTTTTTCGCTTGTGAg                                                                     >  2:653766/1‑150 (MQ=255)
                                                                                   aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCa                                                                   >  2:4078494/1‑150 (MQ=255)
                                                                                    tGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCaa                                                                  <  2:6134773/150‑1 (MQ=255)
                                                                                        ttGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTg                                                              <  2:5760579/150‑1 (MQ=255)
                                                                                         tGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGGGAGCAAAGTGc                                                             >  2:5649804/1‑150 (MQ=255)
                                                                                          ggCGAGCAGTTGCTGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGCTCAAAAAAGGAATGACAGACTCACTTTTTTTCGCTTGTGAGCAATGTGCa                                                            <  2:6169436/150‑1 (MQ=255)
                                                                                          ggCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCa                                                            <  1:3147275/150‑1 (MQ=255)
                                                                                             gAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGGGGCCCGCCTAGGCGCGGCGTTAGCCGGGGAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACAAACATTTTTTCGCTTGTGAGCAAAGTGCAGTc                                                         >  1:4663812/1‑150 (MQ=255)
                                                                                                  gTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGc                                                    <  2:7349200/150‑1 (MQ=255)
                                                                                                   gtGCAGGTTTGCCGGCCTCTTCGCGGCTACAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCt                                                   <  2:2427419/149‑1 (MQ=255)
                                                                                                           tcgACGGCCCTTTCCCGGCTAAAGCCCCGCCTAAGCGGTCCCCGTATGAGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAAACACACATTTTTTCGCGTGGGAGCAAAGTGCAATCCCAGCATTCATGTa                                           >  2:1307203/3‑150 (MQ=16)
                                                                                                              aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGc                                        >  2:5070551/1‑150 (MQ=255)
                                                                                                                     ctTCGCGGCTAAAGCCGCTCCTCCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGaa                                 <  2:5341830/150‑1 (MQ=255)
                                                                                                                     ctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGaa                                 <  2:4796576/150‑1 (MQ=255)
                                                                                                                      ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                                >  2:2837602/1‑150 (MQ=255)
                                                                                                                      ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                                <  1:199707/150‑1 (MQ=255)
                                                                                                                      ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                                <  1:2181680/150‑1 (MQ=255)
                                                                                                                         gcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCa                             >  1:5002142/1‑150 (MQ=255)
                                                                                                                           ggCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                           <  2:3554314/150‑1 (MQ=255)
                                                                                                                              tAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCTAGCTTTCATGTAAGCGTCCGAACCCAACCAg                        >  1:1518045/1‑150 (MQ=255)
                                                                                                                               aaaGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAgg                       <  2:5663967/150‑1 (MQ=255)
                                                                                                                                aaGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                         >  1:3412776/1‑99 (MQ=25)
                                                                                                                                aaGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtg                                                                         <  2:3412776/99‑1 (MQ=25)
                                                                                                                                 aGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                     <  2:4357116/150‑1 (MQ=255)
                                                                                                                                  gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg                    >  2:6595311/1‑150 (MQ=255)
                                                                                                                                  gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg                    >  1:5241178/1‑150 (MQ=255)
                                                                                                                                  gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGCGAAGAGGCCGTCAAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg                    >  2:4481339/1‑150 (MQ=255)
                                                                                                                                   ccGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtttttt                                                                                <  2:5032074/89‑1 (MQ=21)
                                                                                                                                     gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCg                                                                             >  2:7319049/4‑90 (MQ=25)
                                                                                                                                     gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                 >  2:1091902/4‑150 (MQ=39)
                                                                                                                                     gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                 >  1:7026106/4‑150 (MQ=37)
                                                                                                                                     gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                 <  2:1494050/150‑1 (MQ=255)
                                                                                                                                     gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa                 <  2:7067292/150‑1 (MQ=255)
                                                                                                                                      cgccTACGCGGTCCCCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAg                >  1:4044540/3‑150 (MQ=255)
                                                                                                                                      cTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAg                <  1:3330588/150‑1 (MQ=255)
                                                                                                                                       tCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGc               >  1:3844521/1‑150 (MQ=255)
                                                                                                                                        ccTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCg              >  1:6777531/1‑150 (MQ=255)
                                                                                                                                          tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGAt            >  2:815175/1‑150 (MQ=255)
                                                                                                                                          tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTAGGGAGAAAAGTGCATTCACAGCTTTCATGTAAGCGTCCGAACCCAAACAGGAAGCTAGCGAt            >  2:2933060/1‑150 (MQ=255)
                                                                                                                                            cGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGc          <  2:320653/150‑1 (MQ=255)
                                                                                                                                              gggACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGccc        <  1:6712889/150‑1 (MQ=255)
                                                                                                                                                   cgcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCaa   >  2:2788030/1‑150 (MQ=255)
                                                                                                                                                    gcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAAc  >  2:3841836/1‑150 (MQ=255)
                                                                                                                                                    gcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAAc  >  1:5724094/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAAC  >  NC_002947/6009857‑6010154

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: