Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,231,467 G→C 9.3% intergenic (+228/‑198) 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,4670GC 9.3% 295.4 / 17.7 97intergenic (+228/‑198)ygjK/fadHglycoside hydrolase/2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (39/47);  new base C (3/6);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAAT  >  NC_000913/3231318‑3231615
                                                                                                                                                     |                                                                                                                                                    
cAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgg                                                                                                                                                      >  1:436936/1‑150 (MQ=255)
cAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgg                                                                                                                                                      <  1:472860/150‑1 (MQ=255)
      ttGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCt                                                                                                                                                <  2:944785/150‑1 (MQ=255)
            ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                                          <  1:407270/150‑1 (MQ=255)
            ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                                          >  1:200086/1‑150 (MQ=255)
            ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                                          >  2:437487/1‑150 (MQ=255)
             cGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGGCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGc                                                                                                                                         <  2:291494/150‑1 (MQ=18)
              ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                        <  2:214730/150‑1 (MQ=255)
              ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGcc                                                                                                                                        >  1:1041955/1‑150 (MQ=255)
              ggCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGcc                                                                                                                                        <  1:1409682/150‑1 (MQ=21)
               gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGccc                                                                                                                                       <  1:136472/150‑1 (MQ=255)
               gCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGccc                                                                                                                                       <  1:1033424/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)
                   cTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccagtgg                                                                                                                                   <  2:1163507/150‑6 (MQ=18)
                     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:630494/150‑1 (MQ=255)
                                 tGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTa                                                                                                                     <  2:615690/150‑1 (MQ=255)
                                 gggCGGATAAGGCGTTCACGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTa                                                                                                                     >  2:772405/2‑150 (MQ=17)
                                  ggCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTaa                                                                                                                    >  2:14487/1‑150 (MQ=255)
                                     ggataaggcgTTCACGCCGCATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCt                                                                                                                 >  2:985080/7‑150 (MQ=17)
                                      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                                                                                                                                    >  2:168777/1‑113 (MQ=17)
                                                   cGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATAAGGCGTCCACGCCGCATCCGccagtg                                                                                                                                    <  1:168777/117‑5 (MQ=17)
                                                         aTCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGt                                                                                             >  2:1507732/1‑150 (MQ=255)
                                                              cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcg                                                                                                       >  2:288338/1‑135 (MQ=255)
                                                              cccTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAcg                                                                                                       <  1:288338/135‑1 (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                                                                         >  2:1569617/1‑150 (MQ=255)
                                                                             cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                                         >  1:511624/1‑150 (MQ=255)
                                                                               cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                       <  2:436936/150‑1 (MQ=255)
                                                                               cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                       >  2:1609887/1‑150 (MQ=255)
                                                                               cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                       <  1:1609887/150‑1 (MQ=255)
                                                                               cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                       >  2:523439/1‑150 (MQ=255)
                                                                               cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                                       <  1:1569617/150‑1 (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                                                          <  2:568/150‑1 (MQ=255)
                                                                                            aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTg                                                          >  1:766662/1‑150 (MQ=255)
                                                                                            aGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTg                                                          <  1:784320/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                          <  2:607794/150‑1 (MQ=255)
                                                                                                                            ccGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttc                          >  1:743604/1‑150 (MQ=255)
                                                                                                                               cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCCCCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                       >  2:768170/1‑150 (MQ=255)
                                                                                                                               cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                       <  2:856176/150‑1 (MQ=255)
                                                                                                                               cgcCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCcccttttccct                       <  2:1155072/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)
                                                                                                                                     tGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTtctc                 <  1:1365269/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)
                                                                                                                                                     |                                                                                                                                                    
CAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGGCGACGGCTAAAT  >  NC_000913/3231318‑3231615

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: