Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,410,037	T→C	7.3%	intergenic (+31/‑96)	rlmB → / → yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,037	0	T	C	7.3%	158.5 / 3.0	55	intergenic (+31/‑96)	rlmB/yjfI	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (31/20);  new base C (2/2);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGAGTAGTATCAGCAATCCAGATGA  >  NC_000913/4409738‑4410332                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                          cTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGCAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTTGACGCGTTTCTTCCATGAACGTTTCGGTTGCGCCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAACTAATTTATCAGAATTATGAAACGCCATCCAGAgtt                                                                                                                                                                                                                                                                                                          >  1:217803/1‑301 (MQ=255)       tttggcgcACCATGCGTTTGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATgg                                                                                                                                                                                                                                                                                                    <  1:147064/297‑1 (MQ=255)        tGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑t                                                                                                                                                                                                                                                                                                 >  2:519157/1‑301 (MQ=255)                          gTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCt                                                                                                                                                                                                                                                                               >  2:191401/1‑301 (MQ=255)                             gCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGCTCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGCCCGGACTTTGCTTATTTTAAGCGGTGCGCCAGCGCAGCTAATTTATCAGAATTCTGAAAAGCCATCCAGATTTGGCTGG‑‑TTTTTTTTTGTCTATCGCTgtt                                                                                                                                                                                                                                                                            >  1:602491/1‑299 (MQ=255)                             gCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTTCGCCAGCGCAGCTAATTTCTCAGCATTATGAAAAGCCATCCAGCTTTGGCAGG‑‑TTTTTTTTTTTCTAAAGCTGGt                                                                                                                                                                                                                                                                            >  1:94704/1‑301 (MQ=255)                                ggAAGGGAATATCGGGGTCGTCGGTACGGCAGGCGAGGCGGTTCATACACTCTATCAGAGCAAAGTGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                                                                                                                                                                         <  1:359929/301‑1 (MQ=255)                                 gAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTNNNNAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTNNCGGTTGCGACCGGAANNNNCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAAAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCtaacactggaaagt                                                                                                                                                                                                                                                                        >  1:321132/1‑289 (MQ=255)                                  aaGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                                                                                                                                                                                      <  2:170917/301‑1 (MQ=255)                                   agagAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCTCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                                                                                                                                                                                      <  1:266082/301‑1 (MQ=255)                                     tcagacgtgtgctcttccgatctgCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                                                                                                                                                                                                    <  2:166485/278‑1 (MQ=255)                                      gattatTTGGATGGTCGGTACGGCAGGCGAGGGGGATCATACACTCGATCAGGGCAAAATGACCGGACGTCTGGCGCTGGGGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                                                                                                                                                                                                                   <  1:238419/298‑1 (MQ=255)                                             tGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGACCATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCTGACTTTGCTTCTTTGAAGCGGTGCGCCACCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGG‑‑TTTTTTTTTGTCTAAAGCTGGTAAGATAATTACATAtt                                                                                                                                                                                                                                                            >  1:522646/1‑301 (MQ=255)                                                  acgacgctgttccgttctggcgGTTCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                                                                                                                                                                                                                       <  1:504930/283‑1 (MQ=255)                                                     cgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAataa                                                                                                                                                                                                                                                    >  2:228531/1‑300 (MQ=255)                                                       gTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGATAAGATAATTACGTATTGAAAATAttt                                                                                                                                                                                                                                                  >  1:135846/1‑300 (MQ=255)                                                            gCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAtagaacggaagagcgacgagcagg                                                                                                                                                                                                                                             >  1:166288/1‑278 (MQ=255)                                                                ccgatctggATCATACACTCTTTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAGCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                                                                                                                                                                                                         <  1:276012/294‑1 (MQ=255)                                                                    ggcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaaagatcggaagagcacaca                                                                                                                                                                                                                                     >  2:505602/1‑283 (MQ=255)                                                                       ggATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAGCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAagatcgg                                                                                                                                                                                                                                  >  2:276315/1‑298 (MQ=255)                                                                        taTCATACACTCTATCAGAGCAACATGACCGGACGTCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                                                                                                                                                                                                                                 <  2:8433/300‑1 (MQ=255)                                                                                     aTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAACAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                                                                                                                                                                                                                    >  2:378111/1‑301 (MQ=255)                                                                                       caggggaaagtggcgggtcGCCGGGGGCTGGTGATGGGTGCGGAAGTGGAAGGTATGCTTCGCCTGACTCGGGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                                                                                                                                                                                  <  1:518485/283‑1 (MQ=255)                                                                                               aaTGACCGGACGCCTGGCGTTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGACTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                                                                                                                                                                          <  2:135971/301‑1 (MQ=255)                                                                                                            ttgtcgttGGTTATGGGTGCGGAAGGTGAAGGTTTGCGTCGCCTGACTCGTGAACATTGCGATGAGTTTATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                                                                                                                                                                             <  2:218046/294‑1 (MQ=255)                                                                                                            cTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGAAATAATATATGAAATGGaaat                                                                                                                                                                                             >  1:460318/1‑299 (MQ=255)                                                                                                                  cTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATAGAATCAGTTga                                                                                                                                                                                       >  2:16458/1‑300 (MQ=255)                                                                                                                                                cGTTGCCTGACTTGTGGTCATTGTGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTg                                                                                                                                                         <  1:228288/301‑1 (MQ=255)                                                                                                                                                            cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCAAGCGACGGAGCTGAAAACTGTACCTGAAAAAAAAAttt                                                                                                                                             >  1:248185/1‑300 (MQ=255)                                                                                                                                                               ttaCATTGTGATTAGTTTATCAGCGTCCTGGTGGTTGGAAGCGTTTCTTCGTTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGGGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATg                                                                                                                                          <  2:321536/299‑1 (MQ=255)                                                                                                                                                               gAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATg                                                                                                                                          <  2:73125/301‑1 (MQ=255)                                                                                                                                                                   tttGGGATGGGTTGATAAGCATCCTGATGGGTGGAAGCGTTTTTTTACTGAACGTGTCGGTTGCGACTGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAAc                                                                                                                                      <  2:94828/300‑1 (MQ=255)                                                                                                                                                                         tttaGTTGTTCTGCTTTCCTATTGGTGGAAGGGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAg                                                                                                                                <  2:460882/298‑1 (MQ=255)                                                                                                                                                                               tGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGaa                                                                                                                          >  2:537483/1‑301 (MQ=255)                                                                                                                                                                               tGATCAGCATCCCGATGGCTGGAAGCGTTTATTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTCAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGAAACAAATaacaaaaa                                                                                                                          >  2:255468/1‑294 (MQ=255)                                                                                                                                                                                         cccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCattaata                                                                                                                >  2:643393/1‑300 (MQ=255)                                                                                                                                                                                                gCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAAGTTCTCAGAATTATGAGAAGCCATCCAGATTCGGATGG‑‑TTTTTTTTNNNNNATAGCAGGTAAGATNNNNACGTATNNNNNNTATTCATTTNNNNAATATTTAAANNNNNNCTATATGGAATAANNNNTAACATGGAANCCGTTGGCGCTNNNNNNNGCGATGCAAACAGTACCTGAAAACAAACTTAATGTAACAAAAAACAAGAACACAATAAATATTaaaa                                                                                                         >  1:250987/1‑301 (MQ=255)                                                                                                                                                                                                                cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAACAATTGATGTAACAAATCGCAAAAACGCAATAATTATTAAAATAAATAATAATGGCGa                                                                                         >  1:61344/1‑301 (MQ=255)                                                                                                                                                                                                                            tggtttctattgtaaTTTGCTTTTTTTAAGCGGTGCGGCAGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCaa                                                                             <  2:248452/288‑1 (MQ=255)                                                                                                                                                                                                                                  cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTGACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCt                                                                       >  2:533828/1‑301 (MQ=255)                                                                                                                                                                                                                                                    ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGACCGCATTAATTATTCAAATGAATGATTATGGAGATTTGCAAATAAATATTCTTTTTACTTCCCGCCAAAt                                                     >  1:408621/1‑301 (MQ=255)                                                                                                                                                                                                                                                     tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAAta                                                    >  2:251702/1‑300 (MQ=255)                                                                                                                                                                                                                                                           ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATc                                              >  2:41943/1‑301 (MQ=255)                                                                                                                                                                                                                                                              gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCTATTTGCAAATAAATATTATTTTTACTTCCCGCCCAATGATTATCGaa                                           >  1:376214/1‑301 (MQ=255)                                                                                                                                                                                                                                                               ttttttatttttttttttcctatcacaaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGa                                            <  2:552179/271‑1 (MQ=255)                                                                                                                                                                                                                                                               cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTACCAAATAACGAGAACGCATTAATTATTAAAATGAATGATTAAGGCGATTTGCAAATAAAAATTCTTTTTACTTCCCGCAAAATAAATATCaaaa                                          >  1:501310/1‑301 (MQ=255)                                                                                                                                                                                                                                                                 tcagacgtgtgctcttccgatctcaaaaaaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGTTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAAcc                                        <  2:23649/273‑1 (MQ=255)                                                                                                                                                                                                                                                                         gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCTAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTAtttt                                >  2:553639/1‑300 (MQ=255)                                                                                                                                                                                                                                                                               tttttaGAATTATTAAAAGCCATCCAGATTTTGATTG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTGACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGt                          <  1:533155/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                      aaTTATGAAAAGTCTTTCTTATTTGGATGG‑‑TTTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTTATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGagtagt                   <  1:251435/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                        caaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTAACaaaaatcggaaaagagtcgggtagggg                   >  1:551506/6‑278 (MQ=255)                                                                                                                                                                                                                                                                                        caaaaaaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGTTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAAccagatcggaagagcgtcgagtagg                 >  1:23598/6‑278 (MQ=255)                                                                                                                                                                                                                                                                                            gAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTTCGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTTTACCTGAACAAAATATTGCTGTAACAAATACCACGACCGCATTAATTCTTAAACTGAATGATTATGGCGATTTGCAAAACACTCTTCTTTTTACTTCCCGCCAAATAATTCACGAAACCTTTCTTTCTCCGGTGcgtagcacactc             >  1:456662/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                               aaGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCCAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATATGTCCGGTGATTAGTATCAGCaaa          >  1:461233/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                    aTCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTAATTGTCCGGTGAGTAGTAACAGAAATCCAGa     >  1:543758/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                       cAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGAGTAGTATCAGCAATCCAGATGa  >  2:609661/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                          CTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGAGTAGTATCAGCAATCCAGATGA  >  NC_000913/4409738‑4410332

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