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,2640CT96.8% 200.6 / NA 62G192D (GGT→GAT) dhaKdihydroxyacetone kinase subunit K
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/26);  minor base G (0/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGA  >  NC_000913/1250032‑1250509
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
tttggtttgTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTg                                                                                                                                                                                                                                     >  1:26004/1‑251 (MQ=255)
 ttggtttgTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTg                                                                                                                                                                                                                                     >  1:54370/1‑250 (MQ=255)
       tgttgtTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGcc                                                                                                                                                                                                                              >  1:542671/1‑251 (MQ=255)
        gttgttCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAACGCGATACCTATTGAGTGGCCt                                                                                                                                                                                                                             >  2:581242/1‑251 (MQ=255)
                 ccTGCCAACTGCCTTGTTGATAATCCCATAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCAGTGAGAGCGATACGTAGTGAGTGGCCTTGATTATTc                                                                                                                                                                                                                    <  2:374393/251‑1 (MQ=255)
                    gCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGc                                                                                                                                                                                                                 <  2:506422/251‑1 (MQ=255)
                            ccTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGccc                                                                                                                                                                                                         <  2:54370/251‑1 (MQ=255)
                               tgttgATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCa                                                                                                                                                                                                       <  2:166758/250‑1 (MQ=255)
                               tgttgATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCa                                                                                                                                                                                                       <  1:368453/250‑1 (MQ=255)
                                   gATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCAATATACCAGCAGGGTGTCGAACATTTCATCGCCGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCACTACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGCACAAGCAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTcc                                                                                                                                                                                                  >  2:582259/1‑251 (MQ=255)
                                      aaTCCCAGCAACGCAAAGTGCGCTGGTATGAGCCATTTACAAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGATTGCCCGCGGAAGGCACGGTACAGGCATCGAGAGCGATACCTCTTGAGTGGCCTTGATTATTCAGCTTCAGCCCCAGTTCCGc                                                                                                                                                                                                >  2:106420/1‑250 (MQ=255)
                                          ccAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACGCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAgg                                                                                                                                                                                           <  1:582646/251‑1 (MQ=255)
                                               aaCGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTcc                                                                                                                                                                                      <  1:143040/251‑1 (MQ=255)
                                                             tGGTATGAGCCCTTTACCAGCAGGGTGTCGAACATTTCATCGCAGGTTTGATAAAGGGAAGAGCAGGCGCGGCGGTAAATACCCGGCTCCCCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCCGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGCGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGGGTCGCCACg                                                                                                                                                                        >  2:297452/1‑251 (MQ=255)
                                                                  tGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGCTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCCGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTccg                                                                                                                                                                   >  1:552800/1‑251 (MQ=255)
                                                                   gAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGGAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACTCCAAACTCCATCTCAGTTTGCGCCAGGGTAAAAGAAGGTGTGCCCGCGGCAGGAACGGTACATGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTGACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTccg                                                                                                                                                                   <  1:346611/250‑1 (MQ=255)
                                                                   gAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTccgc                                                                                                                                                                  >  1:378308/1‑251 (MQ=255)
                                                                         tttACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTgc                                                                                                                                                            >  1:439990/1‑251 (MQ=255)
                                                                           tACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGattat                                                                                                                                                                                                                      >  2:29117/1‑191 (MQ=255)
                                                                           tACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGattat                                                                                                                                                                                                                      <  1:29117/191‑1 (MQ=255)
                                                                              cagcagGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCCTTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCt                                                                                                                                                        >  1:162761/1‑250 (MQ=255)
                                                                                       gtCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGtttt                                                                                                                                               <  1:177174/250‑1 (MQ=255)
                                                                                          gAACATTTCATCGACGGTTTGATCAAGGGAAGCGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGGTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCa                                                                                                                                            >  1:464791/1‑250 (MQ=255)
                                                                                          gAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCa                                                                                                                                            >  2:482975/1‑250 (MQ=255)
                                                                                               tttCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAAt                                                                                                                                      >  1:106418/1‑251 (MQ=255)
                                                                                                                   aGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGCCGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGCACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACCCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCcgcg                                                                                                                  >  1:344365/1‑251 (MQ=255)
                                                                                                                     ggAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCAGTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTc                                                                                                                <  2:403414/251‑1 (MQ=255)
                                                                                                                      gAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCg                                                                                                               <  1:582259/251‑1 (MQ=255)
                                                                                                                      gAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTGCAATTAATACGGTGTTGGCAACGCCGCGTCg                                                                                                               <  1:155263/251‑1 (MQ=255)
                                                                                                                         gagaAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGccc                                                                                                            >  1:7813/1‑251 (MQ=255)
                                                                                                                         gagaAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGccc                                                                                                            <  2:194089/251‑1 (MQ=255)
                                                                                                                         gagaAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCAGCGAGAGCGATACCTATTGAGTGGCCTTGATTATGCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGccc                                                                                                            <  2:335616/251‑1 (MQ=255)
                                                                                                                           gaAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCCCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGGAACGCCGCGTCGCCCgg                                                                                                          >  2:429052/1‑251 (MQ=255)
                                                                                                                                ggcggcGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCAATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGGTTATTCCGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGt                                                                                                      >  1:466010/1‑250 (MQ=255)
                                                                                                                                        tCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCACCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACt                                                                                             >  1:385824/1‑251 (MQ=255)
                                                                                                                                         cAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTg                                                                                            >  1:226025/1‑251 (MQ=255)
                                                                                                                                           aTACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTc                                                                                          <  2:106418/251‑1 (MQ=255)
                                                                                                                                              cccGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCCTCGAGAGCGATAACTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACCCCGCGTCGCCCGGCAGTTTAAAGACTGCCttt                                                                                       >  2:385251/1‑251 (MQ=255)
                                                                                                                                                cGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCtttt                                                                                      >  2:11154/1‑250 (MQ=255)
                                                                                                                                                 ggCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCAATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTAc                                                                                    >  1:507876/1‑251 (MQ=255)
                                                                                                                                                  gCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTAcc                                                                                   >  1:422485/1‑251 (MQ=255)
                                                                                                                                                               gCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGtcat                                                                      <  1:315880/251‑1 (MQ=255)
                                                                                                                                                                 cGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGtcatc                                                                     >  1:450123/1‑250 (MQ=255)
                                                                                                                                                                    cGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATg                                                                 <  2:450123/251‑1 (MQ=255)
                                                                                                                                                                    cGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGGCGTCATCAATg                                                                 >  2:466383/1‑251 (MQ=255)
                                                                                                                                                                     gCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATg                                                                 >  1:569909/1‑250 (MQ=255)
                                                                                                                                                                                tctcATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAAAGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTACGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGCCATCACTGACCACAGTGGTc                                                     >  2:258768/1‑251 (MQ=255)
                                                                                                                                                                                       aTCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTa                                              <  1:476096/251‑1 (MQ=255)
                                                                                                                                                                                                gttAAAAGAAGGTTTGCCCGCGGCAGGAACTGTACAGGCATCGAGAGCGATACCGATTGAGTGGCCTTGATTATTCAGCTTACGCCCTAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTTTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGACACTTTTACGCCGCTAt                                     <  2:182201/249‑1 (MQ=255)
                                                                                                                                                                                                  aaaaaGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATGCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCTCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTTTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCg                                   <  2:542671/250‑1 (MQ=255)
                                                                                                                                                                                                     aagaagGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTTATTATTCAGCTTACGCCCCAGGTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCTTCATCAATTACCACAGTGGTCACTTTTACGCCGCTATCGTGc                                <  2:76055/251‑1 (MQ=255)
                                                                                                                                                                                                     aagaagGTTTGCCCGCGGCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGc                                <  2:71023/251‑1 (MQ=255)
                                                                                                                                                                                                          ggtttactcgcgGCAGGAACCTTACAGGCATCGAGGACGATACCTATTGAGTGGCCTTGATTATTCAGATTTAGCCCCAATTCCGCACAGACGTCCATTGGGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGTGGTATAAAGCTTGTCTTTTACCGCCACGTCGTCATCAATGGCCACAATGGTCGCTTTTCCGCCGCTATCATGCAGTa                            >  2:569267/9‑250 (MQ=255)
                                                                                                                                                                                                                 ccgaggCAGGAACGGTACATGCATCGAGAGCGATCCCTATTGAGTGGACTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCATTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTGTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGACTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTc                    <  2:408642/247‑1 (MQ=255)
                                                                                                                                                                                                                     ggCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCCCAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCCATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTg                >  2:83066/1‑251 (MQ=255)
                                                                                                                                                                                                                     ggCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTg                <  1:482975/251‑1 (MQ=255)
                                                                                                                                                                                                                      gCAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGt               <  1:200717/251‑1 (MQ=255)
                                                                                                                                                                                                                       cAGGAACGGTACAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGtt              >  2:39577/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGGTACAGGCATCGAGAGCGATACCTATTGCGTGGCCTTGATTATTCAGCTTACGCCCCAGATCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACAACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCa           >  2:55987/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aCAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTaa    <  1:153355/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cAGGCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACAACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTAGGTCGCTGTTTCAAAGTTa     >  1:286972/1‑249 (MQ=255)
                                                                                                                                                                                                                                   cggCATCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCATCGCCGCGGCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGa  <  2:221607/250‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
TTTGGTTTGTTGTTCTTCCTGCCAACTGCCTTGTTGATAATCCCAGAAACGCAAAGTGCGATGGTATGAGCCATTTACCAGCAGGGTGTCGAACATTTCATCGACGGTTTGATCAAGGGAAGAGAAGGGGCGGCGGTCAATACCCGGCTCACCATGAATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAACGGTACAGGCACCGAGAGCGATACCTATTGAGTGGCCTTGATTATTCAGCTTACGCCCCAGTTCCGCACAGGCGTCCAGTGAGTCGCCACGCTCCGCCGCTGCGCCTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGA  >  NC_000913/1250032‑1250509

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