Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,641,131	T→C	L64L (TTG→CTG)	yjtD →	putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,641,131	0	T	C	97.7%	138.4 / ‑4.5	44	L64L (TTG→CTG)	yjtD	putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/18);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

ATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTG  >  NC_000913/4640901‑4641373                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ctCACCCTGTTCGCGGCTATTTGGTCAAGCAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGCGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGa                                                                                                                                                                                                                                 >  1:380128/2‑250 (MQ=255) accACCCGGTGGGCGGTTATTGGGTCAAGAAGAGGGGAAATATGCGTATAACGATTATTCGGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGAt                                                                                                                                                                                                                                <  2:5090/249‑1 (MQ=255)                       ggtCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGccacc                                                                                                                                                                                                           >  1:209184/3‑250 (MQ=255)                         tCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACt                                                                                                                                                                                                        >  2:135665/1‑251 (MQ=255)                               agagGGGCAATATGCGTAGAACAATTATTCGGGTCGCACCCGCCAGAGCAGAAAATATGGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTgcgcgc                                                                                                                                                                                                  <  1:502458/251‑1 (MQ=255)                                       aaTATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGc                                                                                                                                                                                          >  2:65662/1‑251 (MQ=255)                                                            cTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAACACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTacgctac                                                                                                                                                                     >  1:224/1‑248 (MQ=255)                                                                                       atTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTag                                                                                                                                           <  1:51037/250‑1 (MQ=255)                                                                                         tGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTActacgctacg                                                                                                                                                                    <  1:143484/223‑1 (MQ=255)                                                                                         tGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTAATacgctacg                                                                                                                                                                    >  2:143484/1‑223 (MQ=255)                                                                                          ggggCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCTATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACAGTGATCTGGTGATATTATTGATAA‑TATTAAGGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCATTTGAACTGGTGCCGCTGTTAGAGGa                                                                                                                                       <  2:202804/251‑1 (MQ=255)                                                                                            ggCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAggca                                                                                                                                      >  1:334224/1‑248 (MQ=255)                                                                                                 cggcgCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGCCCCTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTAAGCCAGATGAACTGGTACCGCTGTTAGAGGAAAAATCt                                                                                                                                >  2:1428/1‑251 (MQ=255)                                                                                                 cggcgCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTCCGCCAGTAGAACTGGta                                                                                                                                                     >  2:465565/1‑229 (MQ=255)                                                                                                 cggcgCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGta                                                                                                                                                     <  1:465565/230‑2 (MQ=255)                                                                                                      cgGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATg                                                                                                                           >  2:257135/1‑251 (MQ=255)                                                                                                       gggCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGCTAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTATAGGAAAAATCTTCATgg                                                                                                                          >  2:30366/1‑251 (MQ=255)                                                                                                                 aaCGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATg                                                                                                                <  1:173044/251‑1 (MQ=255)                                                                                                                    gATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCcg                                                                                                             >  1:433125/1‑251 (MQ=255)                                                                                                                             gaGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCCCCCGCTGGGTCGCAAATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGgtgt                                                                                                    <  2:141864/250‑1 (MQ=255)                                                                                                                                      gaggtttgacgctAGTCAGGCACACCTGGATCCAGCCACCCGCTGGGCCGCACATGGGTCTGGTGATATTATTGATAA‑TATTAAAGTTTCCCCGACACTGGCTGAATCTTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAAGCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTcgcg                                                                                           <  2:355040/239‑1 (MQ=255)                                                                                                                                      gCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTcg                                                                                             <  1:321744/249‑1 (MQ=255)                                                                                                                                            tGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGAt                                                                                      <  1:207556/250‑1 (MQ=255)                                                                                                                                               cGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCAGCCACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGACACCACTgc                                                                                                                                                                                                      >  1:451413/1‑135 (MQ=255)                                                                                                                                                            caacTGGAGCCAGCCAACCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGTAATATCATTACTACGCTACGCCAGTAGAACTGGTAAAGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACg                                                                      >  1:179066/4‑250 (MQ=255)                                                                                                                                                                           aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGc                                                      >  1:232419/1‑251 (MQ=255)                                                                                                                                                                           aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGAGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAACATCTTCATGGATGAGACATGCCGCGCTGGTGTTTGGTCGCGAAGCTTCCGAGCTGACTAACGAAGAGTTAGCGTTGGc                                                      >  2:333682/1‑251 (MQ=255)                                                                                                                                                                               gtgggggcGCACAGGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACCCTGGCTGAATCGTTACACGATGTCGATTTCCCTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGAc                                                  <  2:491701/244‑1 (MQ=255)                                                                                                                                                                                        cacaTGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACt                                         <  1:3880/251‑1 (MQ=255)                                                                                                                                                                                        cacaTGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACt                                         >  1:337024/1‑251 (MQ=255)                                                                                                                                                                                             ggATCGGGTGATATTAGTGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGGCGCCACCACGGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGGCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGgtgt                                    <  2:70940/251‑1 (MQ=255)                                                                                                                                                                                                   gatgaatttattGATAACTAT‑AACGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTGCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACTCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAt                              <  2:94522/244‑1 (MQ=255)                                                                                                                                                                                                           tattGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGca                       <  1:135665/250‑2 (MQ=255)                                                                                                                                                                                                           tattGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGGTTGGTCGCGAAGATTCCGGGCTGACTCACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAtggttgacg                      >  2:123695/1‑246 (MQ=255)                                                                                                                                                                                                            attGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc                        >  1:108888/1‑248 (MQ=255)                                                                                                                                                                                                             ttGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc                        >  1:94365/1‑247 (MQ=255)                                                                                                                                                                                                              tGATAA‑TATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGAt                    <  1:257135/250‑1 (MQ=255)                                                                                                                                                                                                                          aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTCCTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCAAGGATGAGCCATGCCGCGCTGGTGTTTAGACGCGAAGATTACGGGCAGACTAACGAAGAGTTAGCGTTGGCTGACGATCTTACTGGTGTGCCACTGGGGGAAGATTATCCTTCACTc        >  2:329035/1‑251 (MQ=255)                                                                                                                                                                                                                          aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGa                                                                            >  1:447721/1‑183 (MQ=255)                                                                                                                                                                                                                          aaGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTc        >  1:486673/1‑251 (MQ=255)                                                                                                                                                                                                                          aaGTGTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGa                                                                            <  2:447721/183‑1 (MQ=255)                                                                                                                                                                                                                            gTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTACGGGCTGAATAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCa       >  2:113151/1‑250 (MQ=255)                                                                                                                                                                                                                             ttttCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGCCTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGCTGGTGGCAGATTATCCTTCGCTCAAt     >  2:214402/1‑251 (MQ=255)                                                                                                                                                                                                                                ccccGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTg  <  2:87157/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAA‑TATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTG  >  NC_000913/4640901‑4641373

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