Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	36.0%	55.3 / 17.8	77	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (7/41);  new base G (22/5);  total (31/46)
Fisher's exact test for biased strand distribution p-value = 1.21e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCA  >  NC_000913/433416‑433884                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               gaagtgggagAAGCGAGCGCGCCCCTGGAGGACTAAGCCGTGAAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGtc                                                                                                                                                                                                                                               <  2:893197/226‑2 (MQ=255)    tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAACCAAAAGGGCGAccc                                                                                                                                                                                                                                         >  1:1646990/1‑237 (MQ=255)                aTCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGATTATTACATGGCGCGGGCGATAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAGGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCac‑ac                                                                                                                                                                                                                           <  1:1045321/238‑1 (MQ=255)                  gggccgcCTGGAGGACTAAGCCGTGCAGGAGGATTATTACAGGGCGCGGGCGCTAGAGCT‑GGCGCAACGAGGGCGTTTTACCACG‑CACCCCAACCCGAA‑TGGCGGGTCCGTCATTGACAAAGAGGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACCCGCGTTGCGTAGGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCg                                                                                                                                                                                                                        <  2:231026/235‑1 (MQ=255)                        ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAACCAAAAGGGCGACCG‑CCTATGTCAC‑ACTCGAAcccc                                                                                                                                                                                                                  >  1:1321003/1‑238 (MQ=255)                           gtaggaCTAGGCCGTGCAGTACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGGTTTACCACG‑CACCCCAACCCGAA‑TGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGTTGAACCACATGCCGAAGTACCCGCGTTGCGTATGG‑C‑GGGTAAAAAAGCAAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTa                                                                                                                                                                                                               <  2:2092323/237‑1 (MQ=255)                               gaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGAG‑GGGTGAAAAA‑CCAAAAGGGCGACCC‑CCCATGTCACAACCCGAAcctggaggcc                                                                                                                                                                                                            >  1:2055703/1‑231 (MQ=255)                                  cacgcagTGCAGGGCTGGTATTCGAGGGCGAGGGCG‑TAAAACTGGGCGCACCG‑GGTCCTTTTTCCCCGCCATCCAACCCCAAA‑TTGCGGGGGGGTAATTGCAAAAGAGGGCGAAATGGTGGGTGAGGGTTCCCCCCACGGGGCGGGTGACCCCCATGCCAAAGTCCCCGGGTTGCGTATGG‑G‑GGGTAAAAAAGCCAAGGGTGGGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCcatca                                                                                                                                                                                                        <  2:667163/233‑1 (MQ=255)                                   aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAACCAAAGGAGCGGCCG‑CCTATGCCCC‑CCTCCAACCCTCTAGCcctctt                                                                                                                                                                                                       >  1:223427/1‑237 (MQ=255)                                         tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGGGCGACCG‑CCCATGTCAC‑ACTCGAAACCTGTAGCCCACATGGGCGt                                                                                                                                                                                                 >  1:407712/1‑239 (MQ=255)                                         tGCAGGACGAGTATTACATGGCGCGGGCCCTAAAGCG‑GGCGCACCGAGGACGTTTTACCACG‑CATCCAACCCCAAA‑TTTCGGGTGCGTCATGGCCAAAGATGGCGAAATGGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                     <  2:1321003/235‑1 (MQ=255)                                                tgagttgtACTTGGCCGGGGCGTTAACGCG‑GGCGCAAGGAGGACGTTTACCCTCG‑CATCCCAACCCGAG‑TGGCGGGTGCTTCGTTGCCAAAGATGGCAGAATTTTCGGTGGAGTTGCCCACCAGCGTGGGGTTGACCCACGTGCCGACGTCCCCCTTTTGCGTATGG‑C‑GGGTGAAAAGCCCAAGGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                           <  2:2055703/231‑1 (MQ=255)                                                cGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAACCAAAGGTGGGACCG‑CCTCTGTCAC‑ACCCCAACCCTGGAGCCATCATGGGTGTAAGCCa                                                                                                                                                                                           >  1:1595137/1‑238 (MQ=255)                                                   gTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAACCAAAAGGGGGACCG‑CCTCTGTCAC‑ACCCCAACCCTGTAGCCAACATgg                                                                                                                                                                                                     >  1:42220/1‑225 (MQ=255)                                                   gTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAACCAAAAGGGGGACCC‑CCCATGTCAC‑AACCCAACCCCGTAGCCAACATGGGTGTACCCCCCCgg                                                                                                                                                                                       >  1:183183/1‑238 (MQ=255)                                                   gTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGCAACCATCATGGTCGTACGCCACCGt                                                                                                                                                                                       >  1:1072880/1‑239 (MQ=255)                                                   gTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑AATCGAACCCTGTACCCCACCTGGTCGTACGCCAcccg                                                                                                                                                                                       >  1:42878/1‑237 (MQ=255)                                                    tcttactgggCGCGGGCTTTAAGGCT‑GGCGCACCGAGGACGTTTTACTGCG‑C‑TCCACACCCGAATTGCCGGGGCCTCCTTTGTCAAAGGGGGCAAATGTGTCGTGAAAGTTTCCCCCCACCGTGCGGGCGCCCCCCTGGCCGACGCCCCCGCGTGGTGTTGGG‑C‑GGGAAAAAAAGCCAAGGGTGCGCCCG‑CCTATGTCAC‑CCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                      <  2:1813020/231‑1 (MQ=255)                                                     aTTACATGGCGGGGGCGATAAAGCT‑GGCGCAACGAGGACGTTTTCCCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                      <  1:1796735/238‑1 (MQ=255)                                                        aCATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAACCCAAAGGGCGACCC‑GCTATGGCAC‑CACCGAAcc                                                                                                                                                                                                                    >  1:1988084/1‑205 (MQ=255)                                                              cgcgGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑CGGGGGGAAAAAGCC‑AACGTGGGGCCGCCCCATGCCA‑‑ACTCGAAACCCGTAGCCATCCCCGGCCGACGCCACCCTGCtgttcagc                                                                                                                                                                             >  2:699433/1‑233 (MQ=255)                                                              cgcgGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGGAAAAGCCAAAGGGGCGACCG‑CCCCTGTCAC‑ACTCGCACCCCGTAACCATCCAGGTAGGTAGCCACCGAGCTTGGACGc                                                                                                                                                                             >  2:1256657/1‑238 (MQ=255)                                                               cgggcgtcaaaaGGT‑GGCCCAACGATGCCGTTGTCCCACG‑CATCCCAACCCGAT‑GGTCGGGTGCGTCATGGCCAAAGATGGCAAAATGGACGGGGAGGGTCACCACCAACGTGCGGGTGACCCACATGCCGAATTACCCTCGTTCCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCTACCG‑CCTATGTCAC‑ACTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCAc                                                                                                                                                                           <  2:407712/230‑1 (MQ=255)                                                                   gcgcTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACCCCAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTcc                                                                                                                                                                        >  1:1946252/1‑237 (MQ=255)                                                                             t‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAACCAAAGGGGCGAACG‑CCTAGTTCAC‑CCTCCAACCCCGTAg                                                                                                                                                                                                              >  1:1523039/1‑190 (MQ=255)                                                                                                         cATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:1274966/239‑1 (MQ=255)                                                                                                          gtcacAACCCGAA‑TGTCGGGTGCGTCGTTGGAAAAGAGTGAGAAAGGTGCGGGGAAGTTGCCCACCAACGTGCGGGTGAACCACATGCCGAAGAGCACGCGTGGGGTATGG‑G‑GGGTAAAAACACAAGGGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  2:115077/234‑1 (MQ=255)                                                                                                                      a‑TGGCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGGCGGTGAAGGTTCCCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTAGGG‑G‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGaa                                                                                                                                                                                                                      <  1:2024093/143‑1 (MQ=255)                                                                                                                                   tCATTGGAAAAGAGGCCGAATTTGTGGGTGAAGGTTACCACCAACGGGGGGGTGACCCACATGCCGAAGTACCCGGTTTGCTTATGG‑G‑GGGTGAAAAAGCAAAAGGTGCGACCG‑CCTATGTCAC‑ACTc                                                                                                                                                                                                                         <  2:807205/128‑1 (MQ=255)                                                                                                                                    taTTGTCAAAGATGGCGAAATTGTGGGGGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                         <  1:1256657/238‑1 (MQ=255)                                                                                                                                       tGTCAAAGATGGCGAAATTGTCGGTGAAGTTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGc                                                                                                                                                                                                             <  1:1194100/136‑1 (MQ=255)                                                                                                                                       gcataaatggggCGAAGTTGGCGGGGATGTTGACCCCCACCGGGCGGTGTAACCCCTGGCCGAAGTCCCCGCGTTGCGTGTGG‑G‑GGGAGAAAACGCAAAAGGTGCCACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                       <  2:1119474/228‑1 (MQ=255)                                                                                                                                         tCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑G‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                      <  1:2030007/157‑1 (MQ=255)                                                                                                                                         acAAGGATGGAGAAGTGGCGGGTGAAGTTTCCCCCCACGGGGCGGGGAAACCGCAGGCCGAAGTCCCCGCTTTGCGAAGGG‑C‑GGGTGAAAACGCCAAGGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTc                                                                                                                <  2:1390484/226‑1 (MQ=255)                                                                                                                                             aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAAGGGCGACCG‑CCTCTGTCAC‑ACTCGAACCCTGTa                                                                                                                                                                                                               >  1:807205/1‑128 (MQ=255)                                                                                                                                               aTGGCGAATTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGtt                                                                                                                                              <  1:661711/191‑1 (MQ=255)                                                                                                                                                   cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCCAAGGGGCGACCC‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGCCGCTGGGCggg                                                                                                      >  1:1390484/1‑225 (MQ=255)                                                                                                                                                         tGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGAATGG‑C‑GGGTGGAAAAGCCAAAGGCGGCACCG‑CCGATTTCTC‑ACTCGGCCCCTGGAACCCTCCTGGTCGTAGCTCAACGTGCTGTGACCCCATCACCCCCACCTGTGTACCGCGGCTGGCTGCCTCCCAGCAAGAACCTCACCCCCCGGGCGCCGGGCGTGGGA‑CTTATCGTCTGGAAc                                                                                   >  2:396180/1‑239 (MQ=255)                                                                                                                                                            cGGTGAGGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACCCGCGTTGCGTAGGG‑C‑GGGTAAAAAAGCCAAAGGTGCGCCCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAgg                                                                                <  1:699433/239‑1 (MQ=255)                                                                                                                                                            cGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAACCAAAGGGGCGACCG‑CCTCTGTCAC‑CCTCGAACCCCCTAGCCACCATGGTCGTACGCCACCCCGCTGTGACGCACCCCTCGCCGCTGGCGTAGCGGGCGGGGGTGCCTCGCTGCAAGAGCCCAACCCGCAGGTGCCTGGGCGCGGGC‑TTTACCGTCTGCAACac                                                                                 >  2:287959/1‑237 (MQ=255)                                                                                                                                                             ggTGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTAGGG‑C‑GGGTGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTa                                                                                                                                                                                                               <  2:1211033/112‑1 (MQ=255)                                                                                                                                                              gggAAGGTTCCCCCCAACGGGGGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTAGGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:318831/185‑1 (MQ=255)                                                                                                                                                              gggAAGGTTACCCCCAAGGGGCGGGTGACCCACATGCCGAAGTACACGCTTTGCGTAGGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGc                                                                               <  1:287959/236‑1 (MQ=255)                                                                                                                                                                    gTTCCCACCAACGTGCGGGTAAACCACAGGCCGAAGTCCACGCGTTGCGTAGGG‑C‑GGGTAAAAAAGCAAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                                                            <  2:592353/140‑1 (MQ=255)                                                                                                                                                                    gTTACCACAAAGGGGGGGGGGAAACACATGCCGAAGTCCAGGAGTGGCGTGTGG‑G‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                              <  1:983603/202‑1 (MQ=255)                                                                                                                                                                     ttACCACCAACGTGCGGGTAACCCCCATGCCGAAGTACCCGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCAt                                                                         <  1:1051498/237‑1 (MQ=255)                                                                                                                                                                       accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGTGGGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                     >  1:1211033/1‑112 (MQ=255)                                                                                                                                                                        ccaccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCCAAGGGGGGACCG‑CCCATGTCCC‑ACTCCAACCCCGTAGGCCTCATGGTCGTACGCCACCGTGCTGTGACGCCCTCCTCCCCCCCGGCGTCGCGCGCGCGGGTGCCTCGATGCAAGATCCTAAc                                                                                                                         >  2:318831/1‑187 (MQ=255)                                                                                                                                                                            caACGGGCGGGTGAACCCCATGCCGAATTCCCCGCGTTGCGTAGGG‑G‑GGTTAAAAAACCAAAAGGTGCGCCCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                         <  1:963528/199‑1 (MQ=255)                                                                                                                                                                              aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATcgccg                                                                                                                                                                  >  1:592353/1‑140 (MQ=255)                                                                                                                                                                                tggggggTGAACCACTAGCCGAAGTACTCGGTTGGCGTAGGG‑G‑GGGTGAAAAAGCCAAGGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGcc                                                              <  1:741296/233‑1 (MQ=255)                                                                                                                                                                                gTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGG‑G‑GGGTAAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGa                                                                                                                                                                                                                       <  2:391437/85‑1 (MQ=255)                                                                                                                                                                                 tGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGG‑G‑GGGTGAAAAAGCAAAAGGTGCGTCCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGc                                                                                                                                                                                     <  1:218537/118‑1 (MQ=255)                                                                                                                                                                                  gggggTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                           <  1:708796/236‑1 (MQ=255)                                                                                                                                                                                      gTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTg                                                                             <  1:1034858/216‑1 (MQ=255)                                                                                                                                                                                      gTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCCC‑CCCCGAACCCCGTGGCCATCATGGGCGGACGCCACCGTGCTGTGACGCACCCATCGCCCCCGGCGTCGCGCGCGTGGTTGCCTTCGTTCAAGATCCTTACCCGCAGGGCGCTGGGCGGGGGC‑GTTa                                                                                              >  2:963528/1‑199 (MQ=255)                                                                                                                                                                                       tGAACCACATGCCGAAGTCCCCGCGTTGCTTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtg                                                      <  1:1453170/238‑1 (MQ=255)                                                                                                                                                                                         aaCCACATGCCGAAGTACGCGCGGTGGGTGTGG‑G‑GGGTGAAAAAGCAAAAGGGGCGACCG‑CCTATGTCAC‑TCTCGAACCCGGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGc                                                                               <  1:1891726/211‑1 (MQ=255)                                                                                                                                                                                          aCCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGc                                                                                                                                                                                                             >  1:391437/1‑85 (MQ=255)                                                                                                                                                                                           ccACATGCCGAAGTACACGCGTTGCGTTTGG‑C‑GGGTAAAAAAGACAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCcatcattg                                                                                                                                                                                                     <  2:1723240/92‑3 (MQ=255)                                                                                                                                                                                            cacaGGCCGAAGTACCCGCGTTGCGTAGGG‑C‑GGGTAAAAAACCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                                                                                                                                                                              <  2:1404617/114‑1 (MQ=255)                                                                                                                                                                                                tGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGGCAAAGGGGCGACCG‑CCTATGTc                                                                                                                                                                                                                                >  1:919679/1‑61 (MQ=255)                                                                                                                                                                                                tGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGGGGAAAAACCAAAAGGGGGGCCC‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCAAGGGCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGGCTCGGTGCCAGATCCTAAAACCCCGGGCGCCGGGGGGGGGCTTTTACCG‑CTGCCACAGGCTCGGCTTGACGt                                                                   >  2:1034858/1‑216 (MQ=255)                                                                                                                                                                                                   cgattaacacGCTTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGGCCAACCG‑CCTATGTCAC‑CCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg                                         <  1:1340684/233‑1 (MQ=255)                                                                                                                                                                                                   caagtcacgcgcgTTGTGTAGGG‑C‑GGGTGAAAAACCCAAGGGTCGTCCCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg                                         <  1:668234/230‑1 (MQ=255)                                                                                                                                                                                                     aaGTCCCCGCGTTGCGTAGGG‑C‑GGGTAAAAAAGCAAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                         <  1:396953/221‑1 (MQ=255)                                                                                                                                                                                                     aaGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAACCCAAGGGGCGACCG‑CCCATGTCAC‑ACCCGAACCCTGTAGCCATCATGGGCGGACGCCa                                                                                                                                                                                           >  1:1723240/1‑92 (MQ=255)                                                                                                                                                                                                      aGTACACGCGTTGCGTATGG‑C‑GGGGGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACCCGAAACCTGTAGCCATCATGGTCGTACGCCACCCGGCTGTGACGCACTCATcgc                                                                                                                                                                    >  1:1404617/1‑114 (MQ=255)                                                                                                                                                                                                       gTACAGGCGTTGCGTAGGG‑C‑GGGTAAAAAAGCCAAGGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCcatca                                                                                                                                                                                                        <  1:470955/77‑1 (MQ=255)                                                                                                                                                                                                        tccacGCTTTGGGTTTGG‑C‑GGGTGAAATAGCCAAAGGTGCGACCG‑CCTATGTCTC‑CCTCGAACCCTGTATCcatcat                                                                                                                                                                                                       <  1:1450691/75‑1 (MQ=255)                                                                                                                                                                                                           cgtcgTTGGTTAGGG‑C‑GGTTAAAAAATCCAAAGGTCCGCCCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                 <  1:396180/236‑1 (MQ=255)                                                                                                                                                                                                               gTTGCGTATGG‑C‑GGGTGAAAAAACCAAAGGTGCGACCG‑CCCATGTCAC‑ACCCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCCCCGGCGGAGCGCGCGTGGGTGCCTCGATGCCAGATCCTAACCCGCCGGGCGCTGGGCGGGGGC‑CTTTCCGTCCGCAACCGGCTGGCACTGGCCGCAGGCACGGCCtggtgatga                                               >  2:396953/1‑221 (MQ=255)                                                                                                                                                                                                                   cGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:696489/134‑1 (MQ=255)                                                                                                                                                                                                                     tATGG‑G‑GGGTAAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                              <  1:842392/153‑1 (MQ=255)                                                                                                                                                                                                                               gTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTa                                                                                                                                                                                                               <  2:1061413/48‑1 (MQ=255)                                                                                                                                                                                                                               gTGAAAAAGCCAAAGGGGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTa                                                                                                                                                                                                               >  1:1061413/1‑48 (MQ=255)                                                                                                                                                                                                                                         cAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGCGGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAATCGGATg      >  2:1837795/1‑238 (MQ=255)                                                                                                                                                                                                                                             ggTGCGCCCG‑CCTATGTCAC‑ACTCGAACCCTGTCGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCa  <  2:174852/238‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               GAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCT‑GGCGCAACGAGGACGTTTTACCACG‑CATCCCAACCCGAA‑TGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGG‑C‑GGGTGAAAAAGCCAAAGGTGCGACCG‑CCTATGTCAC‑ACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGAC‑TTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCA  >  NC_000913/433416‑433884

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