Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG34.9% 60.8 / 15.6 87G65G (GGT→GGG) ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (18/38);  new base G (20/10);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 2.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433414‑433874
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
aggAATTGGGCGAAGAGCTCGCGCGCCTGGAGGACTAAGCCGGGCAGGAGGATTTTTACATGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACTTTTTCCCCCGCACCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTCCCCGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCg                                                                                                                                                                                                                                  <  2:2290067/237‑1 (MQ=255)
    aTTTGCGGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGAGGAGTATTACATGGCGCGGGCGCAAAAGCTGGCGCAACGAGGACGTTTTCCCACGCCTCCCACCCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTCCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGGGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑Ac                                                                                                                                                                                                                               <  2:806861/237‑1 (MQ=255)
     tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAGCCAAAAGGGCG‑GCCg                                                                                                                                                                                                                             >  1:2020255/1‑238 (MQ=255)
     tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCAGAAGGACACGCGTTGCG‑TATGGC‑GGGCGAAAAAGCCAAAAGCGAG‑ACCg                                                                                                                                                                                                                             >  2:1935651/1‑238 (MQ=255)
      tgggCGAAGAGATCGCGCGCATGAGGGACGAGGCCGGGCGGGACGGGTATACAAGGGCGGGGGCGCTAAAGCGGGCGCAAGGAGGACGATTTCCCACGCATCCCACCCCGAAGGTGGGGTGCGTCATTGTAAAAGAGGGCAAAATGGTCGGTGAAGTTTACCACCACCGTGCGGGTGAACCACTGGCGGAAGTCCCCGCGTTGCG‑TAGGGG‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGcc                                                                                                                                                                                                                           <  2:1893563/237‑1 (MQ=255)
       tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAACCAAAAGGGCG‑ACCGccc                                                                                                                                                                                                                          >  1:629066/1‑238 (MQ=255)
             agagATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACAACAGGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAACCCAAAGTGGG‑Acccactcgttac                                                                                                                                                                                                                    >  2:2094346/1‑229 (MQ=255)
                      gcgcCTGAGAGACTAAGCCGGGCAGGAGGAGTATTACGTGGCGCGGGCGCTAAAGCTGGCGCAAGGAGTCGGTTTTCCCCCGCACCCCAACCCAAATGTGGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAGGTTTCCCACCAACGGGCGGGTAAACCACAGGCCGAAGTACCCGCGTTGCG‑TATGGG‑GGGTAAAAAACCAAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGaa                                                                                                                                                                                                           <  2:1413410/238‑1 (MQ=255)
                              gggaCTAAGCCGTGCAGGACGATTATTACATGGCGCGGGCGCAAAGGCGGGCGAAACGAGGACGTTTTCCCCCGCACCCCACCCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGACCCACATGCCGAAGTACACGCTTTGCG‑TATGGG‑GGGTAAAAAAGCAAAGGGTGCG‑ACCGCCTATGTCACA‑CTCg                                                                                                                                                                                                             <  2:279387/227‑1 (MQ=255)
                                  cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAAGGGGG‑ACCGCCCATTGTACA‑CTCGAACCCTGTAAccct                                                                                                                                                                                               >  1:1699916/1‑236 (MQ=255)
                                        cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAg                                                                                                                                                                                                   >  1:279387/1‑228 (MQ=255)
                                         gTGCGGGACGGTTAATACAGGGCGCGGGCTCTAAGGCGGGCGCAACGAGGGCGTTTTCCCACACATCCCAACCCAAATGGCGGGGGCGGCATTGTAAAAGATGGGGAAATTGTCGGTGAGGGTTACCCCCAACGTGGGGGTGACCCACATGCCGAAGTACGCGCGTTGCG‑TAGGGC‑GGGTGAAAAACCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                        <  2:1544732/238‑1 (MQ=255)
                                             aGGACGGGTTTTACAGGGCCGGGGCTGTAAAGCTGGCGCAACGAGGACGTTTAACCCCGCATCCAAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAGGGTTCCCCCCAACGTGGGGGTGACCCACATGCCGAAGTACACGCGTTGCT‑TATGGG‑GGGTAAAAAAGCAAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTa                                                                                                                                                                                                    <  2:1836599/222‑1 (MQ=255)
                                               agggAGTATTACATGGCGCGGGCGCAAAAGCGGGCGCAACGAGGACGTTTTACCACGCACCCCAACCCGAATGGCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAGGTGCGGGTGAACCACAGGCCGAAGTCCACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                     <  2:33809/232‑1 (MQ=255)
                                                 cGAGTATTACGGGGCGCGGGCGCTAACGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGACCCACAGGCAGAAGTCCACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                <  2:1828790/238‑1 (MQ=255)
                                                 cGAATATTACTTGGCGGGGGCGCTAAAGCTGGCCAAACGGGGCCTTTTTCCCACCCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTAAAGGTTCCCACCAACGGGGGGGTGAACCACATGCCGAAGTACCCGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                <  2:2020255/238‑1 (MQ=255)
                                                   aGTATTACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                        <  1:1612482/228‑1 (MQ=255)
                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAACCAAAAGGGCG‑ACCGCCTCGGCGCCA‑CTCCAACCCCGGAGCACACATGGTCGTAAGCCACCgc                                                                                                                                                                            >  1:2217828/1‑238 (MQ=255)
                                                     tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTg                                                                                                                                                                                                                                    <  1:2093901/184‑1 (MQ=255)
                                                     tATTACATGGCGCGGGCGCAAAAGCGGGCGCAACGAGGGCGTTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCAAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTAAAAAACCCAAAGGTCCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGc                                                                                                                                                                                                  <  2:348749/216‑1 (MQ=255)
                                                      agcacgtGGCGCGGCCGCAAAGGCGGGCGCAACGGAGCCTTTTTCCCACGCCACCCACCCCGAAGGTCGGGTGCGTCATTGTCAAAGAGGGCGAATTTGTCGGTGAAGGTTACCACCAGCGGGGGGGTGACCCACATGCCGAAGTACACGCGTTGCG‑TATGGG‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                           <  2:2258575/232‑1 (MQ=255)
                                                       ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAAGTGCG‑ACCGCCTATTTCACA‑CTCCAACCCTGTAGCCACAATgg                                                                                                                                                                                          >  1:1836599/1‑222 (MQ=255)
                                                         aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAGCCAAAGGGGGG‑ACCGCCCATGTCACA‑CTCGAACCCCGTAGCCATCATGGACGTACGCCAcccg                                                                                                                                                                            >  1:38809/1‑232 (MQ=255)
                                                             ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCCCCCATGTCCCA‑CCCGAACCCTGTAGCCATCATGGGCGTACCCCACCCTg                                                                                                                                                                           >  1:33809/1‑231 (MQ=255)
                                                             ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCCCCCATGTCACA‑CCCGAACCCCGTAGCCATCAAGGTGGTACGCCCCCCAGCtgtg                                                                                                                                                                      >  1:586813/1‑236 (MQ=255)
                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGCGGGGGGAAAAAACCAAAGGT‑CG‑ACCCCCtg                                                                                                                                                                                                                         >  2:2093901/1‑183 (MQ=255)
                                                                 cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGACAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAAACCTGTAGGCATCATGGTCTGACCCCACCGGGCGGTGGCGCaa                                                                                                                                                                >  1:2197743/1‑237 (MQ=255)
                                                                 cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAACCCAAAGGGGCG‑CCCGCCTATGTCACA‑CTCGAACCCCGTAGCCATCATGGTCGGACGCCACCGTGCTGCGACGCAc                                                                                                                                                                >  1:777467/1‑238 (MQ=255)
                                                                    gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAAGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCCGTAGCCATCATGGGCGGACGCCACCGTGCTGGTACGCACTCa                                                                                                                                                             >  1:1266071/1‑238 (MQ=255)
                                                                                             gTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCCCAGGCCGAAGTACCCGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGt                                                                                                                                                                                                                      <  1:442197/157‑1 (MQ=255)
                                                                                                   cccgccATCCCACCCCGAATGTCGGGTGCGTCAGTGTAAAAGAGGGCGAAATGGGCGGTGAAGGTTCCCACCACCGTGCGGGGAAACCACATGCCGAGGTACCCGGTTTGCG‑TATGGC‑GGGTGAAAAAGCAAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                                                                                                             <  2:629066/234‑1 (MQ=255)
                                                                                                       gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAACCAAAGGGGCG‑GCCGCCCCTGTcaccaccgaa                                                                                                                                                                                                            >  2:442197/1‑150 (MQ=255)
                                                                                                        cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCCCCGGCGTAGCGCGCGGGGTTGCCTCCATGCa                                                                                                                        >  1:847998/1‑239 (MQ=255)
                                                                                                         accccaacccaaATTGCGGGGGGGTATATGTAAAAGATGGGGAAATTGTCGGTGAAGGTTACCCCCAACGTGGGGGGAAACCACAGGCAGAAGTACCCGGGTTGCG‑TATGGG‑GGGGAAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                        <  2:670555/236‑1 (MQ=255)
                                                                                                         aTCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑Acc                                                                                                                                                                                                                              <  1:61692/137‑1 (MQ=255)
                                                                                                              aaCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAACCAAAGGTGCG‑ACCGCCCATGTCACA‑CCCCAACCCCGTAGCCATCCTGGTCGTACGCCCCCCGTCTGTGCCGCACCCATCGCCCCCGGCGTAGCGCGCGCGGTTGCCCCCGTGCCAGATc                                                                                                                   >  2:2119153/1‑238 (MQ=255)
                                                                                                                   gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAACCAAAAGGGGG‑ACCCCCCATGTca                                                                                                                                                                                                                    >  2:61692/1‑137 (MQ=255)
                                                                                                                    aaTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAACCAAAAGGGCG‑ACCCCCCATGTCCCACCCCAAACCCTTAACCCACCA‑GGGCGTCACCCACCGCGCGGTGAGGCCCCCACCCCCGCCTGCGTAGCGGGCGGTGGTGCCTCGATTGACGATccccacc                                                                                                             >  2:1268621/1‑233 (MQ=255)
                                                                                                                     aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAACCAAAGGTGCG‑ACCGCCCATGTCACA‑CTCGAACCCCGTAGCCACCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCGCCCCTGGCGGCGCGCGCGGGGTGTCCTCCATGCCAGATCCCAAccc                                                                                                            >  2:2081115/1‑238 (MQ=255)
                                                                                                                      tgtgGGGTGCGTCATTGGAAAAGTGTGCGAAATGTGCGGTGGAGGTTACCACCAGCGTGCGGGTGACCCACAGGCGGAGGTACACGTGTAGCG‑TATGGG‑GGGTGAAAAACCAAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc                                                                                                          <  2:1326779/239‑1 (MQ=255)
                                                                                                                                   atggacAAAGAGGGCAAATGTGGCGGTGAGGGTCACCGCAAAGGGGGGGGAAAACCACAGGCCGAAGTACCCGCGGTGCG‑TAGGGG‑GGGTAAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                              <  2:38809/233‑1 (MQ=255)
                                                                                                                                    ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCGCATGCCGAAGTACACGCGTTGCG‑TAGGGC‑GGGTGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                            <  1:1268621/239‑1 (MQ=255)
                                                                                                                                    ttGACAAAGATGGAGAAATTGTCGGTGAAGGTTCCCACCAAGGGGCGGGTGAACCACATGCCGAAGTACCCGGGTTGCG‑TATGGC‑GGGTGAAAAAGCAAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                            <  2:847998/239‑1 (MQ=255)
                                                                                                                                     tttCAAGGATGCCGAAATTGTCGGTGAGGGTTCCCGCCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCG‑TAGGGC‑GGGTGAAAAAGCCAAAGGTGCG‑TCCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                               <  1:1529267/201‑1 (MQ=255)
                                                                                                                                       tCAAAGATGGCGAAATTGTCGGTGAGGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAGGGTGCG‑CCCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGt                                                                                                                                       <  2:427380/193‑1 (MQ=255)
                                                                                                                                               ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGt                                                                                                                                                                                                                      >  1:1376321/1‑107 (MQ=255)
                                                                                                                                               ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGt                                                                                                                                                                                                                      <  2:1376321/107‑1 (MQ=255)
                                                                                                                                               ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAACCAAAAGTGCG‑GCCGCCCCTGTCACC‑CCCCAACCCCGTAGACATCACGGGCGTACCCCCCCGTGCTGTGTCGCCCCCACCGCCCCCCGCGCACCGCGCGGGGTTGCCCTCATGCACGAGCCCACACCCCACGTCGCCGGGCgggggggttact                                                                                  >  2:1471686/1‑230 (MQ=255)
                                                                                                                                               ggCGAAATTGTCGGGGAGGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                               <  1:37367/193‑1 (MQ=255)
                                                                                                                                                 cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGGGTAGCGCGCGGGGTTGCCTCg                                                                                                                             >  1:427380/1‑193 (MQ=255)
                                                                                                                                                      ttGTCGGTGAAGGTGACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCG‑TATGGC‑GGGTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                                          <  1:2081115/239‑1 (MQ=255)
                                                                                                                                                         tCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAACCAAAGGGGCG‑ACCGCCCATGTCACC‑CCCGAAACCCGTAGCCACCATGGGCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCCCGCTggcaga                                                                                                                     >  2:37367/1‑188 (MQ=255)
                                                                                                                                                            gTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGAACCCGCTTTGCG‑TATGGC‑GGGTAAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTa                                                                                                                                                                                                    <  1:1740907/111‑1 (MQ=255)
                                                                                                                                                             agAAGGTTACCACAAGCGTGCGGGTGACCCACATGCCGAAGTACCCGCTTGGCG‑TAGGGC‑GGGTGAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGt                                                                                                                                                                                                     <  2:881786/108‑1 (MQ=255)
                                                                                                                                                                      ccaccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGTGCG‑ACCCCCCATGTCACA‑CTCGAACCCCGTAGGCAACATgg                                                                                                                                                                                          >  2:1740907/1‑111 (MQ=255)
                                                                                                                                                                       caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGGAAAAGCCAAAGGTGGG‑ACCGCCCATGTCCCA‑CTCGGACCCTGTAGCCAACATg                                                                                                                                                                                           >  1:881786/1‑109 (MQ=255)
                                                                                                                                                                            aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAACCAAAAGGGCG‑AACGCCCATg                                                                                                                                                                                                                       >  1:1204498/1‑77 (MQ=255)
                                                                                                                                                                             ggggggggTGACCCACAGGCCAAAGTCACCGCGTTGCG‑TAGGGC‑GGGTGAAAAAGCAAAAGGGGCG‑CCCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:350098/234‑1 (MQ=255)
                                                                                                                                                                             cGTGCGGGTGAACACCAGGCCGAAGTCCCCGCGTTGCG‑AATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:1261089/239‑1 (MQ=255)
                                                                                                                                                                                gCGGGAGACCCACATGCCAAAGGACA‑GCGGTGGGTTAGGGA‑GGGTAAAAAACCCAATGGGCCA‑ACCGCCTATTTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                 <  2:178805/238‑1 (MQ=255)
                                                                                                                                                                                  ggggGAACCACAGGCCGAAGTACACGCTTTGCG‑TATGGC‑GGGTAAAAAAGCCAAGGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                               <  1:1280923/238‑1 (MQ=255)
                                                                                                                                                                                  gggTAACCCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTa                                                                                                                                                                                                    <  2:2247044/89‑1 (MQ=255)
                                                                                                                                                                                    gtaaaCCCCAGGCAGAAGTACCCGCGTTGGG‑GATGGC‑GGTAAAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCGCA‑CTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                                                   <  1:2058003/101‑1 (MQ=255)
                                                                                                                                                                                    gTGAACCACATGCCGAAGTACACGCGTTGCG‑TTTGGG‑GGGTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                           <  2:1928937/112‑1 (MQ=255)
                                                                                                                                                                                       aaCCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGt                                                                                                                                                                                                                      <  2:1237890/67‑1 (MQ=255)
                                                                                                                                                                                       aaCCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGTGCG‑ACCGCCTATGt                                                                                                                                                                                                                      >  1:1237890/1‑67 (MQ=255)
                                                                                                                                                                                            caTGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGGGCG‑ACCCCCTATGTCACA‑ATCGAACCCTGTAACCATCATgg                                                                                                                                                                                          >  1:2247044/1‑89 (MQ=255)
                                                                                                                                                                                             aTGCCGAAGTACCCGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCcatc                                                                                                                                                                                              <  2:320081/84‑1 (MQ=255)
                                                                                                                                                                                              tGCCGAAGTCCGCGCGTTGCG‑TAGGGG‑GGGTGAAAAAGCCAAGGGGGCG‑GCCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTg                                  <  1:46666/239‑1 (MQ=255)
                                                                                                                                                                                              tGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGGAAAAGACAAAAGGGCG‑GACCCCTATTTCACA‑CACGAACCCTGGAGCCCTCATTGTCGTACACCACCGTGCt                                                                                                                                                                         >  2:2058003/1‑104 (MQ=255)
                                                                                                                                                                                              tGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGGAGCCATTATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                                                 >  1:1928937/1‑112 (MQ=255)
                                                                                                                                                                                              tGCCGAAGTACACGCGTTGCG‑TAGGGC‑GGGTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                <  1:1108339/97‑1 (MQ=255)
                                                                                                                                                                                              ggCCGAAGTACACGCGTTGGG‑TATGGG‑GGGTGAAAAAGCAAAGGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                   <  2:221903/173‑1 (MQ=255)
                                                                                                                                                                                               gCCGAAGTCCCCGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAGGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                               <  2:144873/145‑1 (MQ=255)
                                                                                                                                                                                                      tACACGGGTTGCG‑TAGGGC‑GGGTAAAAAACCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGt                                                                                                                                    <  1:695502/133‑1 (MQ=255)
                                                                                                                                                                                                      tACACGCGTTGCG‑TATGGG‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                <  1:509826/89‑1 (MQ=255)
                                                                                                                                                                                                       acacGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGGCGTAc                                                                                                                                                                                    >  1:320081/1‑84 (MQ=255)
                                                                                                                                                                                                        cacGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCCCCCATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGTGCtgtg                                                                                                                                                                      >  2:1108339/1‑97 (MQ=255)
                                                                                                                                                                                                        cacGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGGCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGa                                                                                         >  1:221903/1‑174 (MQ=255)
                                                                                                                                                                                                         acGCGTTGCG‑TATGGC‑GGGGGAAAAAGCCAAAGGGGCG‑ACCGCCTATTTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                                                                                                     >  1:144873/1‑145 (MQ=255)
                                                                                                                                                                                                            gtgtggctaggggC‑GGGAGAAACAGCAAAGGGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTGGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                            <  2:777467/221‑1 (MQ=255)
                                                                                                                                                                                                                gCG‑TATGGC‑GGGGGAAAAAGCCAAAGGGGCG‑ACCCCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACCCATCGCCCCTGGCGTAGCGCGCGGGGTTGCCTCGATg                                                                                                                          >  2:695502/1‑133 (MQ=255)
                                                                                                                                                                                                                             gggAAAAAACCAAAAGGGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCCGTGACGCACCCCTCGCCGCTGGCGTAACGCGCCTGGTTGCCTCGAGACAAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTACCGGCTGCAACAGGCTGGCCTTGACGGCAGCCCCGGCCCGGATATGAGTGTAGGCGGGCCATTGGATAAAAGCCTtct      >  2:203783/3‑238 (MQ=255)
                                                                                                                                                                                                                                aaaacGCCAAGGGGCCG‑ACCGCCTATGTCACC‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                       <  2:1437796/218‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaaGCCGAAGGTGCGCCCCGCCTATGTCTAG‑CTCGAACCCTGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg  <  2:591101/239‑1 (MQ=255)
                                                                                                                                                                                                                                          aGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAg             >  1:1437796/1‑218 (MQ=255)
                                                                                                                                                                                                                                           ggTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg  <  2:1789553/228‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
AAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCG‑TATGGC‑GGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433414‑433874

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