Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	433,649	T→G	36.6%	G65G (GGT→GGG)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	36.6%	24.1 / 10.1	57	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base G (15/7);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

AGAATTTGGCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433415‑433871                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  aGAATTTGGGGAAGAG‑TCGCGGCGCCTGGAGTACTAAGCCGGCGGGGAAGAGTATTAAAGGGCGGGGGCGCTAAAGCGGGGGCAC‑C‑GAGGACGTTTCACCACGCACCCCAACCCGAATGG‑CGGGTGCATCATTGTCAAAGATGGCGAAA‑TGGTCGGTGAAGTTTACCACCAA‑CGTGCGGGTGAACC‑AAATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGa                                                                                                                                                                                                                              <  2:276285/239‑1 (MQ=255)     tttGGCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGGGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAAC‑CAAAAGG‑GGGAccc                                                                                                                                                                                                                           >  1:492089/1‑237 (MQ=255)     tttGGCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGAAGAGTTTTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGAGGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑AAAAGGT‑GCGACCGc                                                                                                                                                                                                                          <  1:1032546/239‑1 (MQ=255)       tGGCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGGATGGC‑GGGTGGAAAAGC‑CAAAAGG‑GCGCCGGccc                                                                                                                                                                                                                        >  1:772222/1‑238 (MQ=255)        ggCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAGAAAGC‑CAAAGGG‑GGGACCGcc                                                                                                                                                                                                                         >  1:963111/1‑237 (MQ=255)                          ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGC‑CAAAGGG‑GCGACCGCCTagttccacccccaac                                                                                                                                                                                                         >  1:807831/1‑222 (MQ=255)                              ggggaCTAAGCCGGGCAGGGAGAGTTTTACAAGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GGGGGCGTTTTACCACGCAGCCCAACCCGAATGT‑CGGGTGCGTCATTGTAAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAG‑CGTGGGGGTGAACC‑AAATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑AAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTa                                                                                                                                                                                                  <  1:1241767/236‑1 (MQ=255)                                    taggcggggaggggggggttttACAGGGGGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGA‑TTTTTACCCCGCATCCCACCCCGAATTG‑CGGGTGCGTCATGGCCAAAGAGGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGGGCGGGTAAACC‑ACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCcatca                                                                                                                                                                                           <  2:1182181/218‑1 (MQ=255)                                           tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAATC‑CCAAAGG‑GGGACCCCCTATTGCCCA‑CTCGGACCCTgatgccccactgagcg                                                                                                                                                                                     >  1:454453/1‑223 (MQ=255)                                           tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGC‑CCAAGGT‑GCGACCGCCCATTTCACA‑CTCGAACCCGGTAGCcatccatggcag                                                                                                                                                                                    >  1:117070/1‑231 (MQ=255)                                           tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGC‑CAAAGGG‑GCGACCCCCCATTTCCCA‑CTCGAAACCTGTAGcctgtccctttt                                                                                                                                                                                     >  1:376933/1‑228 (MQ=255)                                           tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGCAAAAGC‑CCAAAGA‑GGGGCCCCCCACAGCACA‑CTCGCAccccataacaaacaacggtga                                                                                                                                                                                    >  2:669240/1‑221 (MQ=255)                                             ggaggtgggtttttactGGCGCGGGGGCTAAAGGTGGGGGAA‑CGGGGGACGTTTTCCCCCGCATCCAACCCCGATGGT‑CGGGTGCGGTATAGTAAAGGGGGACAAAG‑TGGTGTGAGAAGGTTACCCCCAA‑CGTGGGGGTGAACC‑ACATGCGGAAGTCCCCGCGGTGCTTTGGGG‑GGGTAAAAAAGC‑AAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                   <  2:24573/223‑1 (MQ=255)                                             cAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGC‑CAAAGGG‑GCGACCCCCCATGTCCCA‑CCCCAACCCTGTAGCcatccatgacgta                                                                                                                                                                                   >  1:568770/1‑229 (MQ=255)                                                     gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGC‑CAAAGGG‑GCGACCGCCTATGTCCAA‑CTCGAAACCTGTAGCCACCATGGTCGTaacgccccc                                                                                                                                                                           >  1:174216/1‑230 (MQ=255)                                                      ttttACATGCCGCGGTCGTTAAAG‑TGGGACAAGC‑GGGGTCGTTTGACCACGCACCCCACCCCGAATGG‑GGGGTGCGGTATTGACAAAGAGGGAGAAA‑TGGTGGGGGAGGTTTACCACCAC‑GGGGCGGGTAACCC‑CAAGGCCGAATTCCCGGCGTGGCGTATGGG‑GGGTAAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                         <  2:395717/237‑1 (MQ=255)                                                        ttaacaggCGGGGGCGGTAAAGGTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCCCCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCt                                                                                                                                                                       <  1:236545/233‑1 (MQ=255)                                                          aCATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGG‑GCGACCGCCCATGTCCCA‑CTCGAACCCTTGAGCCATCATGGCGGTACGCCCCCGTGCtgt                                                                                                                                                                     >  1:70454/1‑239 (MQ=255)                                                              ggCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAAC‑CAAAAGG‑GCGACCGCCCATGTCACA‑CTCGCACCCCGTAGCCATACTGGTCGTACGCCACCGGGCTGTGAc                                                                                                                                                                  >  1:359011/1‑238 (MQ=255)                                                               gggcgGGCGTAAAAGGTGGCGCAA‑C‑GGGGCGGTTTTACCACGCATCCAAACCCGAATGG‑CGGGGGCGGCATTGTCAAAGATGGCGAAA‑TTGTCGATGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGAc                                                                                                                                                                  <  1:1155546/235‑1 (MQ=255)                                                                cgcgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAGC‑CAAAGGT‑GCGACCGCCCATGTCACA‑CTCGAACCCTGTAGCCACCAGGGTGGTACGCCCCCggcgctggacgc                                                                                                                                                                >  1:19772/1‑227 (MQ=255)                                                                cgcgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGC‑CAAAAGG‑GGGACCGCCTATTTCCCA‑CTCGGACCCTGTAGCCATCATGGGCGTACGCCACCGTGCCGTGACGc                                                                                                                                                                >  1:368718/1‑238 (MQ=255)                                                                cgcgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAAC‑CAAAGGG‑GCGACCCCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTAAGCCACCGGGCTGGGACGc                                                                                                                                                                >  1:1126619/1‑238 (MQ=255)                                                                  cgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGG‑GCGACCGCCTATGTCACA‑CTCGAACCCCGAAACCCTCATGGTCGGACCCCCCCCTGCTGTGACGccc                                                                                                                                                              >  1:1324212/1‑236 (MQ=255)                                                                  cgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAAC‑CAAAAGT‑GCGACCCCCCATGTCACA‑CTCGAACCCCGTAGCCATCATGGGCGTACGCCCCCCGCCTGTGCCGCAc                                                                                                                                                              >  1:508101/1‑238 (MQ=255)                                                                  cgGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGATGAAGGTTACCACCAA‑CGTGCGGGGGAACC‑ACATGCCGAAGTACACGCGTTGCGAATGGCGGGGGGAAAAACA‑CAAAGG‑‑GGGACGCCCTATGTCACA‑CCCCAACCCTCTAGCCCCAATGGTCGTACCCCACCCCGCTGTGGCGAAc                                                                                                                                                              >  2:1155546/1‑238 (MQ=255)                                                                                ggCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGC‑CCAAAGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCAAGGGCGTACCCCCCCGTGCTGTGACGCACTCATCCCCCCTGGCg                                                                                                                                               >  1:530131/1‑239 (MQ=255)                                                                                          ggggtcGTTTTCCCCCCCCTCCCACCCCGAATGG‑GGGGGGGAGCATTGCAAAAGAGGGGAAAG‑TTGCCGGAGAAGTTTACCACCAA‑CGGGGGGGGGAACC‑ACGGGCGGAAGTACCCGCGTTGCGTAGGGG‑GGGTGAAAAAGC‑AAAAGGT‑CCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcgc                                                                                                                                       <  2:508101/234‑1 (MQ=255)                                                                                                   ttccccccccccccACCCCGAATGG‑GGGGTGCGTAATTG‑CAAAAGTGGGGAAA‑TTTGCGGGTGGGGTTTCCCCCAAGCGGGGGGGGAAACC‑CCATGCCAAAGCCCACGCTTGGCGTTGGGG‑GGGTAAAAAAGC‑CAAGGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGcc                                                                                                                              <  2:1369922/228‑1 (MQ=255)                                                                                                           ctccccAACCCGATTGG‑GGGGTGGTTCTTTGTCAAAGGTGGCAAAT‑TTGGCGGTGAGGGTTGCCCCCAC‑CGGGCGGGTGACCC‑ACGGGCGGAAGTACGCGCGTTGCGTAGGGC‑GGGAGAAAAAGC‑AAAGGGG‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                      <  2:70454/236‑1 (MQ=255)                                                                                                           caccccAA‑CCGAATGGCCGGGGGCGCAATTGTCAAGAATGGGAAAA‑TTGCCGGTGAAGGTTCCCACCAA‑CGGGCGGGTAAACC‑CAATGCCGAAGTCACCGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                      <  1:955844/236‑1 (MQ=255)                                                                                                           cATCCAAACCGGAATGG‑CGGGGGCGTCATGGTAAAAGAGGGCGAAA‑TTGTCGGTGAAGGTTACCCCCAA‑CGTGGGGGGAACCC‑CCATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑AAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                      <  1:753640/239‑1 (MQ=255)                                                                                                             ccccccccagatactgggggTGGTTCAGTGTAAAGGAGGGAGAAA‑TTGGTGGTGATGTTTCCCACCAA‑GGTGCGGG‑GAACCGACTGGCCGAAGTCCGCGTGTTGCGAAGGGG‑GGGTGAAAAATC‑AAGAGGG‑CAGACCGCCTATGTCACA‑CTCGAACCCGGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCaa                                                                                                                     <  2:142203/222‑1 (MQ=255)                                                                                                                        aTGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATAATGGTGGTACCCCCCCGGGCTGTTACGCACCCCACGCCCCTGGCGTAGCGCGCGTGGGTGCCCCGAGGCAAAGTCCTAAccc                                                                                                          >  1:1175961/1‑238 (MQ=255)                                                                                                                           tggggggtggtttttGTAAAAAGTGGGGAAATTTGGCGGGGAAGTTTCCCCCCAA‑GGGGCGGGTAACCC‑ACATGCCGAAGTACCCGGGTTGCGTATGGG‑GGGAAAAAAAGC‑AAAGGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc                                                                                                        <  2:174216/226‑1 (MQ=255)                                                                                                                                   cGTCATTGAAAAAGGTGGCGAAT‑TTGTCGGTGAAGGTTACCCCCAA‑GGGGCGGGTAAACC‑CAATGCCGAAGTACCCGGTTGGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCt                                                                                               <  1:1186953/238‑1 (MQ=255)                                                                                                                                      cATTGAGAAAGAGGGAGAAT‑TTGTCGGTGCGGGTTCCCGCAAA‑GGTGGGGGTGACCC‑CCACGCCGAGGTCCACGCGTTGCGTAGGGC‑GGGGGAAAAGGC‑CAAGGGG‑GCAACCGCCTAGGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGc                                                                                           <  2:235245/239‑1 (MQ=255)                                                                                                                                         tGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGTAAAACC‑CAAAATG‑GCGACCGCCTATGT‑ACACCTCGAACCCTGTATCCCTCCTTGTTGCACGCCCCCGTGCTGTTACGCACCCCTCTCCCCCTGCGTAGCACGCGTGGTTGCCTCGATGCCCGAACCCAACCCCCAAGG‑CGCGGGGCgcg                                                                                        >  2:166777/1‑237 (MQ=255)                                                                                                                                          gTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAAC‑CAAAAGG‑GCGGCCGCCTTTGTCACA‑CTCGTAACCTTTAGCCATCATTGTCGACCGCCCCCCTGGTGTGACGCCCTCACCGCCGCTTGCGTAGCGCCGGGGGGTGCCCCGGTCAACGACCCTAACCCCCAGGG‑CGCggggggggg                                                                                       >  2:1305233/1‑234 (MQ=255)                                                                                                                                                     cGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCCCAGGG‑CGCTGGGCGTGGGCTTTACCg                                                                              >  1:538600/1‑237 (MQ=255)                                                                                                                                                      aaaa‑TTCGCGGTGAAGGTTACCACCAG‑CGGGCGGGGGAACC‑CAATGCCGAAGTACCCGCTTTGCGTAGGGG‑GGGTGAAAAGCC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTc                                                                            <  1:1305233/237‑1 (MQ=255)                                                                                                                                                                        ttACCACCAA‑CGGGCGGGTGAACC‑ACATGCCGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                      <  1:435474/180‑1 (MQ=255)                                                                                                                                                                        gtACCCCCAA‑CGGGGGGGGGAACC‑ACATGCCGAAGTACCCGGGTTGCGTATGGC‑GGGTAAAAAAGC‑AAAAGGT‑GCGACCGCCTATGTCACG‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAtt                                                           <  1:951299/237‑1 (MQ=255)                                                                                                                                                                                   ggcgggggTGACCC‑ACAGGCCGGAGTACCCGGGGGGCGGATGGG‑GGGTGAAAACGC‑CAGGGGT‑GCGCCCGCCTATGTCACC‑CTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                <  2:960969/234‑1 (MQ=255)                                                                                                                                                                                   cgggggggTAAACC‑CAAGGCCGAAGTCCCCGGGTGGCTTATGGG‑GGGTAAAAAACC‑CAAAGGG‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                <  1:975689/234‑1 (MQ=255)                                                                                                                                                                                   cgggCGGGTGAACC‑AAGTGGCGAATGACGGGTTTGGCTTAGGGC‑GGGTAAAAACGC‑AAAAGGT‑GCCACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTAGCATTGACGTCAGCCa                                                <  1:799241/236‑1 (MQ=255)                                                                                                                                                                                        gggTGAACC‑ACAGGCGGAA‑TTCCCGGGTTTGGTGAGGG‑GGGGAAAAAACC‑CAAAGGTGGCAACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGGCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                            <  2:538600/238‑1 (MQ=255)                                                                                                                                                                                         gggggaccccgttGCCAAGGTACACGGGTGGCGTTGGGG‑GGGTGAAAAACC‑CAAGGGG‑CCACCCGCCTTGTTCACC‑CGCGAACCCTGTGGCCATCAGGGTCGTACGCCACCGTGCTGTGGCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                            <  2:359011/226‑1 (MQ=255)                                                                                                                                                                                          gTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGC‑CAAAGGT‑GCGACCGCCCATGTCACA‑CTCGAACCCTGTAACCATCATGGTCGTACGCCACCGGGCTGTGACCCCCTCATCGCCCCCGGCGTAACGCGCGGGGTTGCCCCGGTGCCAGATCCTAACCCCCCGGTACCCGGGGCG‑GGGCTTTACCGTCTGCAACAGGCTGGCCTTGAGGGCAGCCACGGcccg                                         >  2:1117436/1‑237 (MQ=255)                                                                                                                                                                                                           aagaagacacgggCGTATGGC‑GGGTGAAAAAGC‑CAACGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                         <  1:1363463/227‑1 (MQ=255)                                                                                                                                                                                                              aaccgggTTGCATATGGC‑GGGT‑AAAAAGCACAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                       <  1:930248/232‑1 (MQ=255)                                                                                                                                                                                                                            tgggggggTAAAAAACC‑CAAAGTG‑CCCACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAgg         <  1:105846/234‑1 (MQ=255)                                                                                                                                                                                                                                   gTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                             >  1:752538/1‑117 (MQ=255)                                                                                                                                                                                                                                   gTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                             <  2:752538/117‑1 (MQ=255)                                                                                                                                                                                                                                   gTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCTGGCGTAGCGCGCGTGGGTGCCTCGCTGCAAGATCCTAACCCCCCGGT‑GGCTGGGCGGGGACTTTCCCGTCTGCACAAGGCGGGCCTTGACGTCCGCCCCGGCCTGAGGGTGGGTGGAGCCGGGGAAATGaaaatgtctctttta  >  2:1363463/1‑226 (MQ=255)                                                                                                                                                                                                                                               aaaGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCcatca                                                                                                                                                                                           >  2:247156/1‑44 (MQ=255)                                                                                                                                                                                                                                               aaaGGG‑GGGACCGCCTATGTCACA‑CTCGAACCCTGTAGCcatca                                                                                                                                                                                           <  1:247156/44‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  AGAATTTGGCGAAGAGATCGC‑GCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAA‑C‑GAGGACGTTTTACCACGCATCCCAACCCGAATGT‑CGGGTGCGTCATTGTCAAAGATGGCGAAA‑TTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACC‑ACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGC‑CAAAGGT‑GCGACCGCCTATGTCACA‑CTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGT‑CGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433415‑433871

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