Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,561 G→A 5.6% A185V (GCC→GTC)  dhaM ← dihydroxyacetone kinase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,5610GA 5.6% 278.0 / 3.4 88A185V (GCC→GTC) dhaMdihydroxyacetone kinase subunit M
Reads supporting (aligned to +/- strand):  ref base G (45/38);  new base A (3/2);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCC  >  NC_000913/1248413‑1248710
                                                                                                                                                    |                                                                                                                                                     
cagcTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGc                                                                                                                                                      >  2:150954/1‑150 (MQ=255)
    tGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCGGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGCTGTGACGTATTTCAAGTTTTTTTACAGCAACATATCGGCATTa                                                                                                                                                  >  2:935964/1‑150 (MQ=255)
     gACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                                                                                                                 >  1:77381/1‑150 (MQ=255)
     gACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                                                                                                                 >  1:722892/1‑150 (MQ=255)
       cGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAAt                                                                                                                                               >  2:158344/1‑150 (MQ=255)
       cGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAAt                                                                                                                                               >  2:1416894/1‑150 (MQ=255)
         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:1170461/150‑1 (MQ=255)
                 tCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                                                                     <  1:203562/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                                                                                                        <  1:205659/150‑1 (MQ=255)
                                              gCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                        >  2:1195470/1‑150 (MQ=255)
                                                     ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGc                                                                                                 <  1:1063956/150‑1 (MQ=255)
                                                       cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                        >  1:368389/1‑141 (MQ=255)
                                                       cgcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACa                                                                                                        <  2:368389/141‑1 (MQ=255)
                                                        gcAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGt                                                                                              <  1:88869/150‑1 (MQ=255)
                                                           gCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTAc                                                                                           >  2:541106/1‑150 (MQ=255)
                                                                  gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGttttt                                                                                    <  2:714429/150‑1 (MQ=255)
                                                                    tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                  <  1:126027/150‑1 (MQ=255)
                                                                    tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                  <  2:1241194/150‑1 (MQ=255)
                                                                    tataGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                  <  2:1246532/150‑1 (MQ=255)
                                                                      taGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATg                                                                                >  1:544665/1‑150 (MQ=255)
                                                                        gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                              >  1:1196435/1‑150 (MQ=255)
                                                                        gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                              >  1:18622/1‑150 (MQ=255)
                                                                        gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                              >  1:8948/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:707153/1‑150 (MQ=255)
                                                                                                                  gACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtgt                                    >  1:415366/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                           <  2:635985/150‑1 (MQ=255)
                                                                                                                           gCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAt                           >  1:441142/1‑150 (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:1328323/1‑150 (MQ=255)
                                                                                                                                ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                      >  1:1258078/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)
                                                                                                                                                    |                                                                                                                                                     
CAGCTGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCC  >  NC_000913/1248413‑1248710

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: