Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,640,924	A→G	intergenic (+153/‑18)	ytjE → / → yjtD	protein YtjE/putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,640,924	0	A	G	97.5%	126.1 / NA	41	intergenic (+153/‑18)	ytjE/yjtD	protein YtjE/putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/14);  minor base T (0/1);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTA‑TTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGT  >  NC_000913/4640679‑4641165                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cTAGTTGCATCATACAACTAATACACGGGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTCAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCagg                                                                                                                                                                                                                                               >  2:204285/1‑249 (MQ=255)         atcatACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTgg                                                                                                                                                                                                                                       <  2:6144/251‑1 (MQ=255)                    aataaACGAGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACACTTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAACTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTCCATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTAGCGGTTA‑TTTGGTCAAGAAGAGTGGCACTATGCGt                                                                                                                                                                                                                             >  1:403553/1‑250 (MQ=255)                                                 tttatttTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGcc                                                                                                                                                                                                  >  1:337615/1‑248 (MQ=255)                                                                        aaGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcggc                                                                                                                                                                        >  2:390702/1‑251 (MQ=255)                                                                                  aaGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCGTTTTTAATCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGTTTA‑TTTGGTCAATAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATg                                                                                                                                                              <  2:165539/251‑1 (MQ=255)                                                                                                 tttttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTCCGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATCTTGGGGCCGCGGCGCGGGCAATGAAAACGATggtattt                                                                                                                                               >  1:390440/1‑246 (MQ=255)                                                                                                    ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGCAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGc                                                                                                                                            >  1:530579/1‑251 (MQ=255)                                                                                                    ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGTGGCAGCGGCTCGGGCaa                                                                                                                                                                <  1:180207/231‑1 (MQ=255)                                                                                                    ttAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCaa                                                                                                                                                                >  2:180207/1‑231 (MQ=255)                                                                                                     tAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGc                                                                                                                                            <  1:303010/250‑1 (MQ=255)                                                                                                         acttCAACCGTTAGTACCATCAGGAAATAGTTTCGCCTTTTTTACGCTAAGTACAGGGATTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACAGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCAGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGAGGGGTTTTAGCGaact                                                                                                                                       >  2:437795/4‑248 (MQ=255)                                                                                                           gtCAACCGTTAGTACATTCAGTAAAGAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGctaactg                                                                                                                                      <  2:28716/249‑7 (MQ=255)                                                                                                              aaCCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCgg                                                                                                                                   >  1:16731/1‑250 (MQ=255)                                                                                                                       tACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGAt                                                                                                                          >  1:50148/1‑250 (MQ=255)                                                                                                                                              ttttttAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGagccagcc                                                                                                  >  1:349534/1‑251 (MQ=255)                                                                                                                                               tttttAAGCAAAGTAAAGGGATTTTTCTGCGACTTACTTTAAGGATTTGTAAATTCGCGCCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑GTTGGTCAAGAAGAGTGGCAATATGCTTATAACGATTATTCTGTTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGATAACGATGGGTTTTAGCGAACTGCGGATTGTCGAGTGGCAGGCACACCTGGaaccagc                                                                                                   <  2:563152/249‑1 (MQ=255)                                                                                                                                                   tAAGCTAAGTAAAGGGCGTTTTCTGCGCCTTACGTTAAGGATTTGTAAATTCGAACCGCGTAATAAGTTGACAGTGACCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAAGATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCAccc                                                                                              <  2:143946/250‑1 (MQ=255)                                                                                                                                                     aGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAACTATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGAGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCt                                                                                           >  2:3815/1‑251 (MQ=255)                                                                                                                                                       cTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTgg                                                                                         >  2:141697/1‑251 (MQ=255)                                                                                                                                                         aaGTAAAGGGATTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCATAGCAGAAACTATTGGGGCAGCGGCGCGGGCAATGAACACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCAGGAGCCAGCCACCCGCTgggg                                                                                       >  1:110734/1‑250 (MQ=255)                                                                                                                                                           gTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                                                                                    <  1:375327/124‑1 (MQ=255)                                                                                                                                                           gTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                                                                                    >  2:375327/1‑124 (MQ=255)                                                                                                                                                              aaGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGcac                                                                                  >  1:403144/1‑251 (MQ=255)                                                                                                                                                                 ggCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATg                                                                               <  1:135122/251‑1 (MQ=255)                                                                                                                                                                 ggCTTTTTATGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATg                                                                               <  2:213796/251‑1 (MQ=255)                                                                                                                                                                      tttCGGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGGTCGCGGTTA‑TTGGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACTCCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATc                                                                           <  2:530579/250‑1 (MQ=255)                                                                                                                                                                       ttCTGCGACTTACGTTAAGGATTTGTAAGTTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTg                                                                         >  2:162424/1‑251 (MQ=255)                                                                                                                                                                        actggtaCTTACGTTAAGTATTTGTAAATTCGCACCGTGTAATAAGTTGACATTGATCCCCCTTTTCGCGGTTA‑TTTGTTCAAGAAGAGTGGTAATATGGTTAGAACGATTATTCGGGGCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAAGGAAAACGATGGGGTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGTAAGCCACCCGCTGGGTCGcacaatgaactg                                                                         <  2:262232/244‑9 (MQ=255)                                                                                                                                                                            cGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGCGCCAGCCACCCGCTGGCTCGCACATGGATCTGGTgct                                                                    >  2:148228/1‑249 (MQ=255)                                                                                                                                                                                     ttAAGGATTTGTAAATTAGCACCGCGTAATACGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGATGGGTCGCACATGGATCTGGTGATATTATTGat                                                           >  1:458450/1‑251 (MQ=255)                                                                                                                                                                                             gtGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGTTTCGCGGTTA‑TTTTGTCAAGAAGAGTGGAAATATGCGTATAACGATTATTCTGGGCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACACGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTa                                                    <  2:148156/249‑1 (MQ=255)                                                                                                                                                                                               gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTaa                                                   >  1:537279/1‑249 (MQ=255)                                                                                                                                                                                                  aaTTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGttt                                              >  2:130806/1‑251 (MQ=255)                                                                                                                                                                                                           cggttaataaGTGGACAGTGATTACCCGGTGCGCGGTTAGTTTG‑TAAAGAAGAGTGTCAATATGGGTATAACGATTATTCTGGTGGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCTGGCAATGAAAACGATGGGTGTTAGCGAACTGCGGAGTGTCGATAGGCAGGCACACCTTGAGCCAGTCTCCCGCTGGGTCGCACATGGAGCTGGTGATATTATTGATAATATTAAAGTTTG‑TCCGACac                                    <  2:389553/247‑2 (MQ=255)                                                                                                                                                                                                                 ataaGTTGACAGTGATCACCCTGTTCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACACTGGCTg                              >  1:148497/1‑251 (MQ=255)                                                                                                                                                                                                                                 cACCCTGTTCGCGGTTA‑TTTGGTCAAAAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACACTGGCTGAATCGTTACACGATg               <  2:16731/250‑1 (MQ=255)                                                                                                                                                                                                                                      tgttcaccgTTA‑ATTGGTCAAGAAGAGAGGCAATATGCGTAGAACGATGATTCTGGTCGCACCCGCCAGAGCAGGAACTATTGGGGCAGCGGCGCGGGCACTGAACACGCTGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCAGGCCACCCGCGGGCGCGAACCTGGATCTGGGGATATTATTGAGAAACTTCCAGTTTT‑CCCGCCACTGTCGGCATCGTTACAAGATGTCgctt         >  2:358617/9‑248 (MQ=255)                                                                                                                                                                                                                                        ttCGCGGTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACACTGGCTGAATCGTTACACGATGTCGATTTc       >  1:253462/1‑251 (MQ=255)                                                                                                                                                                                                                                          cgcgGGTA‑TTGGGTCAAGATGAGTGGCAATTGGCGTAAAGCGATTATTCTGGTCGCACCCGGCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGTAAGTCAGGCACACCTTGAGCGAGCCGCCCGCGTTGTCGCATATGGATCTGGTGATATTATGGAGAGTAGTAA‑GTTTGGCCCGACACTGACTGAATCGTTACGCGAGGGCGATTTCa      <  2:541866/250‑1 (MQ=255)                                                                                                                                                                                                                                             ggTTA‑TTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGt  >  1:438484/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTA‑TTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTT‑CCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGT  >  NC_000913/4640679‑4641165

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