Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,064 T→A 5.5% K351* (AAA→TAA)  dhaM ← dihydroxyacetone kinase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,0640TA 5.5% 348.6 / 7.5 109K351* (AAA→TAA) dhaMdihydroxyacetone kinase subunit M
Reads supporting (aligned to +/- strand):  ref base T (55/48);  new base A (2/4);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 110 total aligned reads displayed.

AGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247916‑1248212
                                                                                                                                                    |                                                                                                                                                    
aGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTa                                                                                                                                                     <  2:696177/150‑1 (MQ=255)
aGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTa                                                                                                                                                     <  2:348680/150‑1 (MQ=255)
    gAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGat                                                                                                                                                 >  2:1079212/1‑149 (MQ=255)
    gAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGaa                                                                                                                                                 >  2:1084900/1‑150 (MQ=255)
       tAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACtt                                                                                                                                              >  1:1538726/1‑150 (MQ=255)
         gTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgc                                                                                                                                            >  2:332517/1‑150 (MQ=255)
           tGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgctg                                                                                                                                          >  1:318064/1‑150 (MQ=255)
            ggCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttagcca                                                                                                                                         <  1:1018881/150‑6 (MQ=255)
            ggCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttagcca                                                                                                                                         <  1:868175/150‑6 (MQ=255)
             gCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGcc                                                                                                                                        <  2:1538726/150‑1 (MQ=255)
              cGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttagccagc                                                                                                                                       <  2:1387748/150‑8 (MQ=255)
                atttAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAgg                                                                                                                                     <  2:295108/148‑1 (MQ=255)
                aGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAgg                                                                                                                                     <  2:974161/150‑1 (MQ=255)
                   tAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCat                                                                                                                                  >  2:109214/1‑150 (MQ=255)
                   tAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCat                                                                                                                                  >  1:1421780/1‑150 (MQ=255)
                      aCTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATAtt                                                                                                                               <  1:1102429/150‑1 (MQ=255)
                          cGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGc                                                                                                                           <  1:1047257/150‑1 (MQ=255)
                            ggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCg                                                                                                                         >  2:1051115/1‑150 (MQ=255)
                            ggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCg                                                                                                                         >  2:1646545/1‑150 (MQ=255)
                             gAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGt                                                                                                                        <  2:318064/150‑1 (MQ=255)
                             gAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGt                                                                                                                        <  1:594439/150‑1 (MQ=255)
                              aGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTg                                                                                                                       <  1:1655885/150‑1 (MQ=255)
                              aGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTg                                                                                                                       >  1:1340819/1‑150 (MQ=255)
                               gTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGc                                                                                                                      <  1:1585687/150‑1 (MQ=255)
                                  cttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   >  2:299458/1‑150 (MQ=255)
                                  cttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCAATTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   >  2:523386/1‑150 (MQ=255)
                                      tttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTc                                                                                                               <  2:1049124/150‑1 (MQ=255)
                                       ttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCa                                                                                                              <  1:1377143/150‑1 (MQ=255)
                                         cGTTTGGGGCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGGTGCTGCTATTGCTGGCTAAGTTCTTTAATAACTTGCTGCCAGGCTTATTCTGCCGTGCAATGTTGATg                                                                                                            <  2:1071285/150‑1 (MQ=255)
                                            ttGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTg                                                                                                         <  2:1113218/150‑1 (MQ=255)
                                              gggTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGaa                                                                                                       <  2:1340819/150‑1 (MQ=255)
                                              gggTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGaa                                                                                                       <  2:1334958/150‑1 (MQ=255)
                                               ggTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      >  2:1607707/1‑150 (MQ=255)
                                                gTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAgg                                                                                                     <  2:1142069/150‑1 (MQ=255)
                                                    ggTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGtt                                                                                                 >  2:699671/1‑150 (MQ=255)
                                                       ggACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGc                                                                                              <  1:1281131/150‑1 (MQ=255)
                                                          ccAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTg                                                                                           <  2:543426/150‑1 (MQ=255)
                                                                tGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccg                                                                                     >  1:1400026/1‑150 (MQ=255)
                                                                 gCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgc                                                                                    >  2:1405400/1‑150 (MQ=255)
                                                                    aTGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcag                                                                                 >  2:413800/1‑150 (MQ=255)
                                                                     tGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagc                                                                                <  2:1223955/150‑1 (MQ=255)
                                                                     tGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagc                                                                                <  2:324899/150‑1 (MQ=255)
                                                                     tGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTACAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagc                                                                                >  1:1353614/1‑150 (MQ=255)
                                                                      gCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               >  1:619107/1‑150 (MQ=255)
                                                                        agaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  1:1521787/1‑150 (MQ=255)
                                                                        agaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  2:209838/1‑150 (MQ=255)
                                                                             aTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCgg                                                                        <  1:457767/150‑1 (MQ=255)
                                                                               cGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      <  1:1646545/150‑1 (MQ=255)
                                                                               cGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      <  2:1521787/150‑1 (MQ=255)
                                                                                gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatc                                                                     >  2:907271/1‑150 (MQ=255)
                                                                                gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatc                                                                     >  1:1147213/1‑150 (MQ=255)
                                                                                gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatc                                                                     <  1:1227288/150‑1 (MQ=255)
                                                                                   cacaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  >  2:902537/1‑150 (MQ=255)
                                                                                    acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCt                                                                 <  1:332517/150‑1 (MQ=255)
                                                                                    acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCt                                                                 <  2:1393933/150‑1 (MQ=255)
                                                                                       tCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAAc                                                              >  2:405559/1‑150 (MQ=255)
                                                                                       tCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAAc                                                              <  1:873407/150‑1 (MQ=255)
                                                                                       tCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAAc                                                              <  1:867760/150‑1 (MQ=255)
                                                                                        cAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACa                                                             <  1:993948/150‑1 (MQ=255)
                                                                                         aaTATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAg                                                            >  2:434138/1‑150 (MQ=255)
                                                                                            ataGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCAGAGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtgt                                                         >  2:1347370/1‑150 (MQ=255)
                                                                                             taGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTa                                                        <  1:699671/150‑1 (MQ=255)
                                                                                                cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATgg                                                     >  1:1181581/1‑150 (MQ=255)
                                                                                                cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATgg                                                     <  1:405559/150‑1 (MQ=255)
                                                                                                cGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATgg                                                     <  1:435871/150‑1 (MQ=255)
                                                                                                 gAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGt                                                    <  1:206893/150‑1 (MQ=255)
                                                                                                  aGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTg                                                   <  2:1400026/150‑1 (MQ=255)
                                                                                                  aGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTg                                                   <  1:654610/150‑1 (MQ=255)
                                                                                                    cTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAc                                                 >  2:1200033/1‑150 (MQ=255)
                                                                                                    cTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGAAGTTCCGGATCATCTAACAGTGTATGGTGAc                                                 >  2:923598/1‑150 (MQ=255)
                                                                                                     ttGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAcc                                                <  1:109214/150‑1 (MQ=255)
                                                                                                      tGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCa                                               >  2:554778/1‑150 (MQ=255)
                                                                                                      tGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCa                                               >  2:542936/1‑150 (MQ=255)
                                                                                                       gTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAg                                              >  2:1546141/1‑150 (MQ=255)
                                                                                                        tAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGa                                             >  1:524010/1‑150 (MQ=255)
                                                                                                            tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAg                                         >  2:664048/1‑150 (MQ=255)
                                                                                                            tatTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAg                                         <  1:413800/150‑1 (MQ=255)
                                                                                                               tcatcaTCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAtt                                      >  1:1499624/1‑150 (MQ=255)
                                                                                                                catcatCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTg                                     >  1:1247556/1‑150 (MQ=255)
                                                                                                                    atcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGc                                 >  2:70898/1‑150 (MQ=255)
                                                                                                                       cAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAcc                                                >  1:344740/1‑132 (MQ=255)
                                                                                                                       cAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAcc                                                <  2:344740/132‑1 (MQ=255)
                                                                                                                           tGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCg                          <  1:907271/150‑1 (MQ=255)
                                                                                                                            gCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGt                         >  1:1628510/1‑150 (MQ=255)
                                                                                                                             cTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTc                        <  2:1147213/150‑1 (MQ=255)
                                                                                                                               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      >  1:1619336/1‑150 (MQ=255)
                                                                                                                               ggTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCaa                      >  1:1625029/1‑150 (MQ=255)
                                                                                                                                   ttGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGccc                  <  1:523386/150‑1 (MQ=255)
                                                                                                                                    tGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCg                 >  2:64017/1‑150 (MQ=255)
                                                                                                                                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGc                >  2:824478/1‑150 (MQ=255)
                                                                                                                                     gCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGACCGc                >  1:388872/1‑150 (MQ=255)
                                                                                                                                      ggcagcaaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCt               >  1:711381/7‑150 (MQ=255)
                                                                                                                                       gcagcaaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTg              >  1:441384/6‑150 (MQ=255)
                                                                                                                                        cagcaaGTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTgg             <  1:339974/146‑1 (MQ=255)
                                                                                                                                         gCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGc            >  2:727643/1‑150 (MQ=255)
                                                                                                                                              gTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTg       >  2:573029/1‑150 (MQ=255)
                                                                                                                                               ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGc      <  2:580217/150‑1 (MQ=255)
                                                                                                                                                tCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCt     >  2:115379/1‑150 (MQ=255)
                                                                                                                                                 cTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  1:1308432/1‑150 (MQ=255)
                                                                                                                                                   ttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTc  <  1:1405400/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247916‑1248212

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