Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,634	0	T	G	36.2%	52.8 / 12.3	60	G60G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (7/30);  new base G (17/4);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.33e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATA  >  NC_000913/433411‑433860                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     tCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGaa                                                                                                                                                                                                                                   <  1:866808/227‑1 (MQ=255)           tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGGATGGCGGGGGGAAAAACCCaaa                                                                                                                                                                                                                          >  1:895904/1‑227 (MQ=255)           tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCaaaaag                                                                                                                                                                                                                       >  1:1838739/1‑226 (MQ=255)               aaagaacgcgcgcgCTGGGAGGACTAAGCCGTGAAGAGCGAGTTTT‑CCATGCCGCGGGCGAA‑AAAGCTGTCGCAACGAGGACTTTTTACCCCGCATCCCACCCCAAATGCCGGGGGCG‑TAATGGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTCCCACCACCGGGCGGGTGACCCACAGGCGGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAGGGT‑GCGa                                                                                                                                                                                                                  <  2:1757200/221‑1 (MQ=255)                          gcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGaa                                                                                                                                                                                                                                   <  2:1834829/202‑1 (MQ=255)                          gcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGaa                                                                                                                                                                                                                                   >  1:1837955/1‑202 (MQ=255)                             ccTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGT‑GAGACCGCCt                                                                                                                                                                                                            >  1:1931443/1‑221 (MQ=255)                                 gaggcttacgCCGTGCAGGACG‑GTATTGACATGGCCGGGGCGCT‑AAAGCGGGCCAAACGAGGACGTTTTACCACGCATCCAAACCGGATGGGCGGGTGGG‑TCATTGTAAAAGATGGCAAAATGGTC‑GTGGAAGTTTCCCCCCAACGTGCGGGTGAACCAAAGGCCGAAGTACCCGCGTTGCGTAGGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                <  2:1835612/220‑1 (MQ=255)                                  gggccTAAGCCGTGCAGGAGGAGTTTT‑ACAGGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGATGTTTTACCACGCATCCAACCCCGAATGTCGGGTGCG‑TCTTTGTCAAAGAGGGCAAAATTGTC‑GGTGAGGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTAGGGG‑GGGTAAAAAAGCCAAGGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                <  2:1890195/224‑1 (MQ=255)                                   ggaCTACGCGGGGCAGGAAGATTATT‑CAATGGCGCGGGCGAT‑AAAGCTGGCGCAACGAGGGATTTTTCCCCCGCACCCCAGCCGGAATGTCGGAGGCG‑TCATTGTCAAAGGTGGC‑AAAATGTCGGGTGAAGGTTCCCACCAGCGTGGGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGG‑GGGTAAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                              <  2:894341/229‑1 (MQ=255)                                        aaGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGG‑GCGGCCGCCCCTGTCCCACTCCAACCCTg                                                                                                                                                                                         >  1:1022155/1‑229 (MQ=255)                                           ccgagcggggagAGTTTT‑ACA‑GGGGCGGGCCCTAAAAGCGGGCGAAACGGGCACGTTTTACCACGCACCCAAACCCAAATGCGGGGTGCT‑TCATTGTCAAAGATGGCGAAATGGTG‑GGTAAAGTTTACCACCACGGTGCGGGTGAACCAAAGGCCGAAGTACCCGCGTTGCGTAGGGC‑GGGTAAAAAACCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                       <  2:634126/217‑1 (MQ=255)                                              tGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAAACAAAGAT‑GCGa                                                                                                                                                                                                                  >  1:654441/1‑198 (MQ=255)                                              ggCAGGACGAGTTTG‑ACAGGGCGGGGGCGGT‑AAAGCGGGCGAAACGAGGAGGTTTTCCCCCGCATCCCAACCCTAATGTGGGGTGCG‑TAATTGTAAAAGATGGAGAAATTGCC‑GGTGAGGGTTACCCCCAAGGTGGGGTTAAACCAAATGCCGAAGTACACGCGTTGCGTATGGG‑GGGTGAAAAAGCCAAAGGT‑GCgg                                                                                                                                                                                                                  <  2:653289/197‑2 (MQ=255)                                                 agaaCGTGTATT‑CCATGGCGCGGTCGCT‑AAAGTTGGCGCAAGAAGGACGTTTTCCCACGCATCCCAACCCGATTGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGTTTACCACCACCGTGGGGGTGAACCACAGGCCGAAGTCCACGCGTGGCTTATGGC‑GGGTAAAAAACCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                <  2:1196728/226‑1 (MQ=255)                                                      ttgTATC‑ACAGGGCGGGGGGGCA‑AAAGCGGGCGCACGGAGGCCGTTTTCCCACGCCTCCAACCCAGAATGCCGGGTGCG‑TCATGTGCAAAGAGGGCGAAATTGTC‑GGGGAGGGTTACCACCACCGGGGGGGTGAACCACAGGCGGAGGTACCCGCGTTGCGTAGGGG‑GGGTAAAAAAGCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                           <  2:726838/227‑1 (MQ=255)                                                        gTATT‑CCAGGGCGCGGGCGCA‑AAGGGTGGCCGAACGAGGCTGTTTTCCCACGCACCCCACCCCGAATGTCGGGGGCG‑TCATTGTAAAAGAGGGCGAAATTGTC‑GGTGAAGGTTCCCACCACCGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGG‑GGGTAAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                            <  2:1059573/210‑1 (MQ=255)                                                        gTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAACCAAAAGG‑GCGACCCCCTAGTGCCCACTCCAAcc                                                                                                                                                                                            >  1:1061441/1‑210 (MQ=255)                                                          tgtcaCAGGGCGCGGACGAT‑AAGGTGGGCGCAACGAGGCATTTTTCCCCCCCCCCCCAGCCCGAATG‑CGGGGGCGTTCAGTGGAAAAGAGGGGGAAATTGTG‑GGTAAAGGTTCCCCCCAGGGGGCGGGAGAACCCAATGCCGAAGTACAGGCGTTGCGTAGGGC‑GGGTAAAAAACCCAAGGGT‑GCAACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                             <  2:5082/220‑1 (MQ=255)                                                                   gcgcGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGCGGAAAAGCCCAAGGT‑GCGGCCCGCCATGTCACTCTCTGACCCCGTCACCACCATAGTccgtc                                                                                                                                                                       >  2:549868/1‑217 (MQ=255)                                                                    cgcgGGGGCT‑AAGGCGGGGGCAAGGGGGACGTTTTACCACGCATCCCACCCCGAATGTGGGGTGCG‑TAATTGTAAAAGATGGCAAAATTGTC‑GGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCGGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGa                                                                                                                                                                                               <  2:372333/196‑1 (MQ=255)                                                                    cgcgGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGG‑GCGACCCCCCATGTCacacccca                                                                                                                                                                                               >  1:373006/1‑192 (MQ=255)                                                                                                                 ccaaccccGAAGGGCGGGTGCG‑TCATTGTCAAAAGGGGCGAAATTGGC‑GGTGAGGGTTACCACCAACGTGCGGGTAACCCACATGCCGAAGAACACGCGTTGCGTAGGGC‑GGGTAAAAAAGCAAAGGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCt                                                                                                                 <  2:690712/224‑1 (MQ=255)                                                                                                                          aaggTCGGGTGCG‑TCATTGTAAAAGAGGGCGAAATTGTC‑GGTGAAGTTTACCACCACCGTGCGGGTAAACCACATGCCGAAGTCCCCGCTTTGCGTAGGGG‑GGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                <  2:1797133/139‑1 (MQ=255)                                                                                                                          aaTGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGG‑GCGACCCCCCATTTCACAATCg                                                                                                                                                                                                >  1:1800190/1‑142 (MQ=255)                                                                                                                                gggTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACtt                                                                                                                                                                                                 >  1:689956/1‑134 (MQ=255)                                                                                                                                gggTGCG‑TCATTGTAAAAGATGGGGAAATTGTC‑GGGGAAGGTTACCACCAAGGTGGGGGTGACCCACATGCCGAAGTACACGGGTTGCGTATGGC‑GGGTGAAAAAGCCAAGGGT‑GCGACCGCCTATGTCACACTc                                                                                                                                                                                                 <  2:688745/135‑1 (MQ=255)                                                                                                                                 ggTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAACCCAAGGG‑GCGCCCCCCTATGTCCCCATCGAACCCTTGAGCCATCATTGCCGTCCGCCCCCGTGCTGTGACGCACCCATCGCCGCTGGCGGAGCGCGCGCGGGTG‑CCTCCATGCAAGACCCTaa                                                                                                 >  2:1254233/1‑229 (MQ=255)                                                                                                                                     cG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAGCCAAAGGG‑GGGAACCCCCATGTCACACTCGGACCCCGTa                                                                                                                                                                                       >  1:40997/1‑140 (MQ=255)                                                                                                                                     cG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCCAAGGT‑GCGAACCCCTATGACACACTCGGACCCTGGAgaccccat                                                                                                                                                                               >  1:1236294/1‑141 (MQ=255)                                                                                                                                     cG‑TCATTGTAAAAGATGGCAAATTGGTC‑GGGGAAGGTTACCACCAACGTGGGGGTGAACCACATGCCGAAGTCCCCGTGTTGCGTATGGC‑GGGGAAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                               <  2:1234145/148‑1 (MQ=255)                                                                                                                                        tCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCa                                                                                           <  1:134101/229‑1 (MQ=255)                                                                                                                                                   ggtGGCGAAATTGGC‑GGTAAAGGTTCCCCCCAACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                              <  1:778339/150‑1 (MQ=255)                                                                                                                                                      ggCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                             >  1:393296/1‑134 (MQ=255)                                                                                                                                                      ggCAAAATTGTC‑GGTGAGGTTTCCCACCAGCGGGGGGGTAAACCACAGGCCAAATTCCACGCGTTGCGTTGGGG‑GGGTGAAAAAGCCAAAGGT‑GCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                             <  2:392585/134‑1 (MQ=255)                                                                                                                                                        cGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGACCAAGGG‑GGGAACGCCCATGTCACACCCGAACCCCGTAGCCACCCTGGCCGGACCCCACCGGGCTGTGACGCACCCATCGCCCCTGGCGGAGCGCCCGCGGGTG‑CCCCGAGTGAAGATCCCAACCCCCCGGTCgacgtgggtgga                                                                           >  2:1365718/1‑218 (MQ=255)                                                                                                                                                        cGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCAAGGTCGGACGCCACCCGGCTGTGACg                                                                                                                                                      >  1:405879/1‑155 (MQ=255)                                                                                                                                                        agAAATTGTC‑GGTGAAGTTTCCCACCACCGTGCGGGTAAACCACAGGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                                                                                                                                                      <  2:405147/154‑1 (MQ=255)                                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGT‑GCGACCGCCCATGGCACACTCGAACCCTGTAGCCATCATGGGCGTACGccc                                                                                                                                                                   >  1:1274973/1‑114 (MQ=255)                                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGG‑GCGACCGCCCATGTCACACCCGAACCCCGTAGCCCTCATGGTCGGACGCCCCCCTGCCGTTACCCACCCACCGCCCCTGGCGTAGCCCGCCGGGTTGCCCGCGTTCAAAGTCCCCAACCC‑CCGGGGGCTGGGCGCGGGCTTTACCGCCCGCAAcctgcgtgcta                                                 >  2:1616900/1‑218 (MQ=255)                                                                                                                                                                                    aCGGGCGGGTAACCCACATCCCGAAGTACACGCGTTGCTTAGGGC‑GGGTAAAAAACCAAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                   <  2:1272757/115‑1 (MQ=255)                                                                                                                                                                                       tGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGaa                                                                                                                                                                                                                                   >  2:472997/1‑49 (MQ=255)                                                                                                                                                                                       tGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGaa                                                                                                                                                                                                                                   <  1:473832/49‑1 (MQ=255)                                                                                                                                                                                               acccACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTAAAAAAGCCAAAGGT‑GCACCCGCCTATGTCACACTCGAATCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                     <  1:1727341/227‑1 (MQ=255)                                                                                                                                                                                               aaCCACGTGCGGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAGGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                     <  1:1475723/229‑1 (MQ=255)                                                                                                                                                                                                   ccaTGCCGAAGTCCCCGCGTGGCGTAGGGG‑GGGGTAAAAACCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                     <  1:1256415/224‑1 (MQ=255)                                                                                                                                                                                                   acaTGCCGAGGTACCCGCGTTGCTTAGGGC‑GGGTGAAAAAGCAAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                     <  1:40124/225‑1 (MQ=255)                                                                                                                                                                                                              tACCCGCTTTGCGTATGGG‑GGTTAAAAAACCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCt                                                                                                                 <  1:1030688/138‑1 (MQ=255)                                                                                                                                                                                                                      ttGCGTTGGGC‑GGGTAAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                   <  1:630278/208‑1 (MQ=255)                                                                                                                                                                                                                      ttGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGcc                                                                                                                                                                                                             >  1:615564/1‑39 (MQ=255)                                                                                                                                                                                                                      ttGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGcc                                                                                                                                                                                                             <  2:614486/39‑1 (MQ=255)                                                                                                                                                                                                                        gCGTATGGC‑GGGTGAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                              >  1:379373/1‑52 (MQ=255)                                                                                                                                                                                                                        gCGTATGGC‑GGGGGAAAAAACCCAAAGG‑GCGACCGCCTATGTCACACCCGAACCCCGTAACCATCAAGGTCGTACGCCACCGTGCTGTGACACACACATCCCCCCTGGCGTAGCGCGAGTGGGTG‑CCCCGGTGCAAGAACCTAACCCACACGGCGCCGGGCGTGGGGTTTTGCGTCTGGAACAGGCTGGGCTTGACCTCAGCCACCGCCTGGTTGTGAGTGAAGCCgcg            >  2:698928/1‑227 (MQ=255)                                                                                                                                                                                                                        gCGTATGGC‑GGGGGAAAAAACCAAAAGT‑GCGACCCCCCATGTCACCCCCGAACCCCGTAGCCATCATGGTCGTACGCCACCCTGCTGTGGCGCACCCCTCGCCGCCGGCGTAGCGCGCCTGGTTG‑CCTCGATGCAAGATCCTAACCCCCAGGGCGCTGGGCGTGGGCTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGGGGAGCCgcg            >  2:438827/1‑227 (MQ=255)                                                                                                                                                                                                                         cGTATGGG‑GGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg            <  1:439610/228‑1 (MQ=255)                                                                                                                                                                                                                           cttagggcggTTAAAAAGCCCAAGGTGGCAGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg            <  1:1725523/220‑1 (MQ=255)                                                                                                                                                                                                                            aTGGC‑GGGAGAAAGGGCCAAGGGT‑CCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt        <  1:1101205/229‑1 (MQ=255)                                                                                                                                                                                                                                ggggggGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCa                                                                                                          <  1:335093/126‑1 (MQ=255)                                                                                                                                                                                                                                    gTGAAAAAGCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                              <  1:514265/152‑1 (MQ=255)                                                                                                                                                                                                                                    gTGAAAAAGCCAAAGGG‑GCGACCCCCTATTTCCCACCCCAACCCCGTAGCCATCATGGTCGTACGCCACCCTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCCCGATGCAAAATCCTAACCCCCAGGTCGCTGGGCGGGGGCCTTTCCCTCTGCAACAGGCTGGCATTGACGTCAGCCCACGCCCGAAGATGAATGAAACCGAGCAATTGaaaa  >  2:905191/1‑226 (MQ=255)                                                                                                                                                                                                                                    gTGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                            <  1:311635/170‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     TCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATA  >  NC_000913/433411‑433860

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