Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,231,460 C→A 8.1% intergenic (+221/‑205) 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,4600CA 8.1% 316.9 / 11.0 98intergenic (+221/‑205)ygjK/fadHglycoside hydrolase/2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base C (41/49);  new base A (2/6);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

Only 100 of 106 total aligned reads displayed.

AACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGG  >  NC_000913/3231311‑3231603
                                                                                                                                                     |                                                                                                                                               
aaCGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgc                                                                                                                                                 >  2:1086834/1‑150 (MQ=255)
  cgcgGCAAGCGTTGCATCCGGCAACTGCCCCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcgg                                                                                                                                               >  1:939076/1‑150 (MQ=255)
  cgcgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcgg                                                                                                                                               <  1:335798/150‑1 (MQ=255)
  cgcgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcgg                                                                                                                                               <  2:48931/150‑1 (MQ=255)
  cgcgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcgg                                                                                                                                               <  2:1312655/150‑1 (MQ=255)
   gcgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGGGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcggc                                                                                                                                              <  1:1070221/150‑1 (MQ=255)
   gcgGCAAGCGGTGCAGCCGGCAACTGCCCCGCGCCCCTGGCGGCGGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcggc                                                                                                                                              <  1:1064694/150‑1 (MQ=21)
    cgGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATgcggcg                                                                                                                                             <  2:180104/150‑1 (MQ=255)
     ggCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGt                                                                                                                                            >  2:365138/1‑150 (MQ=255)
     ggCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGt                                                                                                                                            >  1:573969/1‑150 (MQ=255)
       cAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgg                                                                                                                                          >  1:436936/1‑150 (MQ=255)
       cAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTgg                                                                                                                                          <  1:472860/150‑1 (MQ=255)
             ttGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCt                                                                                                                                    <  2:944785/150‑1 (MQ=255)
                   ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                              >  2:437487/1‑150 (MQ=255)
                   ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                              <  1:407270/150‑1 (MQ=255)
                   ccGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCg                                                                                                                              >  1:200086/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: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)
                          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)
                                        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)
                                             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:1365269/1‑150 (MQ=255)
                                                     gTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAAc                                                                                            <  2:140301/150‑1 (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)
                                                                                    cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                             >  1:900477/1‑150 (MQ=255)
                                                                                    cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                             >  1:511624/1‑150 (MQ=255)
                                                                                    cccGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCg                                                             >  2:1569617/1‑150 (MQ=255)
                                                                                      cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                           <  1:1569617/150‑1 (MQ=255)
                                                                                      cGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCggg                                                           <  2:436936/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                                              <  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                                       <  2:383226/150‑1 (MQ=255)
                                                                                                          aCAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCa                                       <  1:439059/150‑1 (MQ=255)
                                                                                                           cAAGTGTCTCATGCGGCATCTTCACCGCGCCTCTGGCTGATGCGGCGTGTACTCCTTTTCCGCCCTATATGTGCAATTCCGTATCCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCAt                                      <  2:440166/150‑1 (MQ=25)
                                                                                                            aaGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATc                                     <  2:1311567/150‑1 (MQ=255)
                                                                                                              gCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATcac                                   <  1:1352198/150‑1 (MQ=255)
                                                                                                              gCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATcac                                   <  2:78490/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)
                                                                                                                                       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)
                                                                                                                                                     |                                                                                                                                               
AACGCGGCAAGCGTTGCATCCGGCAACTGCACCGCGCCACTGGCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTGTGTTCCCGTAGGTCGGATAAGACGCGACAAGCGTCGCATCCGGCATCTGCACCGCGCCACTGGCGGATGCGGCGTGGACGCCTTATCCGCCCTACATGTGCAATTCCGTAACCTGGATAAAACGCGACAAGCGTCGCATCCGGCGTTATCACCGGGCGTATTCTTTTTGAATCCCATCACAAACCCCGCACTCCCCTTTTCCCTTTTCTCCGG  >  NC_000913/3231311‑3231603

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: