Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 639,562 T→G 5.4% intergenic (+54/‑191) ahpC → / → ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913639,5620TG 5.4% 967.0 / 11.2 301intergenic (+54/‑191)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (100/177);  new base G (21/2);  total (122/179)
Fisher's exact test for biased strand distribution p-value = 1.86e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TACGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAG  >  NC_000913/639413‑639711
                                                                                                                                                     |                                                                                                                                                     
tACGTAGCTTCTCACGCAGGTGAAGTTTGCCAGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGt                                                                                                                                                       <  1:4983650/150‑1 (MQ=255)
tACGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGt                                                                                                                                                       <  1:4009806/150‑1 (MQ=255)
    tAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCact                                                                                                                                                   >  1:2441524/1‑150 (MQ=255)
       cTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttac                                                                                                                                                <  1:798476/150‑1 (MQ=255)
       cTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttac                                                                                                                                                <  2:1397687/150‑1 (MQ=255)
         tctcCCCCAGGTGACGTGGGCCCGGCTAAATGGAAAGACGGTGACGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAAGTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                                                                                                              <  2:1881140/150‑1 (MQ=255)
         tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                                                                                                              <  1:1096140/150‑1 (MQ=255)
         tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                                                                                                              <  2:3549470/150‑1 (MQ=255)
         tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                                                                                                              <  2:3328643/150‑1 (MQ=255)
         tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                                                                                                              <  1:2862846/150‑1 (MQ=255)
         tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCGGTCacttactt                                                                                                                                              >  1:1008075/1‑150 (MQ=255)
           tcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtgt                                                                                                                                            <  1:1154251/150‑1 (MQ=255)
           tcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtgt                                                                                                                                            >  2:3860036/1‑150 (MQ=255)
           tcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCAGGCTCCGTCTCTGGACCTGGATGGTAAATTATAAATTTCCTTCGTATTTCAAGCCATAGCGGCGTTGGCGTCGCCAGCTAACCCCCCTCACTAACAtgt                                                                                                                                            >  1:3297315/1‑150 (MQ=255)
             aCCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTgtgt                                                                                                                                          <  2:4403345/150‑1 (MQ=255)
               ccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTaa                                                                                                                                        <  1:4612654/150‑1 (MQ=255)
               ccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCTCACCCCGGTCACATACTTGTGTaa                                                                                                                                        >  1:3566178/1‑150 (MQ=255)
                cAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACATACTTGTGTAAg                                                                                                                                       >  1:368523/1‑150 (MQ=255)
                 aGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCGGGCACGTAATTGTGTAAGc                                                                                                                                      >  1:4029178/1‑150 (MQ=255)
                    tGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCGCTCACATAAATTTGTAAGcgac                                                                                                                                   >  1:456887/1‑147 (MQ=255)
                    tGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCACACCCCCGGCACATACATGTGTAAGCTcc                                                                                                                                   >  1:4904244/1‑150 (MQ=255)
                     gAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTccc                                                                                                                                  >  2:3442281/1‑150 (MQ=255)
                     gAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCCGTCACTTACTTGTGTAAGCTccc                                                                                                                                  >  2:2031849/1‑150 (MQ=255)
                       aGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCgg                                                                                                                                <  2:3020863/150‑1 (MQ=255)
                       aGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAAGTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCgg                                                                                                                                <  2:3620572/150‑1 (MQ=255)
                        gTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCTCACCCCGGTCACATAATTGTGTAAGCTCCCcgg                                                                                                                               >  1:1291403/1‑150 (MQ=255)
                         tttGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCGGTCACTTACATTTGTAAGCTCCCcggg                                                                                                                              >  1:33168/1‑150 (MQ=255)
                           tGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCCGGCACATACTTGTTTTAGCACCCggggtt                                                                                                                            >  1:2256564/1‑148 (MQ=255)
                                ggCTAAATGGACAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                                                                                                                       <  2:4458792/150‑1 (MQ=255)
                                ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                                                                                                                       <  2:1528520/150‑1 (MQ=255)
                                ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                                                                                                                       <  1:4126803/150‑1 (MQ=255)
                                ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                                                                                                                       <  2:2760300/150‑1 (MQ=255)
                                ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCCGTCACGTAAATTTGTTAGCTCCCGGGGGTTcaca                                                                                                                       >  2:2427001/1‑150 (MQ=255)
                                   tAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGc                                                                                                                    >  2:3167200/1‑150 (MQ=255)
                                     aaTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTa                                                                                                                  <  2:4745869/150‑1 (MQ=255)
                                     aaTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTa                                                                                                                  <  2:1609711/150‑1 (MQ=255)
                                     aaTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCCTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCCGGCCCATACTTGGGTAAGCGCCCGGGGATTCACAGGCTa                                                                                                                  >  2:4321472/1‑150 (MQ=255)
                                       tGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCCGGCACATACTTGTGTAAGCTCCCGGGGCTTCACCGGCTAgc                                                                                                                >  1:3524643/1‑150 (MQ=255)
                                       tGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGTGGGTCGCCCGCGCACCCCCGTCACATACTTTTGTGAGCTCCCGGGGATTAAAAGGCGAgc                                                                                                                >  1:3476042/1‑150 (MQ=255)
                                         gAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGGGCGTCGCCCCCACACCCCCGTCACATACATGTGTAAGCGCCCGGGGAGTAACAGGCGAgccg                                                                                                              >  1:1958238/1‑150 (MQ=255)
                                          aaaGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGGGTCGCCCCCTCACCCCCGTCACATAATTGTGTAAGCTACCGGGGATTCACAGGCGAgccgc                                                                                                             >  1:692154/1‑150 (MQ=255)
                                          aaaGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCCGTCACATAATTTTGTTAGAGCCCCGGGAGTAACAGGCTAgccgc                                                                                                             >  1:2327168/1‑150 (MQ=255)
                                             gaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCtt                                                                                                          <  2:2213448/150‑1 (MQ=255)
                                               agGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCGCCCCGGGGTTAACCGGCTAGCCGCCGTGc                                                                                                        >  2:2323765/1‑150 (MQ=255)
                                                ggTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGct                                                                                                       <  2:405311/150‑1 (MQ=255)
                                                 gTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc                                                                                                      <  2:860757/150‑1 (MQ=255)
                                                 gTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc                                                                                                      <  2:1923731/150‑1 (MQ=255)
                                                 gTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc                                                                                                      <  1:1865241/150‑1 (MQ=255)
                                                 gTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc                                                                                                      <  1:517341/150‑1 (MQ=255)
                                                       cAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGGGACGCCCCCGAACCCCCCACACTTAAAATTGTAATCACCAGGGGATTAACAGGCTATCACCCATGGTCTGaacc                                                                                                >  1:4015445/1‑147 (MQ=255)
                                                         aCTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGa                                                                                              <  2:3127186/150‑1 (MQ=255)
                                                         aCTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGa                                                                                              <  2:1943464/150‑1 (MQ=255)
                                                          ctctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGGGTGGCCCCCCCACCCCCCGCCCATAAATGTGTAATATACCCCGGATTGACACGCTAGCCGCCATTGTCGGACGCGaa                                                                                             >  1:1373418/1‑150 (MQ=255)
                                                                cTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACtt                                                                                       <  2:2556778/150‑1 (MQ=255)
                                                                cTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACtt                                                                                       <  1:4320158/150‑1 (MQ=255)
                                                                 tCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTc                                                                                      >  2:4618064/1‑150 (MQ=255)
                                                                 tCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTc                                                                                      <  2:1965633/150‑1 (MQ=255)
                                                                  ccGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCg                                                                                     <  1:43926/150‑1 (MQ=255)
                                                                  ccGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCg                                                                                     <  1:4696817/150‑1 (MQ=255)
                                                                       tctGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAAt                                                                                >  2:826520/1‑150 (MQ=255)
                                                                       tctGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAAt                                                                                <  2:294204/150‑1 (MQ=255)
                                                                         tGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACATACTTGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACATAGGAAAATc                                                                              >  2:4843975/1‑150 (MQ=255)
                                                                           gACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCAc                                                                            <  2:2901027/150‑1 (MQ=255)
                                                                              cTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTa                                                                         >  2:2743116/1‑150 (MQ=255)
                                                                                     gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                  <  2:2169115/150‑1 (MQ=255)
                                                                                     gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                  <  2:187690/150‑1 (MQ=255)
                                                                                       aaaaTCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCgg                                                                <  1:1778388/150‑1 (MQ=255)
                                                                                        aaaTCTAAATGTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCggg                                                               <  1:850087/150‑1 (MQ=255)
                                                                                          aTCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGtg                                                             <  2:574187/150‑1 (MQ=255)
                                                                                            cTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTCGCGTAGCCCGCTCAACCCGGTCACTTACTTGTGGACGCTCCCGGGGATTCACAGGCAAGCCGCCGTGCTCTGGAGCGAAATAATTAGGAAATTTACCGAATTCTTCGGGggcg                                                           >  2:653947/1‑149 (MQ=255)
                                                                                               aaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCg                                                        <  1:1561473/150‑1 (MQ=255)
                                                                                               aaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCGGTCACGGACGTGTGGCACCGCCCCGGGCGTCACCGGCGCGCCCCCGTGCTCGGACGCGAAATACGTCGGAAAGTCACCTAATTCGTCGGGGgcggcg                                                        >  2:3982757/1‑150 (MQ=255)
                                                                                                   tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                    <  2:133661/150‑1 (MQ=255)
                                                                                                   tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                    <  2:3699068/150‑1 (MQ=255)
                                                                                                      ttCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACGTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCg                                                 >  2:4131417/1‑150 (MQ=255)
                                                                                                      ttCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACATGTGTAAGCTCCCGGGGATTCACCGGCTAGCCGCCTTGCTCTTACGCGAAAAAATTAGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCg                                                 >  2:1573945/1‑150 (MQ=255)
                                                                                                         gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                              <  2:4310286/150‑1 (MQ=255)
                                                                                                         gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                              <  2:1640856/150‑1 (MQ=255)
                                                                                                            tttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTCATTCTTCGGGTGCTGCGGCACCCGATTttc                                           >  2:5075907/1‑150 (MQ=255)
                                                                                                            tttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAtt                                              >  2:1115991/1‑147 (MQ=255)
                                                                                                              tCACGCCATAGCGGCGTTGGCGTCGCCCCCGCAACCCCGTCACTTACTTGTGTAAGCGCCCCGGGATTCACAGGCTAGCCGCCGTGCTCAGCACCGAAATACTTCGGAAAATCACCGAATTCTTCGGGTGCGGCGGCACCCCAATttatt                                         >  2:3039397/1‑150 (MQ=255)
                                                                                                                  gCCATAGCGGCGTTGGCGTCGCCCCCTCACCCCGGTCACGTACTTGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCGTGCTCGGACCCGAAATACTTCGGAAAGTCACCTCATTCGTCGGGTGCGGCGGCACCCCAGTTTGTTcccc                                     >  2:81161/1‑150 (MQ=255)
                                                                                                                      tAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCAc                                 <  1:264934/150‑1 (MQ=255)
                                                                                                                           gcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCatga                            <  2:1731891/150‑1 (MQ=255)
                                                                                                                            cgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCatgat                           <  2:4904244/150‑1 (MQ=255)
                                                                                                                             gTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTTGCCGCCTTGCTCTGAGGGGAAATACTTCGGAAATTCACCTATTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCGCCGCACCatgatg                          <  1:3046457/150‑1 (MQ=255)
                                                                                                                              ttGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGc                         <  1:2718430/150‑1 (MQ=255)
                                                                                                                              ttGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGc                         >  1:1576838/1‑150 (MQ=255)
                                                                                                                                ggCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCaa                       <  1:31818/150‑1 (MQ=255)
                                                                                                                                 gCGTCGCCCGCTCACCCCGGTCACTGACGTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACTTCGGACATTCACCTAAGTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGAGGCAAg                      <  2:246525/150‑1 (MQ=255)
                                                                                                                                   gTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCt                    <  1:1225472/150‑1 (MQ=255)
                                                                                                                                         cGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATcc              <  1:2399811/150‑1 (MQ=255)
                                                                                                                                           cTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg            >  1:2923523/1‑150 (MQ=255)
                                                                                                                                              aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGACATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCAGGATGCAAGCTGCATCCAGGtc         >  2:180142/1‑149 (MQ=255)
                                                                                                                                              aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCACGCTGCTTCCAGGTa         >  1:611244/1‑150 (MQ=255)
                                                                                                                                              aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTa         >  1:3738187/1‑150 (MQ=255)
                                                                                                                                              aCCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTa         >  2:701067/1‑150 (MQ=255)
                                                                                                                                               ccccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAg        >  2:2696202/1‑150 (MQ=255)
                                                                                                                                                cccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGc       <  2:3368769/150‑1 (MQ=255)
                                                                                                                                                 ccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGcc      <  2:1291403/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAG  >  NC_000913/639413‑639711

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