Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,571 T→C 5.6% T182A (ACA→GCA)  dhaM ← dihydroxyacetone kinase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,5710TC 5.6% 270.6 / 5.0 89T182A (ACA→GCA) dhaMdihydroxyacetone kinase subunit M
Reads supporting (aligned to +/- strand):  ref base T (43/41);  new base C (3/2);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTC  >  NC_000913/1248422‑1248718
                                                                                                                                                     |                                                                                                                                                   
gAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGt                                                                                                                                                     >  1:495003/1‑150 (MQ=255)
 aaaGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                                                                    <  2:102215/150‑1 (MQ=255)
 aaaGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                                                                    <  2:722892/150‑1 (MQ=255)
   aGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGa                                                                                                                                                  <  1:964114/150‑1 (MQ=255)
       aTCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGAt                                                                                                                                                 >  1:259616/1‑147 (MQ=255)
       aTCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGAt                                                                                                                                                 <  2:259616/147‑1 (MQ=255)
       aTCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatg                                                                                                                                              <  1:1295858/150‑3 (MQ=255)
        tCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                                                                             <  1:203562/150‑1 (MQ=255)
        tCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                                                                             <  1:1170461/150‑1 (MQ=255)
                  ttCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAg                                                                                                                                   <  1:1091902/150‑1 (MQ=255)
                   tCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGc                                                                                                                                  >  2:1063956/1‑150 (MQ=255)
                    cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGcc                                                                                                                                 >  1:765740/1‑150 (MQ=255)
                     aGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCg                                                                                                                                <  2:286944/150‑1 (MQ=255)
                     aGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCg                                                                                                                                >  2:1051100/1‑150 (MQ=255)
                      gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGGTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCgg                                                                                                                               <  1:847626/150‑1 (MQ=255)
                       cTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCggg                                                                                                                              >  1:698410/1‑150 (MQ=255)
                         tCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAg                                                                                                                            >  1:1005086/1‑150 (MQ=255)
                            ggCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGcc                                                                                                                         <  1:800326/150‑1 (MQ=255)
                              cccTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGGAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCgg                                                                                                                       >  2:553391/1‑150 (MQ=255)
                                cTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                                                                                     >  2:540587/1‑150 (MQ=255)
                                    cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                                                                 >  2:394703/1‑150 (MQ=255)
                                    cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                                                                 >  1:461210/1‑150 (MQ=255)
                                    cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTAc                                                                                                                 >  2:336883/1‑150 (MQ=255)
                                     gCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                                >  2:1195470/1‑150 (MQ=255)
                                     gCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                                <  1:205659/150‑1 (MQ=255)
                                            ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGc                                                                                                         <  1:1063956/150‑1 (MQ=255)
                                              cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                                <  2:368389/141‑1 (MQ=255)
                                              cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                                >  1:368389/1‑141 (MQ=255)
                                               gcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGt                                                                                                      <  1:88869/150‑1 (MQ=255)
                                                  gCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                                                                                   >  2:541106/1‑150 (MQ=255)
                                                         gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGttttt                                                                                            <  2:714429/150‑1 (MQ=255)
                                                           tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                          <  2:1241194/150‑1 (MQ=255)
                                                           tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                          <  2:1246532/150‑1 (MQ=255)
                                                           tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                          <  1:126027/150‑1 (MQ=255)
                                                             taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATg                                                                                        >  1:544665/1‑150 (MQ=255)
                                                               gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                                      >  1:1196435/1‑150 (MQ=255)
                                                               gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                                      >  1:8948/1‑150 (MQ=255)
                                                               gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                                      >  1:18622/1‑150 (MQ=255)
                                                                       gtagtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCa                                                                              >  1:390187/1‑150 (MQ=255)
                                                                          gtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCagag                                                                           <  1:1416894/150‑1 (MQ=255)
                                                                          gtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCagag                                                                           >  1:174116/1‑150 (MQ=255)
                                                                                cAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACggg                                                                     >  1:664097/1‑150 (MQ=255)
                                                                                     tGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttct                                                                >  1:448302/1‑150 (MQ=255)
                                                                                       gTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCttcttc                                                              <  2:1312589/150‑1 (MQ=255)
                                                                                          aaTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTTTTCATc                                                           <  2:461210/150‑1 (MQ=255)
                                                                                             aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTa                                                        >  2:362931/1‑150 (MQ=255)
                                                                                             aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTa                                                        <  1:336883/150‑1 (MQ=255)
                                                                                             aCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTa                                                        >  1:1039696/1‑150 (MQ=255)
                                                                                              ctctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTAc                                                       >  2:1442236/1‑150 (MQ=255)
                                                                                                 tctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCa                                                    <  1:541106/150‑1 (MQ=255)
                                                                                                 tctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCa                                                    <  2:77381/150‑1 (MQ=255)
                                                                                                   tGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCTGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCCCCAGAGAACGGGCTTCTTCATCGTACGCAgc                                                  >  2:976467/1‑149 (MQ=255)
                                                                                                      tgtgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACa                                               <  1:1051100/150‑1 (MQ=255)
                                                                                                       gtgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAt                                              >  1:1020564/1‑150 (MQ=255)
                                                                                                       gtgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAt                                              >  2:697534/1‑150 (MQ=255)
                                                                                                        tgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtg                                             <  1:1195470/150‑1 (MQ=255)
                                                                                                        tgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtg                                             <  2:1346543/150‑1 (MQ=255)
                                                                                                         gACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtgt                                            >  1:415366/1‑150 (MQ=255)
                                                                                                         gACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtgt                                            >  1:707153/1‑150 (MQ=255)
                                                                                                          aCGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGTTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtgtg                                           <  1:643173/150‑1 (MQ=255)
                                                                                                             cATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCa                                        >  1:928464/1‑150 (MQ=255)
                                                                                                                  gCCGTTTTTTTGCAGCATCTTATCGGCATTAAATGTCGATAAGGTTTAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCTGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAt                                   <  2:286771/150‑1 (MQ=255)
                                                                                                                  gCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAt                                   >  1:441142/1‑150 (MQ=255)
                                                                                                                  gCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAt                                   <  2:635985/150‑1 (MQ=255)
                                                                                                                   ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGCCATGTGTCAGAGAtt                                  >  2:590326/1‑150 (MQ=255)
                                                                                                                   ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAtt                                  >  2:301556/1‑150 (MQ=255)
                                                                                                                   ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAtt                                  >  2:860150/1‑150 (MQ=255)
                                                                                                                   ccGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAtt                                  <  1:949764/150‑1 (MQ=255)
                                                                                                                      tttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCa                               >  1:1458488/1‑150 (MQ=255)
                                                                                                                       ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                              >  1:1258078/1‑150 (MQ=255)
                                                                                                                       ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                              >  1:1328323/1‑150 (MQ=255)
                                                                                                                       ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                              <  2:415366/150‑1 (MQ=255)
                                                                                                                       ttttttACAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACTCAGGACATGTGTCAGAGATTTCAg                              <  2:148417/150‑1 (MQ=255)
                                                                                                                         ttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAATCCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAgtg                            <  2:495003/150‑1 (MQ=255)
                                                                                                                             ccAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCgg                        <  1:433183/150‑1 (MQ=255)
                                                                                                                              cAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCgga                       <  1:1442236/150‑1 (MQ=255)
                                                                                                                                gCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCggagg                     >  2:828997/1‑150 (MQ=255)
                                                                                                                                 taccttatCGACATTTAATGCCGATAAGGTATAAACCAGTCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCggagga                    >  2:411971/6‑150 (MQ=255)
                                                                                                                                 taccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCggagga                    >  1:837931/6‑150 (MQ=255)
                                                                                                                                  aaCATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGAc                   <  1:158344/150‑1 (MQ=255)
                                                                                                                                      tatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTa               <  1:362931/150‑1 (MQ=255)
                                                                                                                                      tatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTa               <  1:553391/150‑1 (MQ=255)
                                                                                                                                           gCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGTCATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAccc          >  2:220940/1‑150 (MQ=255)
                                                                                                                                           gCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAccc          >  1:995230/1‑150 (MQ=255)
                                                                                                                                            cATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCa         >  1:998383/1‑150 (MQ=255)
                                                                                                                                              ttAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgt       <  1:935964/150‑1 (MQ=255)
                                                                                                                                               tAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgtt      >  1:1122639/1‑150 (MQ=255)
                                                                                                                                                   tGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgtt      <  1:1422057/146‑1 (MQ=255)
                                                                                                                                                   tGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgtt      >  2:1422057/1‑146 (MQ=255)
                                                                                                                                                   tGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTc  <  2:698410/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTC  >  NC_000913/1248422‑1248718

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: