Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,250,264 C→T G192D (GGT→GAT)  dhaK ← dihydroxyacetone kinase subunit K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,2640CT92.9% 147.1 / NA 58G192D (GGT→GAT) dhaKdihydroxyacetone kinase subunit K
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (26/26);  minor base G (0/4);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAGCGGTTGTTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCA  >  NC_000913/1250020‑1250500
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
cTGGAGCGGTTGTTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCgag                                                                                                                                                                                                                                           >  1:374433/1‑249 (MQ=255)
 tGGAGCGGTTGTTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGAGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGTGTAACAGAAGGTTTGCCCGCGGCAGGAAAGGGACAGGCATCGAGAGc                                                                                                                                                                                                                                        >  1:242591/1‑251 (MQ=255)
  ggAGCGGTTGTTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGc                                                                                                                                                                                                                                        <  1:166056/250‑1 (MQ=255)
           gtttggtttgTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTGCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTAtt                                                                                                                                                                                                                              <  1:28482/251‑1 (MQ=255)
             ttggtttgTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGAGCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGTTACAGGCATCGAGAGCGATTCCTATTg                                                                                                                                                                                                                             <  2:554859/250‑1 (MQ=255)
                      tgttCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTg                                                                                                                                                                                                                   >  1:223703/1‑251 (MQ=255)
                         tataccttccAATTGCCTTGGTGATAATCCTAGAAACTCAAAGTGCATTGGTATGAGCCATTTACCAGCAGGGTGACGAACATTTCATCGACGGTTTGATCAAGGGAAGATAAGGGTCGGCGGTTAATACCCGGCTCAACATGAATGCCGACGCTCAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGTTTCATGCATCGAGAGCGATACCTATTGAGTTGCCTTGat                                                                                                                                                                                                                 <  2:26634/242‑1 (MQ=255)
                            tcCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGattat                                                                                                                                                                                                              <  1:190098/250‑1 (MQ=255)
                                gCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCa                                                                                                                                                                                                           <  1:536224/249‑1 (MQ=255)
                                  cAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTTCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGGGAGCGATACCTATTGAATGGCCTTGATTATTCAGCtt                                                                                                                                                                                                       >  1:212094/1‑251 (MQ=255)
                                   taCTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCCCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCCGGGGAAAAGAAGGTTTGCCAGCGGAAGGGCCGGTCCAGGCATCGCGAGCGAGACCTATTGAGAGGCCTTGATTATTCAACttc                                                                                                                                                                                                      >  2:145549/2‑250 (MQ=255)
                                   aaCTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTa                                                                                                                                                                                                      <  1:200146/251‑1 (MQ=255)
                                    aCTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGACGGCGGGCAATACCCGGCTCACCATGAATGCAGACGCCACACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGCTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGACTTGATTATTCAGCTTAc                                                                                                                                                                                                     >  1:562261/1‑251 (MQ=255)
                                              tgATAATCCCAGAAACGCAAAGGGCGATGGTATGAGCCATGTACCAGCAGGGGGTCGAACATTTCATCGACGGTTTGAGCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACACCATCTCATGATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCAGCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGtt                                                                                                                                                                                            <  2:228822/250‑1 (MQ=255)
                                               gATAATCCCAGAAACGCAAAGTGCGAGGGTATGAGCCATTTACCAGCAGGGTGTCGATCAGTTCAGCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGGTTGCCCGCGGCAGGAACGGTACAGGCAGCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTcc                                                                                                                                                                                          <  1:384612/251‑1 (MQ=255)
                                                 tAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGAAGGAACGGTAACGGCATCGAGAGCGAAACCTAATGAGTGGCCTTGATTTTTCAGCTTACGCCCCAGTTCCGc                                                                                                                                                                                        >  2:7415/1‑251 (MQ=255)
                                                   aTCCCAGAAACTCAAAGTGCGATTGTATGATCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGAGCAAGGGAAGAGAAGGGGCGGCTGTCAATACCCGGCTCACCAAGAATGACGACGCCCACCTCCATCTCATTAGCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCAGCGAGAGCGATCCCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGcac                                                                                                                                                                                      <  2:4458/251‑1 (MQ=255)
                                                            aCGCAAAGAGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCCTCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAACGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCCTCGAGAGCGATACCTATTGCGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCTCACAGGCGTcc                                                                                                                                                                              >  1:71432/1‑250 (MQ=255)
                                                                 aaGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGATTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGACAAACTCCATCTCATTATCCGCCAGGGTAACAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGAGGCCTTGATTATTCAGCTTACGCACCAGGTCCGCACAGGCGTCCAGTg                                                                                                                                                                          >  1:46487/1‑249 (MQ=255)
                                                                             aTGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTcc                                                                                                                                                            <  1:442555/251‑1 (MQ=255)
                                                                                    aTTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTg                                                                                                                                                     >  1:176038/1‑251 (MQ=255)
                                                                                                          aTTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTTCGCCTACGAGTTTTTCAATTa                                                                                                                                >  1:429025/1‑250 (MQ=255)
                                                                                                              cATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTACAAGAAGGTTTGCCCGCGGCAGGAACGGGACAGGCATCGAAAGCGATACCTATTGAGTGGCCTTGATTATTCAGCGTACGCCCCAGTTCCGCCCAGGGCTCCAGTGAGTCGCCCCGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTATTACg                                                                                                                           >  2:73541/1‑251 (MQ=255)
                                                                                                               aTCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACACCAAACTACATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGCACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGACCCAGTTCCGCACAGGCGGCCAGTGCGTCGCCACGCTCCACCGCTGCGCCTACGGGTTTTTCAATTACTACgg                                                                                                                          >  2:571327/1‑251 (MQ=255)
                                                                                                                  gACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGgtgt                                                                                                                       >  1:276593/1‑251 (MQ=255)
                                                                                                                      gTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGc                                                                                                                   <  1:114515/251‑1 (MQ=255)
                                                                                                                             cAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAAcgccg                                                                                                            <  2:176038/251‑1 (MQ=255)
                                                                                                                                 ggAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGAATATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTc                                                                                                        <  1:408590/251‑1 (MQ=255)
                                                                                                                                 ggAAAAGAAGGGTCGGCGGTGAATACCCGGCTGACCATGAATGCCGACGCCAAACTTCATCTCATTATCCGCCAGGGTAAAAGAAGGGTTGCCCGCCGCAGGAACGGTACAGGCATCGAGATCGATTCCTATGGAGTGGCCTTGATTATTCAGCTTACGCCCCCGTTCCGCCCAGGCGTCCAGTGAGTCGCCACGCGCCGCCGCTGCGCCTGCGAGTTTTGTAATTAATACGGTGTTGGCAACGCCGCGTc                                                                                                        <  2:118038/251‑1 (MQ=255)
                                                                                                                                        aaGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGc                                                                                                 >  1:15556/1‑251 (MQ=255)
                                                                                                                                            ggcggcGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAATGGCCTTGATTATTCAGCTTACGCCCCAGTGCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGta                                                                                             <  2:212094/251‑1 (MQ=255)
                                                                                                                                                  ggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGt                                                                                              >  1:105707/1‑244 (MQ=255)
                                                                                                                                                  ggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGa                                                                                       <  1:27407/251‑1 (MQ=255)
                                                                                                                                                  ggTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCAGCTCATTAGCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCGGCGCCTACGAGTTTTTCAATTAAGACGGTGTTGGCAACGCCGCGTCGCCCGGCAGt                                                                                              <  2:105707/244‑1 (MQ=255)
                                                                                                                                                   gTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGAc                                                                                      >  1:110756/1‑251 (MQ=255)
                                                                                                                                                                   ccaggAATGCCGTCGCCAAACCCCATCGCATTACCCGCAAGGGTAAAAGAAGGTTTGCCCGCGGCTGGACCGGTACAGGCATCGAGAGCGATACCTATTGAGTGTCCTTGATTATGCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCTCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCatc                                                                      <  2:13206/247‑3 (MQ=255)
                                                                                                                                                                     aTGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAAcgt                                                                    >  1:194920/1‑251 (MQ=255)
                                                                                                                                                                       gAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAAcgtcg                                                                  >  1:94429/1‑251 (MQ=255)
                                                                                                                                                                         aTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACgtcgtc                                                                <  1:7415/251‑1 (MQ=255)
                                                                                                                                                                                 cccAAACGCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGACTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTGGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGa                                                        <  2:429025/250‑1 (MQ=255)
                                                                                                                                                                                                aTTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACtt                                        <  2:94429/251‑1 (MQ=255)
                                                                                                                                                                                                    tCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCCTCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCaa                                                                                                                                   >  2:428575/1‑157 (MQ=255)
                                                                                                                                                                                                    tCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCaa                                                                                                                                   <  1:428575/157‑1 (MQ=255)
                                                                                                                                                                                                        ccAGTGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCATGATTATTCAACTTACGCCCCAGTTCCGCACAGGCGTCCAGAGAGTCGACACGCTCCGCCGCTGCGCCTACGAGTTTTTAAATTAATACGGTGTTGGCCACGACGCGTCGCCCGGCAGTATAAAGACTGTCTGTTACCACAAAGGCGCCTTCAATGACACCCAGTGTCACTTTTAcgcc                                 >  2:545208/1‑250 (MQ=255)
                                                                                                                                                                                                         cAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTAcgccgc                               >  1:291362/1‑251 (MQ=255)
                                                                                                                                                                                                            ggTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTAt                            >  1:229741/1‑251 (MQ=255)
                                                                                                                                                                                                               aaaaGAAGGTTTGCCCGCGGCAGGAACGGTTCAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCg                          <  1:496223/250‑1 (MQ=255)
                                                                                                                                                                                                                 aagaagGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGc                       <  2:194920/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGGCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTaa                  <  1:259663/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAAc                 <  1:396047/251‑1 (MQ=255)
                                                                                                                                                                                                                             ccGCGGCAGGAACGGTGCAGGCATCGCGAGCGATACCTCTTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGAGTTGGCAACGCCGCGCCGCCCGGCAGTATAAAGACTGTCTTTTACCGCACCGTCGCCATCAATGACCACAGTGGAAACTTTTACACCGCTATCATGCAGTacatcagtc           >  2:334985/1‑243 (MQ=255)
                                                                                                                                                                                                                              cgcgGCAGGAACGGTACAGGCATCGAGAGCGATACCTAGTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGGCTTTTACCGCAACGTCGTCATCAATGACCACATTGGTCACTTTGTCGCCGCTATCGTGCAGTAACTCGGTCg          <  1:226283/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgggAACGGTACAGGAAGCGAGAGCGATCCCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCTTCCAGTGAGTCGCCACGCTCCGCCGCGGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGGCACTTTTTCGCCGCTATCGTGCAGTAACTCGGTCGCTGt      <  2:562261/248‑1 (MQ=255)
                                                                                                                                                                                                                                   caagaccGGTACAGGCAACGAGACCGATCCCTATTGAGTGGCCTTGATTATTCAGCTTACTCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCGGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGtt     <  2:15556/245‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGCGCAGTAACTCGGTTGCTGtt     <  1:552939/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aGGAACGGAACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGACACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGttt    >  1:360977/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGCGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCCGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCCATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTACCTCGGTCGCTGTTTc   >  2:239098/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaCGGAACAGGCACCGAGGGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCATGTGATCCGCCCCGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGGCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGTCCACAGTGGTCACTGTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCa  <  2:374433/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
CTGGAGCGGTTGTTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCA  >  NC_000913/1250020‑1250500

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