Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,388	0	C	T	33.3%	43.8 / 11.4	42	A27T (GCA→ACA)	focA	formate channel FocA
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base T (1/12);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 6.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAAT  >  NC_000913/954138‑954615                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     aaCTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTACAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGc                                                                                                                                                                                                                                      >  1:501864/1‑251 (MQ=255)      aCCCCAGGTGATGCGCCCACTCGCCTTAGCGACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTtat                                                                                                                                                                                                                                 <  1:551598/251‑1 (MQ=255)      aCCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGGCTTAAGCGGAT‑GTTTCGTTGCTTtat                                                                                                                                                                                                                                 >  1:247611/1‑251 (MQ=255)       ccccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTtata                                                                                                                                                                                                                                <  1:370670/251‑1 (MQ=255)               gATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACAccc                                                                                                                                                                                                                         <  1:177934/250‑1 (MQ=255)                          tCGCCTTAGCAACAACAATCAACACGGAGGAAATACAGAGATCGGCACCGCAGACAACACAAAGAATCAGCCCCCGAAAGAAGCAAATGCCGCAAACCAGTTTTGCCATGCCGAAGGG‑AATTGTGCC‑TGTG‑CCAGTGGATGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAACCACC‑GG‑CGGTAATAGCAAGATAGAAAGTCATACGCGGAT‑GTTTTGTTGATTTATAGAAAACCGCCTCTTCGGcc                                                                                                                                                                                                             >  1:517354/1‑251 (MQ=255)                           cgcctacaaaaaaacaaacaACACGGGGAAAATAAAGAGAGCGGCCCCGCAGACAACACAAAAAACCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTGTTGCCATGCC‑AAGGGTCATGGGGCC‑TGTG‑CCGGGGTTTGTGGAAA‑TAAAAAAGACGA‑AGGCG‑ATTGAGATGAAAACCCG‑GG‑CGGAAAACGCCAGATAAAAAGTCTGAGGAGAAG‑TTTTCGTTGGTTTACAGACCCCCGCCTCTTCGGCCa                                                                                                                                                                                                            <  2:501864/239‑1 (MQ=255)                                    aacaacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGcc                                                                                                                                                                                                    >  1:537268/1‑250 (MQ=255)                                         caaTCAACACGGGGGAAGAAAAGAGATCGGCTCCGCAGAAAACACAAAGAATCAGCCCAAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAGGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATGGAGATGAAACCCCC‑GG‑CGGTAATCGCCAGATAGAAAGTTTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgca                                                                                                                                                                                              <  1:126418/251‑1 (MQ=255)                                          aaacAACCCGGGGGAAGAAAAGAGAGCGGCCCCGCAGCCAACAAAAAAACTCACCCCCAGGAAGAAGAAAAGGCCGCCAACCAGTTTTGCCAGGCAAGAGGG‑CTTTGGGCC‑TGTG‑CCAGTGGTTGGTGTTA‑TATGAAAGGCAA‑ATGCG‑ATTAAGTTAAAACCCCC‑CG‑GGATAATCTCAAGATAAAAAATGTAAGGGGGAT‑TTTTCGTTGTTTAAAAAACGCCCGCCTCTTCGGCCACTTTGGCCATTgcag                                                                                                                                                                                             <  1:158721/248‑1 (MQ=255)                                                         agaaaaGAGATCGGCCCCGCAGAAAACAAAAAGAACCAGCCCCAGAGAGAAGAAAAGGCCGCAAACCAGTTTTGCCATGCCGAGGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACAA‑ATGGG‑ATTGAGATAAAAACCCC‑GG‑GGGTAATCGCAAGATAGAAAGTCTTAAGCGGAT‑TTTTCGTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGAt                                                                                                                                                                              <  1:149715/248‑1 (MQ=255)                                                                   ccggcccccaggaaaacaCAAAGAACCAGCCCCAGAAAGAACAAAATCCCCCCACCCATTTTTGCCATGCCAAAGGG‑TATTGTCCC‑GGTG‑CCAGTGGTTGTTGTAA‑TAAAAAAGACAA‑ATGCG‑ATGGA‑AAGAAAACCCCCGG‑CGTTATCCGCCAAATAAAACGTCTTAAGCGGAG‑GTTTCGTTGTTTTATAAACACCCGCCGCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGtt                                                                                                                                                                    <  2:115748/237‑1 (MQ=255)                                                                         ccgaaaaaaacaAAAAGAACCACCCCCAGAGAAAAGCAAAGGCCCCCACCCGTTTTTGCCATGCC‑AAGGG‑GTATGTGCCCTGGC‑CCGGTGTTGGCTGAAA‑TATAGAAGAGAA‑ATGGG‑ATTGAGAAGAAAACCCG‑GG‑CGTTAATCGCCAGAAAAAATGTCT‑AAGGGGTTTGTTTCGTTTCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGc                                                                                                                                                              <  2:113342/243‑1 (MQ=255)                                                                          cgcaaacaccaaaaaGAACCACCCCCAGAGAGAACCAAATGCCGCCAACCATTTTTCCAATGCGGAAGGA‑CATTGTGCC‑TGTG‑CAAGTGGTTGTTGGAA‑TAAAGAAGACGA‑TTGAT‑ATTAAGA‑GAAAACCCCCGG‑GGTTAATCCCGAGATAGAACGTCTTCAGCGGAT‑TTTTTGTTGCTTTATAGACTCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGCTTGTCAGCt                                                                                                                                                             <  1:633486/239‑1 (MQ=255)                                                                             aGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGGAGTAAAAAATCAAAAGGGTTGTCAGCTTT‑c                                                                                                                                                         >  2:588832/1‑251 (MQ=255)                                                                                           aTCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAAAAGCCAGATAGAAAGTCTTAAGCGGAC‑GTTTCGTTGATTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGATTT‑CACACTAActctctc                                                                                                                                           >  2:474215/1‑251 (MQ=255)                                                                                                ccccAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTTGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCCGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CAAACTAACTATCTCTttat                                                                                                                                      >  1:316269/1‑251 (MQ=255)                                                                                                     gagagaAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGTT‑GTTTCGTTGTTTTATAGAAACCCGCCTCTTCCGCCACTTTGGCCATTGCAGCCGGAAGTAAAAGATAAAAAGGGTTGGTAGCCTT‑CACACTAACTCTCTCTTTATTAAGt                                                                                                                                 >  2:276101/1‑251 (MQ=255)                                                                                                     gagagaAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGt                                                                                                                                 >  1:379482/1‑251 (MQ=255)                                                                                                       gagaAGCAAATGCCGCCAACCAGTTTTGCCATGACGAAGGG‑CATTGTGCA‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTCATAGCCAGATAGAAAGGCTTAAGCGGAT‑GTTTCGTTTCTTTATAAACACCCGCCCCTTCGGCCACTTTGGCCATTTCAGAAGGAAGTAAAAAATCGAAAGGGTTGGTAGCTTTTCAC‑CTAACTCTCTCTTTATTaaagcg                                                                                                                               >  2:499770/1‑247 (MQ=255)                                                                                                       gaaaaGCAAATGCCGCCACCCATTTTTCCCTTGCCGAAGGG‑TATTGTGCC‑TGTG‑CCAGTGGTTGTTGTGA‑TATAGAAGACGA‑ATGCG‑TTTGAGATAAAACCCCC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGGGGAT‑TTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCg                                                                                                                               <  1:588832/248‑1 (MQ=255)                                                                                                       gaaaaGCAAAGGCCGCCAACCATTTTTGCCATGCGGAAGGG‑CATTGTGCC‑TGGG‑CCGGTGGGTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACCCC‑GG‑CGGTAATCGCCAGATATAAAGTCTTAAGCGGTT‑TTTTCTTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCg                                                                                                                               <  1:253619/248‑1 (MQ=255)                                                                                                            gCAAATCCCGCCAACCAGTTTTGCCAGGCCGAGGGG‑CTTTGTGCG‑TGTG‑CAAGTGGTTGGTGTGA‑TATAGAAGGCGA‑AGGCG‑ATTGAGATGAAAACCCG‑GG‑GGGGAACCGCCAGATAGAAAGTCTTAAGGGGAT‑TTTTCGTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                                          <  2:153571/251‑1 (MQ=255)                                                                                                                         aCCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTTCTTTATAGACACCCGCCTCTTCGGCCCCTTTGGCCATTGCAGCAGGGAGTAAAAGATCAAAAAGGGTGTCCGCCTT‑TACACTAACTCTCTCTTTATTTAGTCGGCGACGAGATACTAACa                                                                                                              >  2:506642/1‑250 (MQ=255)                                                                                                                              tttggcaatgcgaagggg‑CTTTGGGCC‑TGTG‑CCAGTGGTTG‑TGTGAATAAAGAAGAGGACGTG‑G‑TTTGAAAT‑AAAACCCC‑GGGCGGTAATCGCAAGAAAGAAAGTTTTAGGCGGAT‑GTTTCTTTGTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCa                                                                                                         <  2:306061/235‑1 (MQ=255)                                                                                                                                      ggCCGAAGGG‑TTTTGTGCC‑‑GTGGCCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGGA‑ATTGAGATAAAAACCCC‑GG‑CGGTAATCGCAAGAAAAAAAGTCTAAAGCGGAT‑GTTTCGTTGTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATg                                                                                                <  1:627639/250‑1 (MQ=255)                                                                                                                                       gCCGAAGGT‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTTATCTCCAGGTAGAAAATCTTAAGCGGGT‑GTTTCGTTGCTTTATAGACACCCGCCTTTTCTGCCACTTTGGCCATTTCAGCAGGAAGTAAAATATCCAAAGGGTTGTCCGCTTT‑CACCCTTACTCTCTCTTTATTACTTCGGAGACGACATACTAACAAAACATtatacatt                                                                                                >  2:497389/1‑246 (MQ=255)                                                                                                                                         tgcggagggtgttggggcttggcccAGTGGTTGCTGTTA‑TAAAGAAAACGA‑A‑GCGGATGGAGATGAAAACACC‑GG‑GTGTAATAGCCAGATAGAAAGCCTTAGGCGGAT‑TTTTCGCTGTTTTAAAAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCTGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGAt                                                                                                 <  2:627601/227‑1 (MQ=255)                                                                                                                                               g‑CATTGTGCC‑TGTG‑CCAGGGGTTGTTGTAA‑TATAGAAGAAGA‑ATGCG‑ATTGAGATGAAAACCCC‑GG‑CGGTAATCGCAAGATAAAAAGTCTTAAGCGGTT‑GTTTCGTTGCTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAtgga                                                                                       <  1:151447/251‑4 (MQ=255)                                                                                                                                                  aTTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ACGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAAATCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGat                                                                                      >  1:428252/1‑250 (MQ=255)                                                                                                                                                     gtgcc‑tgtg‑ccAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGCAAACACC‑GG‑CGGTTATCGCCAGATAGCAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGGCGGCGACGAGATACTAACAAAGCCTTATAGATGAGAAATTGatat                                                                                    >  2:113261/1‑249 (MQ=255)                                                                                                                                                      tgtcccggtgcagttGGTTCCTGTGA‑TAAAGAAGAAGA‑AGGGG‑TTTGTAATAAAAACCCC‑GG‑CGGAAACCGCAAGATAAAAAGTCTAAAGGGGAT‑TTTTCTTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                                                    <  2:324020/236‑1 (MQ=255)                                                                                                                                                                                        gccGA‑ATGCG‑ATTGAGATGAAACCCCC‑GG‑CGGTAATCGCCAGATAAAAGGTCTTAAGGGGTT‑GTTTGTTTGCTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAAc                                                  <  1:506642/249‑1 (MQ=255)                                                                                                                                                                                         aCGA‑ATGCG‑TTTGAGAAAAAACCACC‑GG‑CGTTAATCGCCAGATAGAAAGTCTTAAGCGGTT‑GTTTCTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCAACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACt                                                 <  1:4051/251‑1 (MQ=255)                                                                                                                                                                                            a‑aTGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGc                                              >  1:149553/1‑251 (MQ=255)                                                                                                                                                                                              aTGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGCACGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAACATCAAACGGGGTGTCAGGTTT‑AAAACTCACTCTCTCTTTTTTAAGTCGGCGACGAGATAATAACAAAGAATTATAGATGAGAAATTGATATAGATCCAATCCCGCCTGGCTTATAGGCCCGTAACTctac                                             >  1:242671/1‑248 (MQ=255)                                                                                                                                                                                                  g‑ATTGAGATGAAACCCCC‑GG‑GGGTAATCGCAAAAAAAAAATTCTAAAGGGGTT‑GTTTCGTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGt                                         <  1:474821/251‑1 (MQ=255)                                                                                                                                                                                                          aTGAAAACACC‑GG‑CGGTAATCGCCCGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGGTTTTAtt                                  >  2:56416/1‑250 (MQ=255)                                                                                                                                                                                                                  ccc‑GG‑CGGTAACCGCAAGATAGAAAGTCTAAGGCGGTT‑GTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg                          <  1:56416/250‑1 (MQ=255)                                                                                                                                                                                                                  ccacggggggTAATAGCAAGAAAGAAAGTATGAAGGGGAT‑GTTTGGTTGTTTAATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa                            <  2:428252/243‑1 (MQ=255)                                                                                                                                                                                                                   cc‑GG‑GGGTAATCGCCAGATAGAAAGTTTTAAGGGGAT‑GTTGCGTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACgg                         <  1:631932/251‑1 (MQ=255)                                                                                                                                                                                                                                cGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTAGGCCCCTTTGGCCATTGCAGCAGGAAGTAAAAGATCACAAGAGTTGTCTGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACTGGTTAAATAttt              >  2:1320/1‑251 (MQ=255)                                                                                                                                                                                                                                             gtttttaaggggttttttCGTTGTTTTATAACCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCaaat  <  2:149553/236‑1 (MQ=255)                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     AACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGG‑CATTGTGCC‑TGTG‑CCAGTGGTTGCTGTGA‑TATAGAAGACGA‑ATGCG‑ATTGAGATGAAAACACC‑GG‑CGGTAATCGCCAGATAGAAAGTCTTAAGCGGAT‑GTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTT‑CACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAAT  >  NC_000913/954138‑954615

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