Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,541	G→A	S82S (AGC→AGT)	yaaI ←	DUF2541 domain‑containing protein YaaI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,541	0	G	A	87.2%	79.2 / 5.1	39	S82S (AGC→AGT)	yaaI	DUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (21/13);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACG  >  NC_000913/11312‑11778                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              cTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGggcg                                                                                                                                                                                                                          >  1:238904/1‑251 (MQ=255)       aaaaTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTg                                                                                                                                                                                                                    >  1:301696/1‑251 (MQ=255)        aaaTCCAGATCTGAATCATTCAGGTAAATAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGa                                                                                                                                                                                                                   <  2:338026/251‑1 (MQ=255)             cAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGCTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTTTCGCTGTTAATGTTTATCCAGTCAGGGGTTTGCCCTTTTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATa                                                                                                                                                                                                               >  2:190788/1‑250 (MQ=255)                        cATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGTAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCAcc                                                                                                                                                                                                   <  1:62587/251‑1 (MQ=255)                           tCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGCAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGACGGCATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCt                                                                                                                                                                                                >  1:464684/1‑251 (MQ=255)                                 aaTTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGACGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTg                                                                                                                                                                                          >  2:198659/1‑251 (MQ=255)                                       cattcattTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAAGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATc                                                                                                                                                                                    <  1:235346/251‑1 (MQ=255)                                        attcattTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGCTGCACCGCTTAACTGTAAATCg                                                                                                                                                                                   >  1:434008/1‑251 (MQ=255)                                              tttttGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATg                                                                                                                                                                              <  1:108081/250‑1 (MQ=255)                                                         tCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCagaga                                                                                                                                                                   <  1:156886/250‑1 (MQ=255)                                                                         aGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCTATGAGTTCACCGTATGACCGGAAAAGGTGATTGTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTcac                                                                                                                                                   <  2:227001/250‑1 (MQ=255)                                                                            cgtcgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATg                                                                                                                                               <  1:106425/251‑1 (MQ=255)                                                                             gtcgCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGc                                                                                                                                              >  1:375659/1‑251 (MQ=255)                                                                                   gATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGATGGAATATTAAGACTCTGACTGGCGCAACGGGCGGCTTTGAAATAAACAGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACaa                                                                                                                                         >  2:405988/1‑250 (MQ=255)                                                                                   gATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACaaa                                                                                                                                        <  2:305758/251‑1 (MQ=255)                                                                                    aTAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGGGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAAc                                                                                                                                       >  2:221978/1‑251 (MQ=255)                                                                                                               gAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGt                                                                                                            >  1:30140/1‑251 (MQ=255)                                                                                                                       cGTATGACCGGACAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTAc                                                                                                    >  2:106119/1‑251 (MQ=255)                                                                                                                        gtaggaacggAAAAGGTGATTTTTGAGACGCAGAGTTTATTGTCGTTAACACTATTAATGTTGATCGAAACAGTGGTATGCCCTTATTTTATTTCTGAAGGAATATTCAGACTCTGACTGTCGCTACGGGCGGCATTGAAATAAACCGCTGCACCGCTTAACTTTAAATCGACATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCAACAAACAGGAAGGTTCTGCACCAGATCGTTGGATTCGTTa                                                                                                     >  2:394057/10‑249 (MQ=255)                                                                                                                          aTGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTAcgcg                                                                                                 <  1:190788/251‑1 (MQ=255)                                                                                                                                   aaaGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCa                                                                                         >  1:177024/1‑250 (MQ=255)                                                                                                                                   aaaGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCa                                                                                         >  1:427913/1‑250 (MQ=255)                                                                                                                                     aGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGGCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTGCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATg                                                                                      <  2:327282/251‑1 (MQ=255)                                                                                                                                           tttttGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGCAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACAGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCaagaacgccg                                                                                >  2:174847/1‑243 (MQ=255)                                                                                                                                               tGAGACGCAGCGTTTATTGTCGTTATCACTGATAATGTTGCTCAAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGGTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgaagacgagg                                                                             >  1:416853/1‑241 (MQ=255)                                                                                                                                                       ggCGTTTATGGTGGT‑ATCACTGTTAATGTTTATCCCGTCAGTGGGTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGGAAATCGCCCTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGt                                                                   <  1:118942/250‑1 (MQ=255)                                                                                                                                                                 gTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCt                                                          <  1:106119/251‑1 (MQ=255)                                                                                                                                                                   cGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGTCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAACTCGCCATGGTCGGCAGCGAGTTGTATGCTTTTCACAATGCGACAACCAGGAAGTTTCAGCGCCAGATCGTTGGTTGCGTAATACGGCATTGCATTGACGCCGAGGAGTTTCTGATCGTTTGCCTgc                                                        >  2:371128/1‑251 (MQ=255)                                                                                                                                                                   aggtATCACTGTTAATGTTGATCCAGTAAGTGGTTTGCCATTCTTTTATTTCTGAAGGAATATTCAGACTCTGGCCGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATAGCCATGGTCGGCAGGGAGATGTATGAGTTTCACAATGCACCACACAGGAAGTTTCAGCGCCCGCTTGTTGTTTTCGTTACGCGGCATTGAAATGACGCCGAGGAGATTATGGGAGTTTCCCTgc                                                        >  2:253303/4‑251 (MQ=255)                                                                                                                                                                      tATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgacg                                                                                   >  1:367994/1‑218 (MQ=255)                                                                                                                                                                      tATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACGCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTGTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATgacg                                                                                   <  2:367994/221‑4 (MQ=255)                                                                                                                                                                        tcacTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGt                                                   >  1:248501/4‑251 (MQ=255)                                                                                                                                                                                    ttGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGATTATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGAATTGCAATGACGACGAGGGGTTTATGGGAGTTTGCCTGCGCCGTGCGACACAGCAt                                       >  2:373575/1‑251 (MQ=255)                                                                                                                                                                                                              ttttATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGaa              >  1:290778/1‑250 (MQ=255)                                                                                                                                                                                                                 tATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCg          >  2:259878/1‑251 (MQ=255)                                                                                                                                                                                                                        gAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTaaaa     <  1:107019/249‑1 (MQ=255)                                                                                                                                                                                                                        gAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAAGAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGGTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGTTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTTGCGGAAATCGTAAAAAc   <  2:464684/251‑1 (MQ=255)                                                                                                                                                                                                                         aaGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACg  <  1:221978/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              CTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACG  >  NC_000913/11312‑11778

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: ‑