Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,745,089	0	A	C	34.0%	70.0 / 11.9	53	V84G (GTG→GGG)	trmD	tRNA m(1)G37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/9);  new base C (0/17);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 3.52e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  NC_000913/2744843‑2745335                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         gCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑tt                                                                                                                                                                                                                                                     <  2:480232/251‑1 (MQ=255)   cGAGGCTTCATGTACCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACAATTCTTCGTCCATTTCGGTTTGGATCACGAGCTCATCTATACCTTCGTGGCGACCGCACACCAGAATCACTTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTAGTG‑ACAGATAAATCACA‑TTTg                                                                                                                                                                                                                                                   >  1:350811/1‑251 (MQ=255)    gAGGCTTCATGTCCCAGTCCTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTgc                                                                                                                                                                                                                                                  <  2:29590/251‑1 (MQ=255)              gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCAccc                                                                                                                                                                                                                                       >  1:11488/1‑251 (MQ=255)              gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCAccc                                                                                                                                                                                                                                       <  2:326286/251‑1 (MQ=255)              gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCAcc                                                                                                                                                                                                                                        <  1:148988/250‑1 (MQ=255)                           cGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCCCCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTtg                                                                                                                                                                                                                          <  1:63338/251‑1 (MQ=255)                                             ggAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGCTCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTccc                                                                                                                                                                                                        >  2:411566/1‑251 (MQ=255)                                              gAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTgcg                                                                                                                                                                                                                                                 >  2:384642/1‑210 (MQ=255)                                              gAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTgcg                                                                                                                                                                                                                                                 <  1:384642/210‑1 (MQ=255)                                              gAGTCAATCAGCGTCAATGCTGGTAAATAACCAACACTGAGAACGAAATCGACGATTGACCATTCTTCGTCAATTTCGGGTTGGATCAAGAGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCCATTTTTGATTCGTTG‑CCAGCTCGCTGAA‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑AGAGATAAGTCACA‑TTTgag                                                                                                                                                                                                                                                 >  2:384635/1‑208 (MQ=255)                                                   acccaggcgctttGCTGGTAACCCCCCCCCCCTGAGAACTTATCCGCCTATTGACTGTTTTTCGTTACTTTCGTTTTGGATCACGCGCTCATCTATCCCTTCTTAGCCCCCGCCCCCCAGAATCAATTTTTAATCCGTTC‑CAAGTTCGCGGAC‑GCCCGCTGGATCAAGCTTGCGTCCGT‑GTGGGG‑ACAGATAAATCCCC‑TTTCCCCC‑TCCACCCGCCGCTGCTTTTGCTGCATGAATGGCGTCCCGCAAgg                                                                                                                                                                                                  <  2:331668/240‑1 (MQ=255)                                                      cAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGt                                                                                                                                                                                                <  1:13776/250‑1 (MQ=255)                                                       aGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                                              >  1:296327/1‑251 (MQ=255)                                                          gTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                                                                                                                                                                                           <  1:395675/251‑1 (MQ=255)                                                                 cgggTACCTCACCACCACTGAGAACTAAACCGCCGATTGACCATTCTGCGGCAATTTCGGTTTGGACAACGCGCTCATCTATACCTTCGTAGCGACCGCACACAAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GGGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCAtt                                                                                                                                                                                    <  2:223558/249‑1 (MQ=255)                                                                      aaccccccaccccTGAGAACGTAATCGCCGATTGCCCATTCTGCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATGCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCG‑GTGTTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACa                                                                                                                                                                                <  1:451419/244‑1 (MQ=255)                                                                          ccccaccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GGGGTG‑ACAGATAAATCTCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc                                                                                                                                                                           <  1:489903/249‑1 (MQ=255)                                                                              cccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTc                                                                                                                                                                       <  1:168239/250‑1 (MQ=255)                                                                                 acTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgc                                                                                                                                                                    >  1:201913/1‑251 (MQ=255)                                                                                  ccgggagtgtatcccccGATTGACAATTGTTCGCCAATTTCGGTTTGGATCCCTCCCTCATCCATCCCTCCGTACCCCCCCCACACAAGAAACAATTTTTGATTCGTTG‑CCAGTTCCCTGCC‑GCCCGCTGGATCAGGCTGTCGTCCCT‑GTGGTG‑ACAAATAAATCCCC‑TTTGCGCC‑TTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgcc                                                                                                                                                                   <  2:398188/236‑1 (MQ=255)                                                                                               cgccatttaaCCATTCTTCTTAAATTCCGGTTTGGACCACGCCCTAACCTATCCCTTCTTAGCGACCGCACACCAGAACCAATTTTTGATCCGTTG‑CCAGTTCGCTGCC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑TTGGTG‑ACAGAAAAATCCCC‑TTTGCCCC‑TTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGCCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACg                                                                                                                                                       <  1:59019/241‑1 (MQ=255)                                                                                                gctgagtgcctaTTCTTCGTCATTTGCGGTTTGGATCACGCGCTCATTTATTCCTTCGTAGCGACCGCACCCCAAACTCAATTTTGGTTTCTTTG‑CCAGTTCGCTGAC‑GCCCGCTTGAACAGGTTTGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTGGCGCC‑TCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGAt                                                                                                                                                     <  1:267273/240‑1 (MQ=255)                                                                                                          cATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTCTACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGAATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCTCCGGTCCTCCGACGTAAGGCCGATCGTCCacggt                                                                                                                                           >  2:359681/1‑251 (MQ=255)                                                                                                           aTTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCCATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                                           >  1:408382/1‑250 (MQ=255)                                                                                                                  gTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                                          >  2:245426/1‑196 (MQ=255)                                                                                                                  gTCAATTGCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                                          <  1:245426/196‑1 (MQ=255)                                                                                                                    cAATTTCGGTTTGGGTCACGCGCTCATCTATACCTCCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CACGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCACGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCCCGGTCCGGAGCCgg                                                                                                                                 >  2:211732/1‑251 (MQ=255)                                                                                                                        ttCGGTTGGGACCCCGCGCGCATCTATCCCTTCGAAGCGCCCGCCCACAAAAACCAATTTTTGATTCGTTC‑CCGGTTGGCTGAC‑GCCCGCGTGATCAAGCTTGGGGCCGT‑G‑GGTGGACA‑AAAAATCCCCCTTTGCCCC‑TTCACCCGCCGCGGCTTTTGCGGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCa                                                                                                                              <  1:210059/250‑1 (MQ=255)                                                                                                                          cGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTCCGGTGCCGGTCATGc                                                                                                                           >  2:91554/1‑251 (MQ=255)                                                                                                                               tGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGa                                                                                                                       >  1:106024/1‑250 (MQ=255)                                                                                                                                 gttcccgcgCTCATCTGTTCCTTCGTACGGACCGCACGCAAAAATCAATTTTTGTTTCGTTG‑CCAGCTCG‑TGACCGCCCTCTTGACCAAGCTTGC‑TCCCTTGGGGGG‑ACAAATAAACCTCT‑GTTGCCCC‑TGCCCCCGCCGCGGCTTTTGCTGCATGAATGGCTGCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGTTCATGCGTGAAGt                                                                                                                    <  1:489365/246‑1 (MQ=255)                                                                                                                                  accccgcgTTCATCTATCCCTTCGTCGCCCCCGCACAACAGAATCAATTTGTTTTGCGTCG‑CCTGTGCGCGCAC‑GCCCTTTTGACCAAGCTGGCGCCTCT‑GTGGTG‑ACAAATAATTCCCC‑TGCGCCCC‑TCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTc                                                                                                                   <  2:11488/247‑1 (MQ=255)                                                                                                                                    cACGCGCTCATCTATACCTTCGTAGCGACAGCACACCAGCATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAACTCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGCCGATCGTCCACGGTACGGTGCAGGTCATGCGTGAAGTcgc                                                                                                                 >  1:175052/1‑251 (MQ=255)                                                                                                                                            atccTATACCTTCGTAGCGCCCGCACACAAGAATCAATTTTTGATT‑TTTGTCCATTTCGATGAC‑CCCCGCTTGAACAGGTTTGCGTCCTT‑GGGGGA‑ACAGATAAATCAC‑‑TTTGCCCCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGGAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACt                                                                                                          <  1:33014/247‑1 (MQ=255)                                                                                                                                              tCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTAGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGCAGTCGCGAGGACTcc                                                                                                        >  1:47894/1‑250 (MQ=255)                                                                                                                                                           accGACCGCACACCAGAATAAATTTTTGATTCGTGG‑CCAGTTCGCTGAC‑CCCCGCTTGATCAAGCTTGCGCCCCT‑GTGGTG‑ACAAAAAAATCTCC‑TGGGCTCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCt                                                                                           <  1:283362/248‑1 (MQ=255)                                                                                                                                                           aGCGACCGCACACGAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑CCCCGCTTGATCAAGTTGGCGTCCCT‑GTGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCt                                                                                           <  1:251838/250‑1 (MQ=255)                                                                                                                                                                ccgccccccAGAATTAATTTTGGATTGGTTG‑CCAGTCCGCTGCC‑GCCCGCTTGATCAAGCTGGCGTCCCT‑GGGGTG‑ACAGAAAAACCCCC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcagcag                                                                                     <  1:389922/244‑1 (MQ=255)                                                                                                                                                                  gCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGCC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCCCGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTAAGAAGGc                                                                                   >  2:24832/1‑251 (MQ=255)                                                                                                                                                                   cacaCCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGcc                                                                                  >  2:377707/1‑251 (MQ=255)                                                                                                                                                                      cccAGAATAATTTTTTAATGCTTTG‑CCAGTTCGGTGCC‑CCCCGCTGAATCAGCCTTGCGCCCCT‑GGGGTA‑ACAGAAAAACCCCC‑TTCGCGCT‑TTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAtt                                                                               <  2:106024/250‑1 (MQ=255)                                                                                                                                                                                       atgtgTTG‑CAAGTTCGCTGAC‑GCCCGTTTGAACAAGCTTGCGCCCCT‑GTGGAA‑ACAGAAAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGGCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAg                                                              <  1:193540/247‑1 (MQ=255)                                                                                                                                                                                       aTTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                >  1:162405/1‑249 (MQ=255)                                                                                                                                                                                                               cccGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAACTCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTCATCGGAAATTgc                                        >  2:386390/1‑251 (MQ=255)                                                                                                                                                                                                                    ttaatcaggctggcttCCCT‑GGGGTG‑ACAGATAAATCCCC‑TTTGCGCC‑TTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCgg                                    <  1:65503/235‑1 (MQ=255)                                                                                                                                                                                                                        ccaggctggCGTCCTT‑TTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCCCCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACAt                               <  1:126917/243‑1 (MQ=255)                                                                                                                                                                                                                             cTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTc                            >  1:370364/1‑249 (MQ=255)                                                                                                                                                                                                                                  gTCCCT‑GTGGTG‑GCAGATACATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGACGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                     >  2:489496/1‑251 (MQ=255)                                                                                                                                                                                                                                            gtg‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAAtt            >  1:153394/1‑251 (MQ=255)                                                                                                                                                                                                                                              g‑ACAGAAAATTCCCC‑TTTGCCCC‑TCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTAt          <  1:31024/251‑1 (MQ=255)                                                                                                                                                                                                                                                aCAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATg         >  1:160760/1‑251 (MQ=255)                                                                                                                                                                                                                                                       aaTCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATcc  >  1:445668/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         GCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTG‑CCAGTTCGCTGAC‑GCCCGCTTGATCAAGCTTGCGTCCCT‑GTGGTG‑ACAGATAAATCACC‑TTTGCGCC‑TTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCC  >  NC_000913/2744843‑2745335

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