Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 648,290 A→C 9.3% intergenic (+56/+35) acoR → / ← accC acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947648,2900AC 9.3% 418.2 / 19.1 143intergenic (+56/+35)acoR/accCacetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (60/67);  new base C (13/2);  total (73/70)
Fisher's exact test for biased strand distribution p-value = 5.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGT  >  NC_002947/648147‑648437
                                                                                                                                               |                                                                                                                                                   
 gCGGCTGGGGCTTTGGGCTTTTCGCGGGCTCCGGTTTTTGGGCGTGTGCATGAGCACCGCTTCGACATGCCGCGGATGATGGGGTAGAGCCGTGGGGCCGCGACGGGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGcc                                                                                                                                              <  1:2790310/150‑1 (MQ=255)
   ggCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCGGGCAGGCAATAAAAAACCCCGGCACAGGGCCgg                                                                                                                                            <  2:5083263/150‑1 (MQ=255)
    gCACGGGCTTTGGGCATTTTGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCCCATTCTTGCAGGCAAAAAAAAACCCCCGCCCACGGCCggg                                                                                                                                           >  2:5262284/1‑150 (MQ=255)
    gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCCCATAATGGCCCGCCATAAAAAACACCCCCACACGGCCggg                                                                                                                                           >  2:3796355/1‑150 (MQ=255)
    gCACGGGCGTTGGGCATTTTGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCACATTAAGGCAGGCAATAAAAAAACCCCGCACAGGGCCggg                                                                                                                                           >  2:5262283/1‑150 (MQ=255)
     cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGGCGCCGCCCCCCCCTCAGGGAGGCCAAAAAAAAAACCCCCCCACGGCCgggg                                                                                                                                          >  1:1510673/1‑150 (MQ=255)
     cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCgggg                                                                                                                                          <  2:435372/150‑1 (MQ=255)
        gggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGttt                                                                                                                                       <  2:596541/150‑2 (MQ=255)
            tttGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCGCCCCGGCCCCCCAATCAGGCAGGCAAAAAAAAAAACCCCCACAAGGCCGGGGGTCTGGt                                                                                                                                   >  2:4110382/1‑150 (MQ=255)
            tttGGGCATTTCGCGGGCGACGCTTGATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                   <  1:1621229/150‑1 (MQ=255)
                ggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGACGCCGCCCCACATTCTGGCAGGCCAAAAAAAAACCCCGCCCCCGGCCGGGGGTCTGGGACaa                                                                                                                               >  2:298232/1‑149 (MQ=255)
                ggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGCTTCTGGCCGGCAAAAAAAAAACCCCGCCCACGGCCGGGGGTAGGGTACAt                                                                                                                               >  2:3028655/1‑150 (MQ=255)
                 gCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGACCCCAGATTCTGGCAGGCAATAAAAAACCCCCGCACAGGGCCGGGGTTCTGGTACATc                                                                                                                              >  1:3033615/1‑150 (MQ=255)
                 gCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGACCCCAGATTCTGGCAGGCAATAAAAAAACCCCGCACAGGGCCGGGGTTCTGGTACATc                                                                                                                              >  1:3034327/1‑150 (MQ=255)
                                  cTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTc                                                                                                             <  2:147858/150‑1 (MQ=255)
                                        cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGtt                                                                                                       <  1:2133540/150‑1 (MQ=255)
                                         ggCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                      <  2:3752915/150‑1 (MQ=255)
                                                 cATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                                                              <  1:2354176/150‑1 (MQ=255)
                                                  aTGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCGGCGCGGCCCCAGATTCTGGCCGGCAATAAAAAACCCCCGGACAGGGCCGGGGGTCGGGTACATCACACGACACGCTCACGCCGGGGTGGCCCGtgtg                                                                                             >  2:4015140/1‑150 (MQ=255)
                                                  aTGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCCGATTCTGGCAGGCAAAAAAAAAACCCCGCACCGGGCCCGGGGTCGGGTACATCACACGACACCATCACTCCAGGTTGGCCAGTTtg                                                                                             >  2:4015192/1‑150 (MQ=255)
                                                   tGAGCACCGCATCGACATGCCGCGGATGGAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTtgt                                                                                            <  1:287447/150‑1 (MQ=255)
                                                     aGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGc                                                                                          <  2:4489774/150‑1 (MQ=255)
                                                     aGCACCGCATCGACAGGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGc                                                                                          <  2:4343261/150‑1 (MQ=255)
                                                         ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                      <  1:2569280/150‑1 (MQ=255)
                                                         ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                      <  2:2296509/150‑1 (MQ=255)
                                                         ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                      <  2:163579/150‑1 (MQ=255)
                                                           gCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCGACGCGGCCCCAGATTATTGCAGGCAAAAAAAAAAACCCGCACAGGGGCGGGGTTATGGTACATAACACAAAAACATAAAAACAGGTTTGCCCGTGTGTGCTCGAgg                                                                                    >  2:3888665/1‑150 (MQ=255)
                                                                gACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCCCATTCTGGCCGGCCAAAAAAAACCCCCGCCCACGGCCGGGGGTCGGGTACATCACACGACACGATCACTACGGGGTGGGCCGGTTGTGCGCGAGGTAGTg                                                                               >  2:5307179/1‑150 (MQ=255)
                                                                 aCATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCCCCGCCCCAGAGTCTGGCCGGCAAAAAAAAAACCCCGCACACGGCCCGGGGTGGGGGACATCACACAACAAGAGCACACCGGGGTGGCCCGGGTGTGCTCGTGGTAGTgg                                                                              >  2:1976825/1‑149 (MQ=255)
                                                                    tGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCCCGCCGCGGCCCCAGCTTCTGGCCGGCCAAAAAAAAACCCCGCCCAGGGCCGGGGGTCGGGTACAACACACAACAAGAGCACTCCGGGGTGGGCCGGTTGTGGTCGCGGTGGTGTATg                                                                           >  2:5286626/1‑150 (MQ=255)
                                                                        gcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCGCAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGa                                                                       <  1:2498473/150‑1 (MQ=255)
                                                                         cgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                      <  2:1184053/150‑1 (MQ=255)
                                                                         atcttGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                      <  1:543584/146‑1 (MQ=255)
                                                                              gAAGGGGTAGAGCCCTGGGGGCGCGACGCGGCCCCCGGTTCTGGCCGGCAATAAAAAAAACCCGCCCCGGGCCGGGGGTGTGGTAAATCACACAACACGATCACACCGGGTTTGGCCGGGTGTGGTCGAGGTTGGGAAGGTTGAcgccgc                                                                 >  2:799671/1‑150 (MQ=255)
                                                                                 ggggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCtt                                                              <  1:1648684/150‑1 (MQ=255)
                                                                                  gggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCttt                                                             <  2:667747/150‑1 (MQ=255)
                                                                                  gggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCttt                                                             <  2:2474146/150‑1 (MQ=255)
                                                                                   ggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                            <  2:3555494/150‑1 (MQ=255)
                                                                                   ggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                            <  1:672718/150‑1 (MQ=255)
                                                                                        agCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                       <  2:2890716/150‑1 (MQ=255)
                                                                                        agCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                       <  2:953470/150‑1 (MQ=255)
                                                                                         gCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGGTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAg                                                      >  2:5334465/1‑150 (MQ=255)
                                                                                         gCCCTGGCGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCCGGGATCTGGAACATCACACTACAAGATCACACCTGCTTGGCCACTTAGTGCTCGACGTAGTGAAACTTGACGCCGCCTTTGCAGAAg                                                      >  2:4734720/1‑150 (MQ=255)
                                                                                          cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                     >  1:1225682/1‑150 (MQ=255)
                                                                                           tcTGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGcc                                                    >  2:2185972/2‑150 (MQ=255)
                                                                                             tGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                  <  1:3952763/150‑1 (MQ=255)
                                                                                              ggggCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                 >  1:4413770/1‑150 (MQ=255)
                                                                                              ggggCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                 >  1:3275241/1‑150 (MQ=255)
                                                                                              ggggCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                 >  1:4413690/1‑150 (MQ=255)
                                                                                                      gACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCa                                         >  1:1061150/1‑150 (MQ=255)
                                                                                                         gcgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCa                                      >  1:2554733/1‑150 (MQ=255)
                                                                                                            gCCCCAGATTCTGGCAGGCAATAAAAATCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGt                                   >  2:4929899/1‑150 (MQ=255)
                                                                                                            gCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGGGCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGt                                   >  2:5114094/1‑150 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                      >  1:3480156/1‑80 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                      <  2:3480156/80‑1 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                >  2:1776710/1‑150 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                >  2:2867781/1‑150 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                >  1:2348374/1‑150 (MQ=255)
                                                                                                               ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                >  1:3366470/1‑150 (MQ=255)
                                                                                                               ccAGATTCTGACAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                >  2:3711268/1‑150 (MQ=255)
                                                                                                                 aGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCgg                              <  2:2259782/150‑1 (MQ=255)
                                                                                                                     tCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCa                          >  2:1432634/1‑150 (MQ=255)
                                                                                                                       tGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGATCGAGGTAGTGACTGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                        >  1:5203972/1‑150 (MQ=255)
                                                                                                                         gcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCgg                      >  2:1067361/1‑150 (MQ=255)
                                                                                                                         gcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCgg                      <  1:1432634/150‑1 (MQ=255)
                                                                                                                         gcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACGACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCgg                      <  1:2369506/150‑1 (MQ=255)
                                                                                                                           aggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGa                    <  1:5399634/150‑1 (MQ=255)
                                                                                                                            ggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                   >  2:171114/1‑150 (MQ=255)
                                                                                                                            ggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                   <  1:665303/150‑1 (MQ=255)
                                                                                                                             gcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATg                  <  2:461915/150‑1 (MQ=255)
                                                                                                                              caATAAAAAACCCCGGCACAGGTCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                 >  1:2455640/1‑150 (MQ=255)
                                                                                                                              caATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                 <  1:1067361/150‑1 (MQ=255)
                                                                                                                               aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGtt                >  1:3396098/1‑150 (MQ=255)
                                                                                                                                aTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg               >  1:1810917/1‑150 (MQ=255)
                                                                                                                                aTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg               <  1:559683/150‑1 (MQ=255)
                                                                                                                                 tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCCGTTTGTGCTCGAGGTAGTGACTGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTgg              >  2:4659097/1‑150 (MQ=255)
                                                                                                                                 tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTgg              >  2:3646766/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt             <  2:2477293/150‑1 (MQ=255)
                                                                                                                                   cagaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                                                          >  1:2572496/4‑104 (MQ=38)
                                                                                                                                   cagaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc            >  1:1725728/4‑150 (MQ=255)
                                                                                                                                   cagaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCggg                     <  2:278375/138‑1 (MQ=255)
                                                                                                                                   aaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc            >  2:4867457/1‑150 (MQ=255)
                                                                                                                                   aaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc            >  1:4768438/1‑150 (MQ=255)
                                                                                                                                   aaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc            >  1:5277313/1‑150 (MQ=255)
                                                                                                                                    agaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt           >  2:3708401/3‑150 (MQ=255)
                                                                                                                                    aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                 >  2:1688589/1‑144 (MQ=255)
                                                                                                                                    aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                 <  1:1688589/144‑1 (MQ=255)
                                                                                                                                     gaaCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt          >  1:1906072/2‑150 (MQ=255)
                                                                                                                                     aaaCCCCGGCACAGGGCCGGGGTTCTGGTCCATCACCCTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCGTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt          <  2:4073624/150‑1 (MQ=255)
                                                                                                                                     aaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt          >  2:4869908/1‑150 (MQ=255)
                                                                                                                                      aaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTg         <  2:3408524/150‑1 (MQ=255)
                                                                                                                                       aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa        <  2:2485909/150‑1 (MQ=255)
                                                                                                                                       aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa        <  1:41450/150‑1 (MQ=255)
                                                                                                                                       aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa        <  1:2957132/150‑1 (MQ=255)
                                                                                                                                       aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa        <  2:2904191/150‑1 (MQ=255)
                                                                                                                                         cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg      <  1:32055/150‑1 (MQ=255)
                                                                                                                                         cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg      >  1:2332127/1‑150 (MQ=255)
                                                                                                                                             gCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt  <  1:799671/150‑1 (MQ=255)
                                                                                                                                             gCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt  <  1:4944038/150‑1 (MQ=255)
                                                                                                                                             gCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt  <  2:3410061/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGT  >  NC_002947/648147‑648437

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: