Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,252,810 T→G 5.5% intergenic (+22/+85) yeiI → / ← nupX putative kinase/nucleoside permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,252,8100TG 5.5% 364.4 / 5.7 117intergenic (+22/+85)yeiI/nupXputative kinase/nucleoside permease
Reads supporting (aligned to +/- strand):  ref base T (51/59);  new base G (5/2);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCG  >  NC_000913/2252663‑2252959
                                                                                                                                                   |                                                                                                                                                     
ggATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTa                                                                                                                                                     >  1:182973/1‑150 (MQ=255)
 gATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTATCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTaa                                                                                                                                                    >  1:1175988/1‑150 (MQ=255)
  aTACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAc                                                                                                                                                   <  1:1702332/150‑1 (MQ=255)
   tACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAcg                                                                                                                                                  <  2:544045/150‑1 (MQ=255)
    aCAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAcgc                                                                                                                                                 >  2:78692/1‑150 (MQ=255)
     cAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAATCAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAcgcg                                                                                                                                                >  1:288222/1‑150 (MQ=255)
        tttCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCt                                                                                                                                             <  1:1332022/150‑1 (MQ=255)
        tttCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCt                                                                                                                                             <  1:1218488/150‑1 (MQ=255)
         ttCCGCGACCGCGCCCGTTTTGCGGTAGCCTGCGCGGCATTTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCGGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCtt                                                                                                                                            >  2:1635126/1‑150 (MQ=255)
          tCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTa                                                                                                                                           <  1:1688181/150‑1 (MQ=255)
             gcgACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACACCCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGAGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATcc                                                                                                                                        >  1:1205906/1‑150 (MQ=255)
               gACAGCGCCCGTTTTGGGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCgg                                                                                                                                      <  2:1150772/150‑1 (MQ=255)
               gACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCTCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGGTTTCATTAGTGCCAATGGTGTAACTATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCgg                                                                                                                                      <  1:1560270/150‑1 (MQ=255)
               gACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCgg                                                                                                                                      <  1:120935/150‑1 (MQ=255)
               gACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCgg                                                                                                                                      <  1:272312/150‑1 (MQ=255)
                 cAGCGCCCGTTTTTCGGTAGCCTGCGCGGGAATTTCGCGAGCCATCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGGGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGcc                                                                                                                                    <  2:1367700/150‑1 (MQ=255)
                 cAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGcc                                                                                                                                    <  1:5195/150‑1 (MQ=255)
                 cAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGcc                                                                                                                                    <  1:676227/150‑1 (MQ=255)
                  aGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCt                                                                                                                                   >  1:1177356/1‑150 (MQ=255)
                  aGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCt                                                                                                                                   <  2:815756/150‑1 (MQ=255)
                   gcgcCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCGGcaac                                                                                                                                  <  1:1187840/150‑4 (MQ=38)
                    cgcCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTAc                                                                                                                                 >  2:623680/1‑150 (MQ=255)
                       ccGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGGTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCGTAACGCGTCTTATCCGGCCTACata                                                                                                                              >  1:350864/1‑148 (MQ=255)
                         gTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATc                                                                                                                            >  2:1251653/1‑150 (MQ=255)
                           tttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAg                                                                                                                          <  1:871283/150‑1 (MQ=255)
                           tttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAg                                                                                                                          <  1:366929/150‑1 (MQ=255)
                            ttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAgg                                                                                                                         <  2:875119/150‑1 (MQ=255)
                            ttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAgg                                                                                                                         <  1:505114/150‑1 (MQ=255)
                            ttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAgg                                                                                                                         <  2:1822668/150‑1 (MQ=255)
                               cGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTACCGCGTCTTATCCGGCCTACAAATCAGGCaa                                                                                                                      >  1:1634465/1‑150 (MQ=255)
                                 gTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCGAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTAGAAATCAGGCAAtt                                                                                                                    >  1:996210/1‑150 (MQ=255)
                                    gCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCa                                                                                                                 <  2:1288870/150‑1 (MQ=255)
                                       tGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCa                                                                                                              <  2:1000406/150‑1 (MQ=255)
                                       tGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCa                                                                                                              <  2:182973/150‑1 (MQ=255)
                                         cgcgGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAg                                                                                                            >  2:896138/1‑150 (MQ=255)
                                           cgGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTa                                                                                                          <  1:150116/150‑1 (MQ=255)
                                               aaTTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGAGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGAGGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATg                                                                                                      <  1:331727/150‑1 (MQ=255)
                                                aTTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGt                                                                                                     >  1:852351/1‑150 (MQ=255)
                                                 tttCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTa                                                                                                    >  2:1373996/1‑150 (MQ=255)
                                                  ttCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAg                                                                                                   >  2:617688/1‑150 (MQ=255)
                                                   tCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAACTTCAGCAAGTAAATGTAgg                                                                                                  <  2:1355026/150‑1 (MQ=255)
                                                     gcgcgcCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGcc                                                                                                <  2:540651/150‑1 (MQ=255)
                                                         gcCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGAt                                                                                            <  1:1583365/150‑1 (MQ=255)
                                                          ccAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACATTGTTGCCGGATGCGACGCTTAACGCGTCTTATCCTGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATa                                                                                           <  1:103692/150‑1 (MQ=255)
                                                          ccAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATa                                                                                           >  2:427774/1‑150 (MQ=255)
                                                          ccAGCGGCAGCTTAAACAACCCTACACTGTCTGCCGATATCGCGCTTTCATTAGTTCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCTGCCTACAAATCAGGCAATTTCAGCAATTAAATGTAGGCCTGATa                                                                                           <  1:275214/150‑1 (MQ=255)
                                                             gaggCATCTTAAAGAACCCTACCCTGTCTGCCGATATCGCGCTTTCATTAGTGCCAATGGGGTTACAATGTTGACGGATGCGACTCTTAACGCGTCTTATCATGACTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGc                                                                                        <  2:1524197/148‑1 (MQ=255)
                                                              cGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCATCAAGTAAATGTAGGCCTGATAAGCt                                                                                       >  2:678808/1‑150 (MQ=255)
                                                                gCAGCGTAAACAACCCTACCCTGTCTGCCGATAACGCACTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTg                                                                                     <  1:129126/150‑1 (MQ=255)
                                                                 cAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTgc                                                                                    <  2:1227881/150‑1 (MQ=255)
                                                                    cTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGAGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCTATTTCAGCAAGTAAATGTAGGACTGATAATCTTGCGCa                                                                                 >  1:919778/1‑150 (MQ=255)
                                                                       aaaCAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCa                                                                              >  2:766092/1‑150 (MQ=255)
                                                                        aacaacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAg                                                                             <  1:78692/150‑1 (MQ=255)
                                                                            accctcccctGTCTGCCGATAACGCGCTTTCATTAGTGACAATGGTGTAACAATGTTGACGGATGCGACGCTTAACGCGTCTTATACGGCCTACAAAACAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTAGCGCATCAGGCaa                                                                         >  2:1397944/1‑150 (MQ=255)
                                                                            accctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCaa                                                                         >  1:1678391/1‑150 (MQ=255)
                                                                             ccctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCAAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAAt                                                                        >  1:1082953/1‑150 (MQ=255)
                                                                              cctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATg                                                                       <  1:1747499/150‑1 (MQ=255)
                                                                                        cTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGc                                                             <  1:528628/150‑1 (MQ=255)
                                                                                         tGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCt                                                            <  1:24781/150‑1 (MQ=255)
                                                                                          gCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTa                                                           <  2:1610840/150‑1 (MQ=255)
                                                                                          gCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTa                                                           >  1:1377068/1‑150 (MQ=255)
                                                                                          gCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTa                                                           >  2:570052/1‑150 (MQ=255)
                                                                                             gATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAAc                                                        <  2:263956/150‑1 (MQ=255)
                                                                                              aTAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAAcc                                                       >  2:563940/1‑150 (MQ=255)
                                                                                                 aCGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAAt                                                    <  2:674268/150‑1 (MQ=255)
                                                                                                  cgcgCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATa                                                   >  1:562739/1‑150 (MQ=255)
                                                                                                  cgcgCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATa                                                   >  1:981698/1‑150 (MQ=255)
                                                                                                    cgcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCCGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTATCCCAATaaa                                                 >  2:68152/1‑150 (MQ=255)
                                                                                                    cgcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATaaa                                                 <  1:439908/150‑1 (MQ=255)
                                                                                                    cgcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATaaa                                                 <  2:207658/150‑1 (MQ=255)
                                                                                                    cgcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATaaa                                                 >  2:308239/1‑150 (MQ=255)
                                                                                                     gcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAg                                                >  1:808655/1‑150 (MQ=255)
                                                                                                           aTTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCg                                          >  1:596795/1‑150 (MQ=255)
                                                                                                            ttAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCgg                                         >  1:388326/1‑150 (MQ=255)
                                                                                                            ttAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCgg                                         >  1:779527/1‑150 (MQ=255)
                                                                                                              aGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCa                                       >  2:1566555/1‑150 (MQ=255)
                                                                                                                 gCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATg                                    >  1:1075249/1‑150 (MQ=255)
                                                                                                                        gtgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCAc                             <  1:1251653/150‑1 (MQ=255)
                                                                                                                        gtgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCAc                             <  1:896138/150‑1 (MQ=255)
                                                                                                                        gtgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCAc                             <  2:1283597/150‑1 (MQ=255)
                                                                                                                           tAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACtca                          >  2:1248601/1‑150 (MQ=255)
                                                                                                                                    tGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTgg                 >  2:1659339/1‑150 (MQ=255)
                                                                                                                                      ccGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaa               <  2:852351/150‑1 (MQ=255)
                                                                                                                                      ccGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaa               <  1:692855/150‑1 (MQ=255)
                                                                                                                                      ccGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaa               <  2:1177356/150‑1 (MQ=255)
                                                                                                                                       cGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCATGGTAAACCTATATAGAACCCGGCTATGGTCGCACTCATCAAGTTGGaaa              <  1:1635126/150‑1 (MQ=255)
                                                                                                                                       cGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaaa              <  2:10596/150‑1 (MQ=255)
                                                                                                                                       cGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaaa              <  2:1037642/150‑1 (MQ=255)
                                                                                                                                       cGGATAAGACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaaa              >  2:932067/1‑150 (MQ=38)
                                                                                                                                        ggataagACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAg             >  2:463970/7‑150 (MQ=255)
                                                                                                                                        ggataagACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAg             >  1:359843/7‑150 (MQ=255)
                                                                                                                                          aTGCGACGCTTAACGCGTCTTATCCTGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAAGGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCg           <  1:71823/150‑1 (MQ=255)
                                                                                                                                            gCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTc         >  2:1504609/1‑150 (MQ=255)
                                                                                                                                            aagACGCGTTAAGCGTCTTATCCGGCCTCCAACTCAGGCAATATCAGCAAGTAACTGTAGGCCTGATAACCTTGCGCATCAGGCAATGACACACAAGCTAAACCAATAAAGCAACCGGCAATGGTCGCACTCATCAACTTGGAAAGCGTc         >  2:879428/3‑150 (MQ=255)
                                                                                                                                             cGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcg        >  1:657555/1‑150 (MQ=255)
                                                                                                                                             cGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcg        >  2:448108/1‑150 (MQ=255)
                                                                                                                                             cGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcg        >  2:1585166/1‑150 (MQ=255)
                                                                                                                                              gACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcgc       >  2:1363073/1‑150 (MQ=255)
                                                                                                                                                cGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcgccg     >  1:145958/1‑150 (MQ=255)
                                                                                                                                                   ttAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcgccgccg  <  2:152952/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCG  >  NC_000913/2252663‑2252959

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