Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,231,473 T→G 8.4% intergenic (+234/‑192) ygjK → / → fadH glycoside hydrolase/2,4‑dienoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,231,4730TG 8.4% 318.6 / 19.4 107intergenic (+234/‑192)ygjK/fadHglycoside hydrolase/2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base G (3/6);  total (48/59)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 108 total aligned reads displayed.

TTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGAAC  >  NC_000913/3231324‑3231621
                                                                                                                                                     |                                                                                                                                                    
ttGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCt                                                                                                                                                      <  2:944785/150‑1 (MQ=255)
      ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                                                <  1:407270/150‑1 (MQ=255)
       cGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGGCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGc                                                                                                                                               <  2:291494/150‑1 (MQ=18)
        ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                              >  1:1041955/1‑150 (MQ=255)
        ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                              <  2:214730/150‑1 (MQ=255)
        ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGcc                                                                                                                                              <  1:1409682/150‑1 (MQ=21)
         gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGccc                                                                                                                                             <  1:1033424/150‑1 (MQ=255)
         gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGccc                                                                                                                                             <  1:136472/150‑1 (MQ=255)
          cAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccag                                                                                                                                            <  1:34349/150‑3 (MQ=18)
          cAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTAGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccag                                                                                                                                            <  2:499874/150‑3 (MQ=18)
            aCTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTAc                                                                                                                                          >  2:39893/1‑150 (MQ=255)
               gCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCACCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACAtg                                                                                                                                       >  1:3035/1‑150 (MQ=255)
                cACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACAtgt                                                                                                                                      >  1:298118/1‑150 (MQ=255)
                  ccGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGc                                                                                                                                    <  2:1446776/150‑1 (MQ=255)
                   cgcgCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATTCGGCGTTTACGCCTTATCCGCCCTACATGTGCa                                                                                                                                   <  2:463512/150‑1 (MQ=255)
                       ccACTGGCGGATGCGGTGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTc                                                                                                                               <  2:1364164/150‑1 (MQ=255)
                        cACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTcc                                                                                                                              >  2:116451/1‑150 (MQ=255)
                        cACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTcc                                                                                                                              >  2:784320/1‑150 (MQ=255)
                         aCTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCg                                                                                                                             <  2:1612830/150‑1 (MQ=255)
                           tGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCTTCTCATCCTGCATCTGCACCTCGCCACTTGCGGATGCTGCGTCGACGCCTTATCCGCCCTACATGTGCAATTCCGTa                                                                                                                           <  2:1304757/150‑1 (MQ=17)
                           tGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTa                                                                                                                           <  2:615690/150‑1 (MQ=255)
                           gggCGGATAAGGCGTTCACGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTa                                                                                                                           >  2:772405/2‑150 (MQ=17)
                            ggCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                                                                                                          >  2:14487/1‑150 (MQ=255)
                                gATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTg                                                                                                                      <  1:838260/150‑1 (MQ=255)
                                  tGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGa                                                                                                                    <  1:1086834/150‑1 (MQ=255)
                                   gcggcgTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGAt                                                                                                                   <  1:243912/150‑1 (MQ=255)
                                    cggcgTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATa                                                                                                                  >  2:1181894/1‑150 (MQ=255)
                                        gTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAc                                                                                                              <  2:140301/150‑1 (MQ=255)
                                        gTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAc                                                                                                              >  2:1365269/1‑150 (MQ=255)
                                          gAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcgc                                                                                                            <  1:116451/150‑1 (MQ=255)
                                           aaCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcgcg                                                                                                           <  1:14487/150‑1 (MQ=255)
                                            aCGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGa                                                                                                          >  1:140955/1‑150 (MQ=255)
                                             cGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccagtg                                                                                                                                          <  1:168777/117‑5 (MQ=17)
                                             cGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccagtg                                                                                                                                          >  2:168777/1‑113 (MQ=17)
                                                   aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGt                                                                                                   >  2:1507732/1‑150 (MQ=255)
                                                        cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcg                                                                                                             <  1:288338/135‑1 (MQ=255)
                                                        cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcg                                                                                                             >  2:288338/1‑135 (MQ=255)
                                                           tACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCg                                                                                           <  2:1062259/150‑1 (MQ=255)
                                                            aCATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCgg                                                                                          >  2:720044/1‑150 (MQ=255)
                                                                gtgtgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGtt                                                                                      >  1:1187816/1‑150 (MQ=255)
                                                                 tgtgtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTa                                                                                     <  2:507561/150‑1 (MQ=255)
                                                                   ggtTGCCGTAGGTCGGATTTGCCTCGACTAGCGTCGCATCCTGCATCTGCACGGCGCGTCTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGGGTCGCATCCGGCGTTATc                                                                                   <  2:1381298/149‑1 (MQ=18)
                                                                    gtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCa                                                                                  <  2:714773/150‑1 (MQ=255)
                                                                    gtTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCa                                                                                  >  1:1311567/1‑150 (MQ=255)
                                                                     ttCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATccc                                                                                 >  1:486843/1‑148 (MQ=255)
                                                                       cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                                               >  1:900477/1‑150 (MQ=255)
                                                                       cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                                               >  1:511624/1‑150 (MQ=255)
                                                                         cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                             >  2:523439/1‑150 (MQ=255)
                                                                         cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                             <  2:436936/150‑1 (MQ=255)
                                                                         cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                             <  1:1569617/150‑1 (MQ=255)
                                                                         cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                             <  1:1609887/150‑1 (MQ=255)
                                                                         cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                             >  2:1609887/1‑150 (MQ=255)
                                                                           tAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGACGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCg                                                                           >  1:279424/1‑150 (MQ=255)
                                                                            aGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGt                                                                          >  2:910216/1‑150 (MQ=255)
                                                                             ggTCGGATAAGACGCGACAAGCGTCGCACCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTa                                                                         <  2:3035/150‑1 (MQ=255)
                                                                                 ggATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCt                                                                     >  2:190252/1‑150 (MQ=255)
                                                                                  gATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtt                                                                    >  1:958187/1‑150 (MQ=255)
                                                                                    tAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACTCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtttt                                                                  <  2:769495/150‑1 (MQ=255)
                                                                                    tAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCtttt                                                                  <  2:140955/150‑1 (MQ=255)
                                                                                     aaGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCttttt                                                                 <  1:987959/150‑1 (MQ=255)
                                                                                      aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTg                                                                <  1:784320/150‑1 (MQ=255)
                                                                                      aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTg                                                                <  2:568/150‑1 (MQ=255)
                                                                                            gACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATccc                                                          <  1:772405/150‑1 (MQ=255)
                                                                                             aCAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCa                                                         <  1:439059/150‑1 (MQ=255)
                                                                                             aCAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCa                                                         <  2:383226/150‑1 (MQ=255)
                                                                                              cAAGTGTCTCATGCGGCATCTTCACCGCGCCTCTGGCTGATGCGGCGTGTACTCCTTTTCCGCCCTATATGTGCAATTCCGTATCCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCAt                                                        <  2:440166/150‑1 (MQ=25)
                                                                                               aaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATc                                                       <  2:1311567/150‑1 (MQ=255)
                                                                                                 gCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATcac                                                     <  2:78490/150‑1 (MQ=255)
                                                                                                 gCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATcac                                                     <  1:1352198/150‑1 (MQ=255)
                                                                                                       cATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAAccc                                               >  1:1613131/1‑150 (MQ=255)
                                                                                                        aTCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAACCGCGACAAGCGTCGCATCCGGCGATATCCCCGGGCGTATTCTTATTGAATCCCATCCCCAAcccc                                              >  2:1480183/1‑150 (MQ=21)
                                                                                                         tCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCg                                             <  1:365138/150‑1 (MQ=255)
                                                                                                                   gCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCccctt                                   >  1:1503735/1‑150 (MQ=255)
                                                                                                                    cACCGCGCCACTGGCGGATGCGGCGTCGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttt                                  >  1:1516252/1‑150 (MQ=255)
                                                                                                                      ccGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttc                                >  1:743604/1‑150 (MQ=255)
                                                                                                                      ccGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttc                                <  2:607794/150‑1 (MQ=255)
                                                                                                                         cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCCCCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                             >  2:768170/1‑150 (MQ=255)
                                                                                                                         cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                             <  2:1155072/150‑1 (MQ=255)
                                                                                                                         cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                             <  2:856176/150‑1 (MQ=255)
                                                                                                                          gcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccctt                            <  1:473158/150‑1 (MQ=255)
                                                                                                                            cACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccctttt                          >  1:1631245/1‑150 (MQ=255)
                                                                                                                              cTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTtct                        <  2:573969/150‑1 (MQ=255)
                                                                                                                                ggCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTATCCCTTTTCTcc                      >  1:944782/1‑150 (MQ=255)
                                                                                                                                 gCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGTAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCg                     <  1:867792/150‑1 (MQ=255)
                                                                                                                                 gCGGATAAGGCGTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCg                     >  2:55770/1‑150 (MQ=39)
                                                                                                                                  cGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCgg                    <  2:1041955/150‑1 (MQ=255)
                                                                                                                                    gataaggcgTCCACGCCGCATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCg                  >  2:517932/6‑150 (MQ=38)
                                                                                                                                         ggcgTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGc             >  1:1311451/1‑150 (MQ=255)
                                                                                                                                           cgTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTa           <  2:279424/150‑1 (MQ=255)
                                                                                                                                              ggACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTGTTCTTTTTGCATCCCATCACGACCCCCGCACTCCCCTTTTCCCTGTTCTCCGGCGACGGGTaact        >  2:682281/1‑148 (MQ=37)
                                                                                                                                               gccGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAAtt       >  1:1521236/3‑150 (MQ=255)
                                                                                                                                                 cGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAg     <  1:190252/150‑1 (MQ=255)
                                                                                                                                                  gCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGa    <  2:460435/150‑1 (MQ=255)
                                                                                                                                                   ccTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGaa   >  1:123008/1‑150 (MQ=255)
                                                                                                                                                   ccTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGaa   >  2:834635/1‑150 (MQ=255)
                                                                                                                                                   ccTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGaa   >  1:1259214/1‑150 (MQ=255)
                                                                                                                                                   ccTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGaa   >  1:1271138/1‑150 (MQ=255)
                                                                                                                                                    cTTATCCGCCGTACATGTGCAATTGCGTAACCTGGATAAAACGCGAGAAGCGTCGCATCCGGCGTTATCAGCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGAAc  <  2:605807/150‑1 (MQ=255)
                                                                                                                                                    cTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCAACGGCTAAATTAGAAc  <  1:1011345/150‑1 (MQ=255)
                                                                                                                                                    cTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCATCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGAAc  <  1:914687/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAATTAGAAC  >  NC_000913/3231324‑3231621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: