Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,988 T→C 9.4% intergenic (+595/+120) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,9880TC 9.4% 162.5 / 6.6 53intergenic (+595/+120)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (28/20);  new base C (2/3);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCA  >  NC_000913/2304826‑2305211
                                                                                                                                                                  |                                                                                                                                                                                                                                
gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAAGCCGCATCCCACAACTAAAGCCCGAGCCCGTTGCCCGATTGGACGCttg                                                                                                                                                                  >  1:1046023/1‑225 (MQ=11)
                         tgtgAAGTGCCACACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTTTGAAGTGCTCCCCCCCGTAGGGTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGcc                                                                                                                                                                                             <  1:1432274/175‑1 (MQ=255)
                         ggtgAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACg                                                                                                                                               <  1:1260684/220‑1 (MQ=37)
                                               gTCGGATGAGGCGGTTCCGCCGCGTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                  <  1:1343719/228‑1 (MQ=255)
                                                   gATAAGGCGGTTACGCCGCATCCGACATCTACCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTa                                                                                                              <  1:1456962/228‑1 (MQ=255)
                                                             ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTa                                                                                                        <  1:148821/224‑1 (MQ=255)
                                                                   cgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGGCCGGCTTATCACCATAAGCCCCATCCGGCCAACCATGCCCGCTGTGGCGGTTCCTTACAAGGCaa                                                                                                >  2:911442/1‑226 (MQ=14)
                                                                     cATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                                                                         >  2:1370450/1‑167 (MQ=35)
                                                                      aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCCACATCTAACGCCCCAGCCGGGTGCCTGATGCGAAGCTGGCGCGGCTTATCAGCATACGCCACACTCGCGAACTCATGCCGGATGTGGCGGATTATTTCAAccgaat                                                                                               >  2:736175/1‑219 (MQ=14)
                                                                      aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                        <  1:1372920/167‑1 (MQ=35)
                                                                      aTCCGACAGCTAGCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCCGGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATc                                                                                           <  1:585632/227‑1 (MQ=255)
                                                                       tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACCCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGCGGCGTATCATTTCAACGCACTATCCg                                                                                         >  2:725385/1‑228 (MQ=255)
                                                                       tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATATAACCCCCGAGCCGGTTGCCTGATGCGACCCTGGCGCGTCTTTTCCGCATACGCCCCATTCGGCATACCATCTCGGGTGTGGGGTATCCtt                                                                                                         >  2:614983/1‑213 (MQ=14)
                                                                        ccGACATCTAACGCCCAAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTa                                                                                                        <  1:615991/213‑1 (MQ=255)
                                                                                  aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCCGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCtttt                                                                              >  1:1156121/1‑228 (MQ=255)
                                                                                    gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTg                                                                            <  1:726620/228‑1 (MQ=255)
                                                                                       cgaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                                                                                                                                                            >  1:719404/3‑130 (MQ=11)
                                                                                       cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtgtc                                                                                                                 >  2:601101/1‑187 (MQ=255)
                                                                                        gaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATg                                                                                                                                                                           <  2:718184/129‑1 (MQ=14)
                                                                                              gTTGCCTGAGGCGACGCTGGCGCGTCTTGTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGc                                                                  <  1:72362/228‑1 (MQ=255)
                                                                                               ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTTCCTGATGCGACGCTGGCCCGGCTTATCCGCATACGcccc                                                                                                                                           >  2:1959421/1‑153 (MQ=255)
                                                                                               ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCCCAACCGGCCTACCCTGCCGGGTGGGGGGTATCCTTTCAACCCAATAATCGCCCACTACTTTTTAACCGGGTGGc                                                                  >  2:1181806/1‑228 (MQ=17)
                                                                                               ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                       >  2:434451/1‑175 (MQ=255)
                                                                                                tGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGTAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                      <  1:435173/175‑1 (MQ=255)
                                                                                                tGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCcaca                                                                                                                                          <  1:1962915/155‑1 (MQ=255)
                                                                                                      aTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCCTACCATTCCGGATGTGGCGGATCATTTCAACCCAATATCCGCCACTTCTTTTTGAACCGGTTGCGGTTTCaa                                                          >  2:1622296/1‑228 (MQ=255)
                                                                                                      aTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATTCCGGATGTGGCGTTTCATTTCAACGCAATATCCGCACCTTCCTTTTTAACGGGTTTCGGGTTCaa                                                          >  2:1148138/1‑228 (MQ=21)
                                                                                                        gCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTa                                                                                                        <  1:1492227/181‑1 (MQ=255)
                                                                                                          caCGTTGGCGCTTCTTACCGGGCTGCCACGCGTGAGAAGGCCTCCCCCCCGGAGTTTGGATAGGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTg                                                       <  2:1044138/226‑1 (MQ=25)
                                                                                                                      tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCCCTTCTTTTTGAACGGGTTGCGGTTTCAACTTCGGTTTCg                                               >  1:1411653/1‑223 (MQ=255)
                                                                                                                          tatcggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCg                                                                                                                                                                                           >  2:1715419/5‑80 (MQ=255)
                                                                                                                           atcggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGa                                                                                                                                                                                          <  1:1718449/77‑1 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTg                                                                                                                                                                   >  1:1555942/1‑89 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATa                                                                                                                                 >  1:1576015/1‑123 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTTTAACTACCACGACATATCCGCCACTTCTTTTTCAACCGGTTGCGGTTTTAACTGCGGTTTCGGCGTAACACctgctactctcggcc                       >  2:1352340/1‑218 (MQ=21)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCCCTTCTTTTTGAACGGGGTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCt                       >  1:407147/1‑228 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAAAGGGTTGCGGTTTCAACAGCGGGTTCGCCGTACTATATGCTGCCTGCGGCt                       >  2:1636704/1‑228 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGCATACGCCACATCCGGCATACCCTGCCGGATGTGGCGTATCATTACAACCCCATATCCGCCACTTCTTTTTGAACGGGGTGCGGTTTCAAATGCGGTTTCGGCGGACTATCTGCTGCCTGCGGcc                       >  2:1923322/1‑227 (MQ=255)
                                                                                                                                         cTGTGAAGTGCTCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAAc                                                                                                    >  2:1427034/1‑152 (MQ=18)
                                                                                                                                          tgtgAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                  <  2:1553175/89‑1 (MQ=255)
                                                                                                                                          tgtgAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                                                                                                <  2:1573218/123‑1 (MQ=255)
                                                                                                                                          tgtgAAGTGCTCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACg                                                                                                   <  1:1429591/152‑1 (MQ=18)
                                                                                                                                            taaaGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGt                    <  1:1298144/226‑1 (MQ=255)
                                                                                                                                               aGTCCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTTTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGGGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGt                 <  1:1412177/228‑1 (MQ=255)
                                                                                                                                                 tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCCCATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTTCGGCTTGTCGtc                >  1:1481663/1‑226 (MQ=255)
                                                                                                                                                 tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGGATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGGTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTa                >  1:1727854/1‑227 (MQ=255)
                                                                                                                                                 tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGCATCATTACAACGCACTATCCGCCACCTCTTTTTGAACGGTTTGCGGTTTTACCTCCGGTTTCGGCGTACTATCTCATGCCTGCGGCTTGtacct                >  2:902205/1‑223 (MQ=37)
                                                                                                                                                   cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGGTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGt              >  1:198949/1‑227 (MQ=255)
                                                                                                                                                    tcaaCCCCGTAGGTCGGATAAGCCGGTTACGCGGCTGCCGACATTTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTc            <  1:485795/225‑1 (MQ=38)
                                                                                                                                                      cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCACTATCCGCCACTTCTTTTTGAACGGGTTGCGGGTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGGAGTTCa           >  1:397850/1‑227 (MQ=255)
                                                                                                                                                           cGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGGAGtttcacc        >  1:653478/1‑220 (MQ=255)
                                                                                                                                                             tAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGc   >  1:133084/1‑228 (MQ=255)
                                                                                                                                                              aGGTTGGAAAAGGCGGTTCCCCCCCCTCCAACATCTAACCCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCa  <  2:1154052/228‑1 (MQ=255)
                                                                                                                                                                  tGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACAt                                                                                                                                         >  1:530303/1‑90 (MQ=255)
                                                                                                                                                                  |                                                                                                                                                                                                                                
GCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCA  >  NC_000913/2304826‑2305211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: