Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,096,090	0	A	G	36.2%	72.7 / 13.8	47	I83V (ATT→GTT)	hisF	imidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base A (4/26);  new base G (17/0);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 2.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCTGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAAC  >  NC_000913/2095844‑2096319                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           aTGCTGGCAAATCGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTaa                                                                                                                                                                                                                                        <  1:187658/251‑1 (MQ=255)         aaaaCGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGCACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑AGCTCTgg                                                                                                                                                                                                                               >  2:200401/1‑251 (MQ=255)          aaaCGCATAATCCCATGTCTCGACGTTCGTGATGGGCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑a                                                                                                                                                                                                                             <  2:208994/251‑1 (MQ=255)              gCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GGCTCCGG‑a                                                                                                                                                                                                                             >  1:47569/1‑247 (MQ=255)              gCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGAt                                                                                                                                                                                                                         <  1:382663/251‑1 (MQ=255)               catcatcccaggttgcGACGTTCGTGATGGTCCGGTGGTGAAAGGC‑GTACAGTTTCACAACCATGAAATCAGTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCGCCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATg                                                                                                                                                                                                                        <  2:221220/236‑1 (MQ=255)                aTAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAA‑Gtctatgttgagagc                                                                                                                                                                                                                        >  2:208586/1‑240 (MQ=255)                  aaTCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGcc                                                                                                                                                                                                                      <  1:78962/250‑1 (MQ=255)                          gTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGATACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTTG‑aaaaaggccccaaaatt                                                                                                                                                                                                             >  1:305821/1‑236 (MQ=255)                             tcGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGATTAA‑GTCTCTGG‑GAGATGcccccgaaatattt                                                                                                                                                                                                          >  2:265556/1‑239 (MQ=255)                                 cGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCt                                                                                                                                                                                                   <  1:116391/251‑1 (MQ=255)                                 cGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCt                                                                                                                                                                                                   <  1:205137/251‑1 (MQ=255)                                      gttttggtgcggtggtGAAAGGCTGT‑CAGTTTCGCAACCAACAAATCATGGGCGATTTCGGGCCGCTGGCAACACGCTACGCTGAAGAAGGCGCTGACGAACTGGGGTTCTACGATATCACCGCTTCCAGCGGTGGCCGTGTGGTATATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACAATCCGTTTTGTGTGGCGGTTTGGATTAA‑GTCTCTGG‑AATATGCCG‑AG‑AAAATTCTTT‑CCTTTGgc                                                                                                                                                                                              <  2:173178/241‑1 (MQ=255)                                      gTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGgc                                                                                                                                                                                              <  2:86591/251‑1 (MQ=255)                                      gTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGATAA‑GGCTTTGG‑GAGAGGGCG‑CG‑CAAATTCTTTTCCTTTgg                                                                                                                                                                                               >  1:5834/1‑251 (MQ=255)                                           ggTCAGGTTGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATTATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTTAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCTGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCgg                                                                                                                                                                                          <  2:348852/250‑1 (MQ=255)                                                   ggtgAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTGA‑GACTCTTG‑GAGATGGCG‑CG‑CAAATTTTTT‑TCTTTGgccgggggaaaaatt                                                                                                                                                                                 >  2:225928/1‑238 (MQ=255)                                                      gAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTcc                                                                                                                                                                               <  1:201393/250‑1 (MQ=255)                                                       aaaGGC‑GAACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAAGT‑TTTGG‑AAGATTCCGGCG‑AAAATTTTTT‑CTTTTggtcgggtaaaaatttcat                                                                                                                                                                              >  1:46580/1‑234 (MQ=255)                                                         aatgcgtgctgcttCGCAACCATGAAATTATTGGCTATATCTTGCCGCTGGCAAAAGGCTACGCTGAAGAACGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGGTCGACATTCCGTTTTGTGTTGCGGTTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGTCTCGGATAAAATTTCCATc                                                                                                                                                                            <  2:34695/239‑1 (MQ=255)                                                                aCAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACtctc                                                                                                                                                                     <  1:15654/251‑1 (MQ=255)                                                                  aGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCt                                                                                                                                                                   <  2:239636/251‑1 (MQ=255)                                                                   gTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑GAAATGCCC‑CG‑AAAAATCTTT‑TCTTTTGCGCGGGTAAAATTTTCATCAACTATCCtt                                                                                                                                                                  >  1:192810/1‑250 (MQ=255)                                                                    tttCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGG‑AAGATGCCC‑CG‑AAAATTCTTT‑TCTTTTGCGCGGGTAAAATTTTCATCAACTCTCCTgc                                                                                                                                                                 >  1:36160/1‑251 (MQ=255)                                                                             cATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTA‑GGCTCTGG‑CAGAAGCCG‑CG‑GAAATTCTTT‑CCTTTGGTGCGGGTACAATTTTCATCAACTCACCTgacctgggcg                                                                                                                                                         >  2:104984/1‑241 (MQ=255)                                                                                 aaaTCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCa                                                                                                                                                    <  1:204621/251‑1 (MQ=255)                                                                                         ggCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTA‑GTGTCTGG‑AAGATGCCC‑CG‑AAAATTTTTT‑TCTTTTGGGCGGGTAAAAATTCCCTCAAACCCCCTGCGCTgggcgaccccacaataa                                                                                                                                             >  2:385360/1‑235 (MQ=255)                                                                                               atCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCg                                                                                                                                      <  1:28122/251‑1 (MQ=255)                                                                                                    gCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTg                                                                                                                                  <  2:306813/250‑1 (MQ=255)                                                                                                           gCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATAGCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAGCTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCg                                                                                                                          <  2:152086/251‑1 (MQ=255)                                                                                                                         tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTgg                                                                                                                                                                                               <  2:150981/168‑1 (MQ=255)                                                                                                                         tGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTgg                                                                                                                                                                                               >  1:150981/1‑168 (MQ=255)                                                                                                                          gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGG‑GAGATGGCG‑CGGAAAATT‑TTT‑CCTTTGGCGCGGGTAAAATTTTCCTCAACTCCCCCGCGCCGGCAGACCCCACACTAATTACTCCCCTGGCCGCTCGCCTTGGGGTGccggg                                                                                                           >  2:89749/1‑247 (MQ=255)                                                                                                                                                    cTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATAGACATTCCGTTTTGTGTGGCGGGTGGGGGTAA‑GTCTCTGG‑AAAATGGCG‑CG‑AAAAATTTTT‑CCTTTGGAGCGGGTAAAATTTTCAACAAATCTCCTGGGCCTGCAGACCCCACATTAATTACCCGCCTGGCCGCACGCTTTGGCGTGCAGGGTATTGTGGTGGGTATTGATAcccggg                                                                                 >  1:209926/1‑247 (MQ=255)                                                                                                                                                      aCGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTA‑GTCTCTGG‑GAGAAGCCC‑CG‑AAAATTTTTT‑CCTTTTGGGGGGATAAAAATTTCATCAACTCTCCCGCCCTGGCCGACCCCACACTAAATAATTCCCTGGCCGATCGCTTTTGCGGGCAGTGGCTTTTGGTTGGTATTggttccctgtac                                                                               >  2:238175/1‑240 (MQ=255)                                                                                                                                                        aatatCACCGCTTCCATCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGAGTAA‑GGTTCTGGGAAG‑TGACG‑CG‑AAAATTATTT‑TCTTTTGCGCGGGTACAATTTTCATCAAATCTCCTGCGCTTGCAGACCCAACATTAATTAATCGCCTGGCCGGTCGCGTTGGGGGGCGGTGACTTGTGGTCGGTATGGAtgaccgggtagt                                                                             >  2:264464/2‑240 (MQ=255)                                                                                                                                                                       aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAATCCGGTaa                                                              <  1:187075/251‑1 (MQ=255)                                                                                                                                                                                     ggTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCTG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATc                                                <  2:47569/251‑1 (MQ=255)                                                                                                                                                                                               aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAA‑GTCTCTGG‑GAGATGCCG‑CG‑GAAATTCTTT‑TCTTTTGCGCGGGTAAAATTTCCATTAACTCTCCTGCGCTGGCAGACCCAACCTTAATTACTCGCCTGGCCGCTCGCTTTGGCGTGCCGTGTATTTTGGGCGGGATTGATACCTGGGACGACGCCGAAACCCGGAAATATCATGTGAATCAATATAACgg                                      >  1:135277/1‑251 (MQ=255)                                                                                                                                                                                                       aTCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGG‑AAGATGCCC‑CG‑CAAATTCTTT‑TCTTTGGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCCACATTAATTACTCGCCCGGCCGATCGCTTTGGCGTGCGGTGTATTGTGGTCGGGATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGaa                               >  1:33144/1‑250 (MQ=255)                                                                                                                                                                                                           cgcgTGGCGGAAGTGTTCGACATTCCATTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCTCTGGCAGACCCATCATTAATTTCTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGGCGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGc                            <  1:149146/249‑1 (MQ=255)                                                                                                                                                                                                           cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGa                                                   <  1:25889/226‑1 (MQ=255)                                                                                                                                                                                                           cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAA‑GTCTCTGG‑GAGATGCCG‑CG‑AAAATTCTTT‑CCTTTTGCGGGGGTAAAAATTTCATCAACTCTCCCGCGCTGGCAGAACCAACATTAATTAATCGCCTGGCCGATCGCTTTGGCGTGCAGGGTATTGTGGGCGGTATTGATACCTGGTACGACGCCGAAAACGGGAAATATTATGTGaagatcggaagagcgtcggggagggg                          >  2:25889/1‑227 (MQ=255)                                                                                                                                                                                                             cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGCGGGGGTAA‑ATCCCTGG‑GCGAAGCCC‑CC‑GAAATTTTTT‑CCCTTTGCGCGGGGAAAAATTTCATCCACCCTCCCGCGCTCGCGGGCCCAAACTTTATTAATCGCCCTGCCGATCGCTTTTGGGTGCAGGGGATTGGTGTCGGGATTTATACATGGGACGACCCCGAAACCCGTCAATCTCATTTGAATCAATATACCCGCGATGGAAGcctgc                        >  2:128429/1‑248 (MQ=255)                                                                                                                                                                                                                    gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGt                 <  2:86289/251‑1 (MQ=255)                                                                                                                                                                                                                     aaGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTc                <  1:246541/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAAc  <  1:200401/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           ATGCTGGCAAAACGCATAATCCCATGTCTCGACGTTCGTGATGGTCAGGTGGTGAAAGGC‑GTACAGTTTCGCAACCATGAAATCATTGGCGATATCGTGCCGCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAA‑GTCTCTGG‑AAGATGCCG‑CG‑AAAATTCTTT‑CCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTGGTACGACGCCGAAACCGGTAAATATCATGTGAATCAATATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAAC  >  NC_000913/2095844‑2096319

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