Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 648,290 A→C 7.6% 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 7.6% 750.1 / 17.8 233intergenic (+56/+35)acoR/accCacetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (90/121);  new base C (18/4);  total (108/125)
Fisher's exact test for biased strand distribution p-value = 5.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

CTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCG  >  NC_002947/648143‑648439
                                                                                                                                                    |                                                                                                                                                     
cTGATCCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTAACGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGACCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGcaca                                                                                                                                                     <  1:1538312/150‑1 (MQ=255)
   aggCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAggg                                                                                                                                                  <  1:5166959/148‑1 (MQ=255)
      cGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAACACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCg                                                                                                                                               <  1:6605412/150‑1 (MQ=255)
        gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                                                             <  2:1662835/150‑1 (MQ=255)
        gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                                                             <  2:582095/150‑1 (MQ=255)
         cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCgggg                                                                                                                                            <  2:5224096/150‑1 (MQ=255)
            gggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCCACCCGGCCCCCCAATCTGGGAGGGCAAAAAAAAAAACCGCCCACGGCCgggggtg                                                                                                                                         >  1:993109/1‑147 (MQ=255)
             ggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCt                                                                                                                                        <  2:1530981/150‑1 (MQ=255)
              gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                       <  1:336505/150‑1 (MQ=255)
              gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                       <  2:3784215/150‑1 (MQ=255)
               cTTTGGGCTTTTCGCGGGCGACGCTTTTTCGGCGAGTGCATGAGCACGGAATCGAGATGCCGCGGATGAAGGGGTATAGCCGTGGGGCCGGGAGGCGGCCCCTGATTCTTGAAGGCAGTAAAAAACCCCGGCACAGGGCCGGGGTTCttg                                                                                                                                      <  2:2501019/150‑3 (MQ=255)
               cTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCCACCCGGCCCCCCATTGTTGCAGGCCATAAAAAAACCCCGCCCCCGGCCgggggtgggg                                                                                                                                      >  1:535745/1‑144 (MQ=255)
                tttGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                     <  1:4900728/150‑1 (MQ=255)
                  tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTATGGCAGGCAATAAAAAAACCCCGCACACGGCCGGGGTTGTGGTAc                                                                                                                                   >  2:4597310/1‑150 (MQ=255)
                    ggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGGGTCGGGGGCCCAGATTGTTGGAGGGAATAAAATTAACCGGCACAGGGGCGGGGTTGTTGTTCAt                                                                                                                                 >  2:5462508/1‑150 (MQ=255)
                     gCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGGCCCGGCCCCAGATTCAGGCAGGGCATAAAAAAACCCCGCACACGGCCGGGGTTCGGGTTCATc                                                                                                                                >  1:6128406/1‑150 (MQ=255)
                        tttCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                                                             <  2:4530485/150‑1 (MQ=255)
                        tttCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                                                             <  1:7036726/150‑1 (MQ=255)
                             cgGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACGTCACACTaca                                                                                                                        <  1:1701617/150‑1 (MQ=255)
                             cgGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTaca                                                                                                                        <  1:2272876/150‑1 (MQ=255)
                                gCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACCCGGCCCCAGATTCTGGCAGGCCAAAAAAAAAACCGGCCAAGGGCCCGGGGTCTGGTACATCACACAACACGa                                                                                                                     >  2:6070501/1‑150 (MQ=255)
                                  gACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGACGCGGCCCCACATTATGGCAGGCAATAAAAAAAACCCGCACACGGCCGGGGGTGTGGTACATCACACAACACAAta                                                                                                                   >  1:5989417/1‑149 (MQ=255)
                                    cGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCAc                                                                                                                 <  2:738290/150‑1 (MQ=255)
                                      cTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAAACCCCGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTc                                                                                                               >  2:7436971/1‑150 (MQ=255)
                                        ttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTGTGGCAGGCAATAAAAAAACCCCGCACAGGGCCGGGGTTCTGGTACAAAACACAACACGATCACACCt                                                                                                             >  2:3672486/1‑150 (MQ=255)
                                        ttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCt                                                                                                             <  1:4602279/150‑1 (MQ=255)
                                         tATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTg                                                                                                            <  2:6528272/150‑1 (MQ=255)
                                            cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGtt                                                                                                         <  2:4429590/150‑1 (MQ=255)
                                            cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGACACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGtt                                                                                                         <  1:1789186/150‑1 (MQ=255)
                                             ggCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGATAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                        <  2:6112477/150‑1 (MQ=255)
                                              gCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCCGCCCCAGATTCTGGCAGGCCAAAAAAAACCCCCGCACAGGGCCGGGGGGCGGGTACATCACACAACACCCGCACGCCGGGTTgg                                                                                                       >  2:1804296/1‑150 (MQ=255)
                                                 aGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCa                                                                                                    <  1:6838531/150‑1 (MQ=255)
                                                 aGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCCGCCCGGCCCCCCAATCTGGCCGGCCAAAAAAAAAACCCCCCCACGGCCGGGGGGGGGGGACAACACACAAAACCATCACACCCGGGTGGccc                                                                                                    >  2:3416344/1‑149 (MQ=255)
                                                   tGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGt                                                                                                  <  2:3476317/150‑1 (MQ=255)
                                                             ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                        <  1:7160228/150‑1 (MQ=255)
                                                             ccGCATCGACAGGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                        <  2:7117887/150‑1 (MQ=255)
                                                               gCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                      <  1:4759626/150‑1 (MQ=255)
                                                                  tCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAg                                                                                   <  1:1739553/150‑1 (MQ=255)
                                                                       aTGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAAt                                                                              <  1:2554304/150‑1 (MQ=255)
                                                                        tGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACCCGGCCCCCGATTCTGGGAGGCCATAAAAAAACCCCGCACAGGGCCGGGGGGCTGGTACAACACACAACACGATCCCGCCCGGGTGGCCCGGTTGTGCGCGAGGTAGTGAATg                                                                             >  2:4859363/1‑150 (MQ=255)
                                                                         gCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAAACCCGGCACAGGGCCGGGGTTCTGGTAAATCACACTACACGATCACTCCGGGGTGGGCCGGTTGTGGTCGGGGTTGTGAATGt                                                                            >  2:6635020/1‑150 (MQ=255)
                                                                          ccGCGGATGCAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAAGAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTCCACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                           <  2:5666409/150‑1 (MQ=255)
                                                                              ggATGACGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcg                                                                       <  1:1492726/150‑1 (MQ=255)
                                                                                aTGAAGGGGTAGAGCCCTGGGGCCGCGCCGCCGCCCCACATTCTGGCCGGCCATAAAAAAACCCCGCCCAGGGCCGGGGGTGGGGTAAATCACACGACACCATCACACCGGGTTGGCCCCGGTGTGCGCGAGGGAGTGACTGGTGAcgcc                                                                     >  2:5345619/1‑150 (MQ=255)
                                                                                  gAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccgc                                                                   <  2:3927686/150‑1 (MQ=255)
                                                                                       ggTAGCGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                              <  1:2544510/150‑1 (MQ=255)
                                                                                            agCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                         <  2:7088306/150‑1 (MQ=255)
                                                                                            agCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                         <  2:1437584/150‑1 (MQ=255)
                                                                                              cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                       <  2:7341266/150‑1 (MQ=255)
                                                                                              cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                       >  1:5642814/1‑150 (MQ=255)
                                                                                               tcTGGGGCCGCGTCGCGGCCCCCGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGcc                                                      >  1:3789119/2‑150 (MQ=255)
                                                                                                cTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGGGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCt                                                     >  1:777042/1‑150 (MQ=255)
                                                                                                 tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGa                                                                                                                     >  1:12878/1‑85 (MQ=255)
                                                                                                 tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGa                                                                                                                     <  2:12878/85‑1 (MQ=255)
                                                                                                 tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTGTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                    >  2:2899030/1‑150 (MQ=255)
                                                                                                  ggggCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                   >  1:205941/1‑150 (MQ=255)
                                                                                                  ggggCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                   >  1:4395065/1‑150 (MQ=255)
                                                                                                  ggggCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                   <  1:3196501/150‑1 (MQ=255)
                                                                                                   gggCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATc                                                                                                                   >  1:1068513/1‑85 (MQ=255)
                                                                                                   gggCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtca                                                  >  1:4239233/1‑150 (MQ=255)
                                                                                                   gggCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGGTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtca                                                  >  2:6272902/1‑150 (MQ=255)
                                                                                                       cgcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                      >  2:3496983/1‑110 (MQ=255)
                                                                                                       cgcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCAc                                              >  2:4834731/1‑150 (MQ=255)
                                                                                                          gtcgcgGCCCCAGGTTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCa                                           >  1:2361045/3‑150 (MQ=255)
                                                                                                          gACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGGTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCa                                           >  2:1387915/1‑150 (MQ=255)
                                                                                                           aCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAc                                          >  1:1948079/1‑150 (MQ=255)
                                                                                                             gcgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCCCTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCa                                        <  2:4792277/150‑1 (MQ=255)
                                                                                                              cgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAg                                       >  1:2133869/1‑150 (MQ=255)
                                                                                                               ggCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAgg                                      >  2:1090042/1‑150 (MQ=255)
                                                                                                                 ccccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTc                                    >  2:2002082/1‑150 (MQ=255)
                                                                                                                 ccccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTc                                    <  2:4179/150‑1 (MQ=255)
                                                                                                                  cccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcg                                   >  2:670000/1‑150 (MQ=255)
                                                                                                                   ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGt                                                                                  >  2:904228/1‑102 (MQ=255)
                                                                                                                   ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                  >  1:359470/1‑150 (MQ=255)
                                                                                                                   ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                                  >  2:6008084/1‑150 (MQ=255)
                                                                                                                     aGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCgg                                >  2:3641892/1‑150 (MQ=255)
                                                                                                                         tCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCa                            <  2:6333621/150‑1 (MQ=255)
                                                                                                                            ggcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCg                         >  1:2876897/1‑150 (MQ=255)
                                                                                                                               aggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGa                      >  2:6886754/1‑150 (MQ=255)
                                                                                                                                 ggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                     <  2:5183063/150‑1 (MQ=255)
                                                                                                                                  gcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATg                    <  2:6575041/150‑1 (MQ=255)
                                                                                                                                   caATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                   <  1:6790801/150‑1 (MQ=255)
                                                                                                                                    aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGtt                  >  1:5723675/1‑150 (MQ=255)
                                                                                                                                    aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGtt                  >  1:1339922/1‑150 (MQ=255)
                                                                                                                                    aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGtt                  >  2:646481/1‑150 (MQ=255)
                                                                                                                                     aTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg                 <  1:2812442/150‑1 (MQ=255)
                                                                                                                                     aTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg                 <  1:2811359/150‑1 (MQ=255)
                                                                                                                                      tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAGGTTgg                >  1:834150/1‑150 (MQ=255)
                                                                                                                                       ccagaaCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGCGCTCGAGGTAGTGAATGTTGa                                                                         <  2:1915963/88‑1 (MQ=38)
                                                                                                                                         aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt             <  1:2489555/150‑1 (MQ=255)
                                                                                                                                           aaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTg           >  1:6359286/1‑150 (MQ=255)
                                                                                                                                           aaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTg           >  1:2611295/1‑150 (MQ=255)
                                                                                                                                            aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa          <  2:2817382/150‑1 (MQ=255)
                                                                                                                                            aCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa          <  1:1912342/150‑1 (MQ=255)
                                                                                                                                               ccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                                 >  2:2341978/1‑124 (MQ=255)
                                                                                                                                                cGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGcc      >  2:2108753/1‑150 (MQ=255)
                                                                                                                                                 ggCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                                                            >  2:5426734/1‑95 (MQ=255)
                                                                                                                                                 ggCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                                                            <  1:5426734/95‑1 (MQ=255)
                                                                                                                                                    acaGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTcg  >  2:4379817/1‑150 (MQ=255)
                                                                                                                                                    acaGGGCCGGGGTTCTGGTACATCACAATACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGTGGTAGTGAATGTTGAAGCCTCCTTTGCAGAAGCCTTCATCACGCACCAGGTATCGGTGCAGCGGGATGTTGTTCTTGATGCCGTcg  <  2:5591998/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCG  >  NC_002947/648143‑648439

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: