Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 705,375 G→A A478T (GCA→ACA)  nagE → N‑acetylglucosamine‑specific PTS enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,3750GA98.3% 186.4 / NA 58A478T (GCA→ACA) nagEN‑acetylglucosamine‑specific PTS enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/31);  minor base C (0/1);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAGCAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCT  >  NC_000913/705130‑705612
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
agatTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACt                                                                                                                                                                                                                                            <  2:271207/249‑1 (MQ=255)
 tatTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTc                                                                                                                                                                                                                                           >  2:108641/1‑251 (MQ=255)
      ctgcGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGc                                                                                                                                                                                                                                       >  2:437198/1‑250 (MQ=255)
       tgcGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCa                                                                                                                                                                                                                                      >  2:71654/1‑250 (MQ=255)
          ggTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTCACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTg                                                                                                                                                                                                                                  >  2:40979/1‑251 (MQ=255)
          ggTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGGGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTg                                                                                                                                                                                                                                  >  1:87490/1‑251 (MQ=255)
              ggcggcACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTGGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCcg                                                                                                                                                                                                                               <  2:27525/250‑1 (MQ=255)
                                 aaGCGATTGACGCCTTTATAACCGGTCTGCGCATTACAGTGGGTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTGCTGGGGTAGTGAAACTGAACAAACAGACTATTAAGGTGATTGTGGGCGCAAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCAGGGTCCGGTAGACGGTGCGCCAGCGGAAACAACTCCGGCAGCTGCCGTGCCTGTGGCAAAGCCgtag                                                                                                                                                                                                           <  2:270057/251‑4 (MQ=255)
                                    cGATAGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGGCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGGCGTTGCCCGTGGTCCGGTAGCCGCTGCGGCAGCTAAAACACCTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAgact                                                                                                                                                                                                        <  1:40979/251‑4 (MQ=255)
                                                    cACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGtat                                                                                                                                                                                        <  1:285230/251‑1 (MQ=255)
                                                          tCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATcg                                                                                                                                                                                  >  2:148835/1‑251 (MQ=255)
                                                             gcgcCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATcg                                                                                                                                                                                  <  2:117820/248‑1 (MQ=255)
                                                                   tACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTgg                                                                                                                                                                         >  1:266868/1‑251 (MQ=255)
                                                                     cAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTa                                                                                                                                                                       <  1:188001/251‑1 (MQ=255)
                                                                         ggCTGACTCTGCCCGCGTTAACGATACGATGCGAAAACGTCTGTGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTAGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCCGACTGTACCAACCGCGGCATCTATCGCGGAGCTGGTATCGc                                                                                                                                                                   >  1:382053/1‑251 (MQ=255)
                                                                              aCTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGAtt                                                                                                                                                              <  1:437198/251‑1 (MQ=255)
                                                                                     ccGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTg                                                                                                                                                       >  2:216731/1‑251 (MQ=255)
                                                                                      cgcgTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTg                                                                                                                                                       >  1:191787/1‑250 (MQ=255)
                                                                                           tAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTc                                                                                                                                                  <  2:112408/250‑1 (MQ=255)
                                                                                                atacgatGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCa                                                                                                                                            <  1:258583/251‑1 (MQ=255)
                                                                                                  acgatGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACt                                                                                                                                          >  1:355951/1‑251 (MQ=255)
                                                                                                   cgatGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTg                                                                                                                                         <  2:336456/251‑1 (MQ=255)
                                                                                                       gtgtAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATc                                                                                                                                     >  2:444311/1‑251 (MQ=255)
                                                                                                                   gggTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGGGATGTCGTGGCACTGGATCAGGTTCCTGACg                                                                                                                         >  1:103438/1‑251 (MQ=255)
                                                                                                                          cctgggggcgggactcTGAACAAACAGACTGTTCAGGTGATTGTGGGCGCGAACGCAGCATCCATCGGCGATGCGATGAAGAAAGTCGTTCCCCGTGGTCCCGTCTCCGACGCGTCAGCTGAAACAACTCCGGCAACTGCCGAGCCTGAAGCAAAACCGCAGACTGTACCAAACGCCGTATCTATCGCGGAGCCGGTATCGCAGATTACCGGTGATGTCGTGCCACTGGATCAGGTTCCTGACGAAGCAtt                                                                                                                  >  2:136770/16‑251 (MQ=255)
                                                                                                                              gggTAGTGAAACGTAAATAACAGAGTATTGAGGTGATGGTTGGCGCGAAAGCAGAATCCATCGGCCATGGGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGTTGCGTCAGCTGAAACATCTCCGGCATCTGCCGCGCGTGTAGCAAAACCGCAGACGGTACCAGACGCGGTATCTATCGCTGAGCTCGTTTCGCCGATTACTGGTCATGTCGTTGCACTGGATCATGCTCGTGACGAACCATTCGcc                                                                                                              <  2:224820/251‑1 (MQ=255)
                                                                                                                                gTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAGGCAGAAGCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAg                                                                                                            <  2:376847/251‑1 (MQ=255)
                                                                                                                                           aaacaaacaGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTATCCGCTGCGTCAGCTGAAACAACTCCGGCATCTGCCGCGCCTGTAGCAAATCCGCATACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTGCCGGTGATGTCGTGGCACTTGATCAGGTTCCTGACGAAGCATGCGCCAGCAAAGCggtt                                                                                                  <  2:223941/249‑2 (MQ=255)
                                                                                                                                               aaacaGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCggtgggtg                                                                                              >  2:206236/1‑250 (MQ=255)
                                                                                                                                                 acaGACTATTCAGGTGATTGTTGGCTCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTAGCTATCGCGGAGCTGGTATTGCCGATTACCGGTGATGTCGTGGCACTGGTTCAGGTTCCTGACGAAGCATTCGCCAGGAAAGCGGTGGGTGACg                                                                                           <  1:363064/251‑1 (MQ=255)
                                                                                                                                                 acaGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACg                                                                                           >  2:18674/1‑251 (MQ=255)
                                                                                                                                                  caGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACgg                                                                                          >  1:235789/1‑251 (MQ=255)
                                                                                                                                                   aGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACggg                                                                                         >  2:214119/1‑250 (MQ=255)
                                                                                                                                                          tCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg                                                                                   <  1:18835/250‑1 (MQ=255)
                                                                                                                                                            aGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTg                                                                                <  2:191787/251‑1 (MQ=255)
                                                                                                                                                               tGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCAATGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGGGACGGTGTGGCGGtcaa                                                                              >  2:116464/1‑247 (MQ=255)
                                                                                                                                                                       ggCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCggggcaaaccacaa                                                                      >  1:263679/1‑238 (MQ=255)
                                                                                                                                                                       ggCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAg                                                                      <  2:145401/250‑1 (MQ=255)
                                                                                                                                                                         cgcgAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATa                                                                   <  1:113360/251‑1 (MQ=255)
                                                                                                                                                                           cgAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATaaa                                                                 <  1:108641/251‑1 (MQ=255)
                                                                                                                                                                            gAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAACCGACACATaaa                                                                 >  1:156139/1‑250 (MQ=255)
                                                                                                                                                                              aaGCAGAATCCATCGGCGATGCTATGAAGAAAGTCGTTTCCCGTGGTCCGGTCGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCTCGCCTGTAGCAAAACCGCAGACTGTACCAATCGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc                                                              <  1:279907/251‑1 (MQ=255)
                                                                                                                                                                                gCAGAATCCAGCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGGCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcgt                                                            <  2:263679/251‑1 (MQ=255)
                                                                                                                                                                                 cAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcgtc                                                           <  1:148835/251‑1 (MQ=255)
                                                                                                                                                                                         aTCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGGCAGATAAAATCGGCGTATAAcc                                                  >  1:419588/1‑251 (MQ=255)
                                                                                                                                                                                                atgcgatgAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa                                           <  2:156139/251‑1 (MQ=255)
                                                                                                                                                                                                 tgcgatgAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGGATcccccaccgc                                            >  1:204914/1‑240 (MQ=255)
                                                                                                                                                                                                 ggcgatgAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAg                                          <  2:235789/250‑1 (MQ=255)
                                                                                                                                                                                                   cgatgAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAggg                                        <  1:135173/251‑1 (MQ=255)
                                                                                                                                                                                                            aaGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAACCGGCAAGTAAAATTGTCGTAATCCCAACCGCAGGGACaaacggg                               >  1:162208/1‑246 (MQ=255)
                                                                                                                                                                                                              gTCGTTGCCCTTGGTCCGGTGGCCGTTGCGTCAGCTGAAACAACTCCGGCATCTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGTATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCTACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaa                             <  1:446923/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccGTGGTCCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg                                                                                   <  2:339056/191‑1 (MQ=255)
                                                                                                                                                                                                                     cccGTGGACCGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCg                                                                                   >  1:339056/1‑191 (MQ=255)
                                                                                                                                                                                                                             ccGGTAGCCGCTGCGTCAGCTGAACCAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACacca              <  2:103438/250‑1 (MQ=255)
                                                                                                                                                                                                                              cGGTAGCCGCTGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaac            <  2:383796/251‑1 (MQ=255)
                                                                                                                                                                                                                                 caGCCGCTGTGTCAGCTGAAACAACACCGGCAACGGCCGCGCTTGTAGCAAAGCCGCAGACTGTTCCAAGCGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGGGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAGCACCAACCAc         <  2:266513/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCGCTGCGACAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCAcgc       <  2:39533/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCGTCAGCTGAAACAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGACTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCCGCCGCAGGGGCAAACGTGAAAATCTTACAAACCAACCAcccggtct  >  1:355874/1‑244 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGGGTAGTGAAACTGAACAAACAGACTATTCAGGTGATTGTTGGCGCGAAAGCAGAATCCATCGGCGATGCGATGAAGAAAGTCGTTGCCCGTGGTCCGGTAGCCGCTGCGTCAGCTGAAGCAACTCCGGCAACTGCCGCGCCTGTAGCAAAACCGCAGGCTGTACCAAACGCGGTATCTATCGCGGAGCTGGTATCGCCGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCT  >  NC_000913/705130‑705612

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: