Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,202	A→C	T4T (ACT→ACG)	yjcH ←	DUF485 domain‑containing inner membrane protein YjcH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,202	0	A	C	98.3%	193.4 / NA	60	T4T (ACT→ACG)	yjcH	DUF485 domain‑containing inner membrane protein YjcH
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (38/21);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCT  >  NC_000913/4284963‑4285443                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    cAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATtcat                                                                                                                                                                                                                                          <  1:317914/251‑1 (MQ=255)  aGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATtcatc                                                                                                                                                                                                                                         >  1:284763/1‑251 (MQ=255)       tAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCCACACCACTTGGCATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGCTCCGTAAAATAAAGCCGCTTTGAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATcac                                                                                                                                                                                                                                       >  1:465970/1‑248 (MQ=255)              cGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATTGTGGCAAACCGTTGCCGGTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTATATCCACTGATAAATCGTACCATTCATCACAGATTCTcc                                                                                                                                                                                                                              >  2:250388/1‑250 (MQ=255)              cGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTcc                                                                                                                                                                                                                              >  1:107429/1‑250 (MQ=255)                    gaACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACCGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGt                                                                                                                                                                                                                        >  1:284398/1‑250 (MQ=255)                        cacggAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCCGCCGGGCGCGAAGTCGATCAGTAAAATAAAGCCGATATAAACTGACAGCATAATAATCGACCGGATGGTGACAAACCGTTGCCCTTTTTCGACTAACTCCCTGAACTGCGCATTGTACTCTATTCGCTGATAAATCGTACAATTCATCACAGCATCTCCAGAGGTAATg                                                                                                                                                                                                                    >  2:248936/4‑250 (MQ=255)                        aaaGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGt                                                                                                                                                                                                                   >  1:90166/1‑251 (MQ=255)                         aaGGAGATCACAATCACTCAAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGCTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTTTATCCGCTGATAAATCGTACCATTCATAACAGATTATCCAGAGGAAATga                                                                                                                                                                                                                   >  2:38988/1‑249 (MQ=255)                            gagaTCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCAGTCATCACAGATTCTCCAGAGGTAATGTAggg                                                                                                                                                                                                               <  1:202579/251‑1 (MQ=255)                                      gcccTCCAACACCAAGTGGAATACCGCGTGGGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTaa                                                                                                                                                                                                     <  2:90166/248‑1 (MQ=255)                                            cAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGAAAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCg                                                                                                                                                                                               >  1:300353/1‑251 (MQ=255)                                                              cgcgTGTGACGCGGGTGTTCGTATTCAGCGGCGTGCCCAGCCAGCCGGGCTCGAAGGCGATCAGTAAAATAAAGCCTATATAAACTGCCAGCATAATAATCGACAGTATGGTGGCAAACCGTTGCCGTTTTGCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCGGATAAATCGTACCATTCATCACATATTCTCCAGAGGTAAGGTAGGGACTGTTTTAATTCCCGCTCCCTGATGGGAGAAgg                                                                                                                                                                             <  2:272917/251‑1 (MQ=255)                                                                   gtgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAc                                                                                                                                                                        >  1:154148/1‑251 (MQ=255)                                                                        gCTGGTGTTCGGATACAGCGGCGTGCCCCGCCAGCCGTGCGCGAAGGCGATCAGTAAAATAATGCCGATATAAACTGTCAGCATAATAATCGGCAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCTCATTGTGTTCTATCCGCTGATAAATCGTACCATTCATCACATATTCTCCAGAGTTAAGGTAGGGATTGTGGTAATTCCCGCTCCCGTATTGAAGAAGGTGAACGCTCg                                                                                                                                                                   <  2:320276/251‑1 (MQ=255)                                                                          tGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCACACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGCTAAATCGTACCATTCATCCCCGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGGTACCGCTCggg                                                                                                                                                                 >  2:98674/1‑251 (MQ=255)                                                                          gggTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCggg                                                                                                                                                                 <  2:154148/250‑1 (MQ=255)                                                                                                   cagccagcCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCgg                                                                                                                                         >  2:141309/1‑250 (MQ=255)                                                                                                                    aGGCGATCAGTAAAAGAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGACGAATGTAGGCCGGATAAGGCG‑TTATAcgccg                                                                                                                      <  2:257988/251‑1 (MQ=255)                                                                                                                     ggCGATCAGTAAAATAAAGCCGCTATAAACTGCCAGCATAATAATCGACAGGATGGAGGCAAACCGTTGCAGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCACTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATAcgccgc                                                                                                                     >  2:453887/1‑251 (MQ=255)                                                                                                                     ggCGATCAGTAAAATAAACCCGATATAAACTGCCAGCCTAATCATCGCCAGGATGGTGGCAAACCGCTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCCGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTATCCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATAcgccgc                                                                                                                     >  2:444964/1‑251 (MQ=255)                                                                                                                       cGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCAt                                                                                                                   >  1:99643/1‑251 (MQ=255)                                                                                                                         aTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCagcc                                                                                                                 <  2:75365/251‑4 (MQ=255)                                                                                                                            aGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTACTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGCATGTAGGCCGGATAAGGCG‑TTATACGCAGCCTCCGGc                                                                                                              >  2:330929/1‑251 (MQ=255)                                                                                                                             gTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCTACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACGCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGca                                                                                                             <  2:284398/251‑1 (MQ=255)                                                                                                                             gTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGc                                                                                                              >  1:96625/1‑250 (MQ=255)                                                                                                                                aaaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAATCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGcaatc                                                                                                          <  1:336277/251‑1 (MQ=255)                                                                                                                                     aaGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGc                                                                                                     >  2:149829/1‑251 (MQ=255)                                                                                                                                           atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATg                                                                                               >  1:284643/1‑251 (MQ=255)                                                                                                                                             ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATTGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTAGACGCCGCATCCGGCAATCAATGCCTGATGc                                                                                              <  2:443750/250‑1 (MQ=255)                                                                                                                                             ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACAGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATCAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCg                                                                                             >  1:467275/1‑251 (MQ=255)                                                                                                                                               aaaCTGCCAGCATAATAAGCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTATGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                                           <  1:473483/251‑1 (MQ=255)                                                                                                                                                  cTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCCATCAATGCCTGCTGCGACGCt                                                                                        >  2:209540/1‑251 (MQ=255)                                                                                                                                                    gCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                      >  2:390377/1‑251 (MQ=255)                                                                                                                                                         cATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGGTA‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                                                 >  2:369551/1‑251 (MQ=255)                                                                                                                                                           taataaTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                               >  1:239084/1‑251 (MQ=255)                                                                                                                                                                   gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCTCATTGTCTTCTATCCGCTGATAAATCGTACCAGTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATGCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCTTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                                       <  2:325745/251‑1 (MQ=255)                                                                                                                                                                      aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                     >  1:131871/1‑250 (MQ=255)                                                                                                                                                                      aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                                     >  2:440296/1‑250 (MQ=255)                                                                                                                                                                        gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                   >  1:148106/1‑250 (MQ=255)                                                                                                                                                                              gtcAAACCGTTGCCGTTTGCCGACTAACTCCCTGAAATGCGCAGTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCCCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCGGATGCGACGCTGTCGCGTCTTATCAGGCCTTCAAACCg                                                             <  1:254342/248‑1 (MQ=255)                                                                                                                                                                                      gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACTCCTCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                           <  2:314178/228‑1 (MQ=255)                                                                                                                                                                                      gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                           >  1:314178/1‑228 (MQ=255)                                                                                                                                                                                       ttGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                      >  2:154910/1‑232 (MQ=255)                                                                                                                                                                                       ttGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                      <  1:154910/232‑1 (MQ=255)                                                                                                                                                                                        tGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTACCATTCATCAAAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCAGCTCACTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCAGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCg                                                  >  1:457631/1‑251 (MQ=39)                                                                                                                                                                                              ttttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCGCAGATTCTCCAGAGGTACTGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGcctccg                                            >  2:435750/1‑246 (MQ=39)                                                                                                                                                                                                  cGACTAACTCCCTGAAATGTGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACCGCTTCTCCAGAGGTAATGTCGGGCTTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACAGACTCGCATCCGGCa                                         >  2:434770/1‑250 (MQ=38)                                                                                                                                                                                                  cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCa                                         <  1:253701/250‑1 (MQ=255)                                                                                                                                                                                                        aCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAAGCGTACCATTCATCACAGATTCTCCAGAGGGAATGTAGGGATTGGTTTAATTCCCGCTCCCTTAAGGGAGAAGGTTAACGCTCGGGTAACCCTGGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGAGGCGACGCTGTCGCGTCTTATCAGTCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTgg                                  <  2:300353/251‑1 (MQ=255)                                                                                                                                                                                                           cccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTGAATGCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAACCGTTACCGACTCGCATCCGGGAATTGTggg                                 <  2:465970/249‑1 (MQ=255)                                                                                                                                                                                                           cccTGAAATGCGCATTGTCGTCTATCCGCTGATAAATCGTACCATTCAACACAGATTCTCCAGAGGTAATTTAGGGATTGTTTTAATTCCCGCCCCCTGATGGGAGAAGGGTAACGCTCGGGTATCCCGGGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCAGCCGGCAATCAAGGCCTGATTCGACGCTGTCGTGTCTTATCAGGCCTACAAACCGTTACCGACTCGCACCCGGCAATTGGGGGtt                               <  2:362542/251‑1 (MQ=255)                                                                                                                                                                                                             cTGAAATGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAAcc                                                                                                                                                            >  1:469338/1‑125 (MQ=255)                                                                                                                                                                                                             cTGAAAGGCGCATTGTCTTCTATCCGCTGATAAATCGTACCATTCAGCACAGATTCTCCAGAGGTAATGTAGGGATTGTTGTAATTCCCGCTCCCTTATGGGAGAAGGTGACCGCTCGGGTAAcc                                                                                                                                                            <  2:469338/125‑1 (MQ=255)                                                                                                                                                                                                                        aTTGTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCCATCAATGCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATcgc                  >  2:67666/1‑251 (MQ=255)                                                                                                                                                                                                                           gTCTTCTATCCGCTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGAt               >  1:360992/1‑251 (MQ=255)                                                                                                                                                                                                                                ctATCCACTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCt          >  1:84912/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATACCCTGCTtctc   >  2:256339/1‑251 (MQ=255)                                                                                                                                                                                                                                        tGATAAATCGTACCATTCATCACAGATTCTCCCGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct  >  1:281521/1‑251 (MQ=255)                                                                                                                                                                                                                                        tGATAAATCGTACCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TTATACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCCGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTctct  >  2:30172/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCG‑TT‑TACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCT  >  NC_000913/4284963‑4285443

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