Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 648,250 A→T 7.9% intergenic (+16/+75) 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,2500AT 7.9% 299.4 / 14.8 102intergenic (+16/+75)acoR/accCacetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (41/53);  new base T (5/3);  total (46/56)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

Only 100 of 104 total aligned reads displayed.

TGTTGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAC  >  NC_002947/648104‑648399
                                                                                                                                                  |                                                                                                                                                     
tGTTGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGAcgc                                                                                                                                                    <  2:560936/150‑1 (MQ=255)
    gCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGcc                                                                                                                                                >  2:1697001/1‑150 (MQ=255)
         aCGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGTCGCGGCCCCAgg                                                                                                                                           <  1:345834/150‑2 (MQ=255)
         aCGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGa                                                                                                                                           >  2:532687/1‑150 (MQ=255)
            cgcTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTc                                                                                                                                        >  1:243982/1‑150 (MQ=255)
             gcTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCt                                                                                                                                       >  1:1043567/1‑150 (MQ=255)
              cTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTg                                                                                                                                      <  2:1544380/150‑1 (MQ=255)
              cTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTg                                                                                                                                      >  2:451378/1‑150 (MQ=255)
               tGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTgg                                                                                                                                     >  1:1829217/1‑150 (MQ=255)
                  tCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggcag                                                                                                                                  >  1:1002180/1‑150 (MQ=255)
                   ccGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggcagg                                                                                                                                 <  2:243982/150‑1 (MQ=255)
                   ccGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggcagg                                                                                                                                 <  1:2125947/150‑1 (MQ=255)
                    cGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggcaggc                                                                                                                                >  1:1445685/1‑150 (MQ=255)
                     gCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggcaggca                                                                                                                               >  2:1321162/1‑150 (MQ=255)
                      ccATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAATGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCaa                                                                                                                              >  1:7677/1‑150 (MQ=255)
                        aTCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGTGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATa                                                                                                                            >  1:466587/1‑150 (MQ=255)
                        aTCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGTCGCGGCCCCa                                                                                                                                             >  2:2022165/1‑133 (MQ=255)
                             tGGAAGCCTGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGCAGGGGTAGAGCCCTGGGGCCGCGACGCGGCACCAGATTCTGGCAGGCAAGaaaaaa                                                                                                                       >  2:378377/1‑150 (MQ=255)
                             tGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGGAGGCAATaaaaaa                                                                                                                       >  1:211386/1‑150 (MQ=255)
                              ggAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggc                                                                                                                                    <  1:1572677/136‑1 (MQ=255)
                              ggAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggc                                                                                                                                    >  2:1572677/1‑136 (MQ=255)
                               gAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAATGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAAcc                                                                                                                     <  2:7677/150‑1 (MQ=255)
                                  gCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                                                  <  1:451378/150‑1 (MQ=255)
                                  gCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCACGGATGAAGGGGCAGACCCCTGGGGCAGCGACGCGGCCCCAGAGTCAGGCAGGCAATacaatgggtcg                                                                                                                  <  2:704861/150‑8 (MQ=255)
                                   ccGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggca                                                                                                                                   >  1:1288653/1‑132 (MQ=255)
                                   ccGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTggca                                                                                                                                   <  2:1288653/132‑1 (MQ=255)
                                       ctgctgCGGCACGGGCTTTGGGCATTTCGCGTGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGCTTCTGGCAGGCAATAAAAAACCACGGcaca                                                                                                             <  1:1109451/150‑1 (MQ=255)
                                       cTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCgg                                                                                                                 <  2:779609/146‑1 (MQ=255)
                                       cTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGCAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCgg                                                                                                                 >  1:779609/1‑146 (MQ=255)
                                        tGATGCGGCCCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAg                                                                                                            <  2:1425280/150‑1 (MQ=255)
                                               gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGGGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                     >  2:2131073/1‑150 (MQ=255)
                                               gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                     <  2:1749577/150‑1 (MQ=255)
                                               gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                     <  1:1604780/150‑1 (MQ=255)
                                               gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATAGGAGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                     <  1:1935003/150‑1 (MQ=255)
                                                    ggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCt                                                                                                <  1:104437/150‑1 (MQ=255)
                                                    ggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGAATCGACATGCCGCGGATGAAGGGGTAGAGTCCTAGGGCCGCGACGCGGCCCCAGATTCAAAAAAAAAAAAAAAAACCCCGGCACAGGGCCGGGGTTCt                                                                                                <  2:712037/150‑1 (MQ=255)
                                                     gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                               <  1:2074711/150‑1 (MQ=255)
                                                      cTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                              <  1:1984218/150‑1 (MQ=255)
                                                      cTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                              <  1:593207/150‑1 (MQ=255)
                                                       tttGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                             >  1:90466/1‑150 (MQ=255)
                                                        ttGGGCATTTCGCGGGCGACGCTGTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGttttttat                                                                                            <  2:211386/150‑7 (MQ=255)
                                                         tGGGCATTTCGCGGGCGACGCTTTAGCGGCGAGGGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGCAGACCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGttttttatt                                                                                           <  2:466587/150‑8 (MQ=255)
                                                          gggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                                          <  2:1043567/150‑1 (MQ=255)
                                                            gCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATc                                                                                        <  2:1204333/150‑1 (MQ=255)
                                                                ttCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCacac                                                                                    <  2:599901/150‑1 (MQ=255)
                                                                    cgGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTaca                                                                                <  2:2049526/150‑1 (MQ=255)
                                                                        cGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAAc                                                                                                                      >  2:1544621/1‑108 (MQ=255)
                                                                        cGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAAc                                                                                                                      <  1:1544621/108‑1 (MQ=255)
                                                                             cTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTc                                                                       >  1:1294522/1‑150 (MQ=255)
                                                                              tttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGt                                                                                                   >  2:1406819/1‑121 (MQ=255)
                                                                              tttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGt                                                                                                   <  1:1406819/121‑1 (MQ=255)
                                                                              tttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTcc                                                                      >  1:233821/1‑150 (MQ=255)
                                                                              tttATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTcc                                                                      >  1:2285127/1‑150 (MQ=255)
                                                                                   cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGtt                                                                 <  2:1002180/150‑1 (MQ=255)
                                                                                   cGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGtt                                                                 <  1:1697001/150‑1 (MQ=255)
                                                                                     gCGAGTGCATGCGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGTGGTTCTGGTACATCACACTACACGATCACTCCTGGTTgg                                                               <  1:2323278/150‑1 (MQ=255)
                                                                                            cATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                        >  2:116817/1‑150 (MQ=255)
                                                                                            cATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGACTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                        <  2:2241598/150‑1 (MQ=255)
                                                                                             aTGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTtg                                                       >  2:2153852/1‑150 (MQ=255)
                                                                                                aGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                        >  1:378682/1‑146 (MQ=255)
                                                                                                aGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                        <  2:378682/146‑1 (MQ=255)
                                                                                                 gCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                              >  1:739032/1‑107 (MQ=255)
                                                                                                 gCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                              <  2:739032/107‑1 (MQ=255)
                                                                                                    ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGTCGCGGCCCCa                                                                                                                                             <  2:1690461/57‑1 (MQ=255)
                                                                                                    ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGTCGCGGCCCCa                                                                                                                                             >  1:1690461/1‑57 (MQ=255)
                                                                                                    ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                <  2:1829217/150‑1 (MQ=255)
                                                                                                     cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAg                                               <  1:651613/150‑1 (MQ=255)
                                                                                                            aCATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAgg                                                                                                           >  1:1942447/1‑83 (MQ=255)
                                                                                                            aCATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAgg                                                                                                           <  2:1942447/83‑1 (MQ=255)
                                                                                                               tGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATg                                     <  2:822115/150‑1 (MQ=255)
                                                                                                                gCCGCGGATGAAGGGGTATAGGCCTTGGTCGGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAACACCGGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGt                                    <  1:532687/150‑1 (MQ=255)
                                                                                                                   gcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATaa                                                                                                                           >  2:25867/1‑60 (MQ=255)
                                                                                                                   gcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATaa                                                                                                                           <  1:25867/60‑1 (MQ=255)
                                                                                                                    cgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                <  1:116817/150‑1 (MQ=255)
                                                                                                                         gAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccgc                           <  2:2285127/150‑1 (MQ=255)
                                                                                                                            ggggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                     <  1:1028643/89‑1 (MQ=255)
                                                                                                                            ggggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                     >  2:1028643/1‑89 (MQ=255)
                                                                                                                            ggggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCtt                        <  2:213809/150‑1 (MQ=255)
                                                                                                                              ggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGa                                                                             <  1:1241545/95‑1 (MQ=255)
                                                                                                                              ggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGa                                                                             >  2:1241545/1‑95 (MQ=255)
                                                                                                                                tAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                    <  1:2131073/150‑1 (MQ=255)
                                                                                                                                   agCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                 >  1:587600/1‑150 (MQ=255)
                                                                                                                                     ccttGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc               >  2:1581038/4‑150 (MQ=255)
                                                                                                                                     cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc               <  2:2260193/150‑1 (MQ=255)
                                                                                                                                     cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc               <  1:2082674/150‑1 (MQ=255)
                                                                                                                                       cTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCt             >  1:225311/1‑150 (MQ=255)
                                                                                                                                        tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt            >  2:2176803/1‑150 (MQ=255)
                                                                                                                                        tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt            >  2:2176892/1‑150 (MQ=255)
                                                                                                                                        tGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTGGCAGAAGCCtt            >  1:534971/1‑150 (MQ=255)
                                                                                                                                            gCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTCCACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAg                                               <  1:1170431/111‑1 (MQ=255)
                                                                                                                                            gCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAg                                               >  2:1170431/1‑111 (MQ=255)
                                                                                                                                            gCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcatc        >  1:2014764/1‑150 (MQ=255)
                                                                                                                                              cgcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCAc      >  1:1122414/1‑150 (MQ=255)
                                                                                                                                               gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGAt                                                                            <  1:2288863/79‑1 (MQ=255)
                                                                                                                                               gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGAt                                                                            >  2:2288863/1‑79 (MQ=255)
                                                                                                                                               gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTgg                                                               >  1:2334323/1‑92 (MQ=255)
                                                                                                                                               gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTgg                                                               <  2:2334323/92‑1 (MQ=255)
                                                                                                                                                cgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGc    >  1:1189267/1‑150 (MQ=255)
                                                                                                                                                  aCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAc  <  1:1923876/150‑1 (MQ=255)
                                                                                                                                                  aCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAc  >  2:1923876/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGTTGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAC  >  NC_002947/648104‑648399

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: