Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,410,859	C→T	L103F (CTC→TTC)	yjfJ →	PspA family protein YjfJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,859	0	C	T	93.1%	154.7 / 3.2	59	L103F (CTC→TTC)	yjfJ	PspA family protein YjfJ
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (29/25);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

TCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGG  >  NC_000913/4410613‑4411105                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       accATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCaa                                                                                                                                                                                                                                                    <  2:159450/250‑2 (MQ=255)      tGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTAcc                                                                                                                                                                                                                                                <  1:347886/250‑1 (MQ=255)        aagaaACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGGCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCt                                                                                                                                                                                                                                             <  2:509418/251‑1 (MQ=255)              cccAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAag                                                                                                                                                                                                                                         >  1:215105/1‑249 (MQ=255)               ccAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGAGCGCCTTGAGAATCTCATTACCGCTGAAGAGc                                                                                                                                                                                                                                      >  2:195682/1‑251 (MQ=255)                 aaggCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCaa                                                                                                                                                                                                                                    <  2:201599/249‑1 (MQ=255)                   ggCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAgt                                                                                                                                                                                                                                  <  2:13819/251‑1 (MQ=255)                         cgaaTGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTAACGCTGAAGAGCAAGTGTTGTc                                                                                                                                                                                                                            <  2:133958/248‑1 (MQ=255)                          gcATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCg                                                                                                                                                                                                                           >  1:245241/1‑251 (MQ=255)                           cATGCTGGACCAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCCTTACCGCTGACGAGCAAGTGTTGTCg                                                                                                                                                                                                                           >  2:466354/1‑250 (MQ=255)                            aTGCTGGAACAGCATATTCCTGACGCTCAAGCTGAACTCGGTCAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAATCACGATAAGCTGAAAGCTTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGACGTTGCTGAAGAAATCTCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGa                                                                                                                                                                                                                          >  1:486105/1‑250 (MQ=255)                                    aCAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGTAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGGAGAGCAAGTGTTGTCGAATCTGGagg                                                                                                                                                                                                                 >  2:257424/1‑249 (MQ=255)                                           atTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCTGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCg                                                                                                                                                                                                          <  1:484943/251‑1 (MQ=255)                                                tGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATg                                                                                                                                                                                                     >  1:342075/1‑251 (MQ=255)                                                     cTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTg                                                                                                                                                                                                >  1:243159/1‑251 (MQ=255)                                                              aaCTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATGAAGGCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGGTGAAGAAATCGCGCGCCTGGAGAATTACATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCa                                                                                                                                                                                       <  2:235563/251‑1 (MQ=255)                                                                       aaGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCCAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTacagcg                                                                                                                                                                              >  1:608495/1‑251 (MQ=255)                                                                          ccGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGCAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCt                                                                                                                                                                                                                      >  1:108477/1‑208 (MQ=255)                                                                          ccGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCt                                                                                                                                                                                                                      <  2:108477/208‑1 (MQ=255)                                                                           cGGAAAATCTCGCGTTGATCTGCTGGCGCGGTTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGCAGCGTTAAGCGAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGACGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAGAAAGCAGTTgcagcgcaa                                                                                                                                                                           >  2:609048/1‑247 (MQ=255)                                                                                aaTCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGCTGGCGTGGAAAAAGCAGTTACAGCGAcagcgc                                                                                                                                                                       >  2:109341/1‑249 (MQ=255)                                                                                 aTCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTTTCCCGTGATGGCGTGGAAAAACAAGGTACAGCGTcagcgcagg                                                                                                                                                                    >  2:112875/1‑250 (MQ=255)                                                                                     cgcgTTGATCTGCTGGCGCGGGTGAAATTAAGGCACGATAATCTGAAAGATTTACGTGAGCGCAAAGCCATTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAAGCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTTCAGCGAcagcgcagcg                                                                                                                                                                   <  2:608495/248‑1 (MQ=255)                                                                                                    gcgcGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGagca                                                                                                                                                 <  1:292676/251‑1 (MQ=255)                                                                                                        gggTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCaa                                                                                                                                             <  2:203281/251‑1 (MQ=255)                                                                                                                  aaGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGAGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCGCAGTTTGAGCAGCAAATGGAAGTCg                                                                                                                                   >  2:230538/1‑251 (MQ=255)                                                                                                                          aTAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGGAGCCCGTGCGCTGGCAGCGTTAAGCAAGAACGTTAATCAGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGAGACAGCGCAGCGTATCGCTCAGTTTGGGCAGCCAATGGAAGACGGTAAAGcc                                                                                                                           >  2:626038/1‑251 (MQ=255)                                                                                                                            aaGCTGAAAGATTTACGTGAGCGCAAAGCCAGACTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCAc                                                                                                                         >  1:504760/1‑251 (MQ=255)                                                                                                                              gCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTAGCGCTCAGTTTGAGCAGCAAATGGAGATCGTTAAAGCCACTg                                                                                                                       >  2:142280/1‑251 (MQ=255)                                                                                                                                  aaaGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGGAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGCATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGAAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGc                                                                                                                   >  2:545763/1‑251 (MQ=255)                                                                                                                                  aaaGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGc                                                                                                                   >  1:674141/1‑251 (MQ=255)                                                                                                                                      aTTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATg                                                                                                               >  1:527088/1‑251 (MQ=255)                                                                                                                                        ttACGTGAGCGCAAAGCCAGACTGGAAGCCCTTGCGCTGGAAGCGTTAAGCACGAGCGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATATCATTACCGCTGCAGAGCAAGCGTTGTCGAATCTGGAAGTATCCAGTGATGGCGTGGAAGAAGCAGTTACAGCGACAGCGCAGCGTATCGCTAAGTTTGATCAGCAAATGGAAGTCGTTAAAAACACTGGAGCCATGc                                                                                                              >  2:666927/1‑250 (MQ=255)                                                                                                                                           cGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTAGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCGACGCATCACTCAGTTTGCGCAGCAAATGGCAGTCGTTAAAGCCACTGAAGCCATGCAGCg                                                                                                          >  2:215040/1‑251 (MQ=255)                                                                                                                                            gTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGt                                                                                                         <  1:195682/251‑1 (MQ=255)                                                                                                                                              gAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGc                                                                                                       >  1:19326/1‑251 (MQ=255)                                                                                                                                                   cAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCacaac                                                                                                  >  1:312541/1‑251 (MQ=255)                                                                                                                                                      aGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAgg                                                                                               <  2:342075/251‑1 (MQ=255)                                                                                                                                                       gCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGc                                                                                              <  1:646280/251‑1 (MQ=255)                                                                                                                                                                    cccGGGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTGACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                                                                 <  2:616073/251‑1 (MQ=255)                                                                                                                                                                     ccGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                                                                <  1:148966/251‑1 (MQ=255)                                                                                                                                                                                taGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTAGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCCGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAAAGGAAGTCGTTCAAGCCACTGAAGCCAAGACGCGTGCACAACAGGCTGTAACAAACACTACCGTTGgcgcc                                                                     >  2:360139/2‑250 (MQ=255)                                                                                                                                                                                 atcggtAAGCAAGAACTGTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCTCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAAGGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCAt                                                                    <  2:10149/246‑1 (MQ=255)                                                                                                                                                                                 aGCGTTAAGTAAGCACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCAAGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCAt                                                                    <  2:245241/251‑1 (MQ=255)                                                                                                                                                                                   cGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                                                  >  1:679104/1‑251 (MQ=255)                                                                                                                                                                                      tAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGAAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                                                               >  1:266265/1‑251 (MQ=255)                                                                                                                                                                                         gCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                                                            >  1:18992/1‑251 (MQ=255)                                                                                                                                                                                               aCGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCg                                                      >  1:545424/1‑251 (MQ=255)                                                                                                                                                                                                        cgtcgtTGATTAACGAAGTGGCGGAAGAAATCGCGCGCCTTGAGAAGTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCCCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCa                                             <  2:80260/251‑1 (MQ=255)                                                                                                                                                                                                           cgtTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGa                                           <  1:197627/250‑1 (MQ=255)                                                                                                                                                                                                                aTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGACGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCAACAGCGGCAGacccctt                                     >  2:307385/1‑245 (MQ=255)                                                                                                                                                                                                                    aCGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGGTGTCGAATCTGGAAGGCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTaaaa                                 <  1:662257/251‑1 (MQ=255)                                                                                                                                                                                                                                 agaaATCGCGCGCCTTGAGCATTTCATTACCGCTGAAGAGCAAGTGTGGTCGAATCTGTAAGTCTCCCGTGATTGCGTTGAAAAAGCAGTTCCAGCGACAGCGCAGCGTATCGCACAGTTTGAGCAGCAAATGGAAGTCGTTAAATCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAACTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAACCTCCTGCAAAcgc                    <  2:586581/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGATGTCGAATCTGGACGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACATCGGCAGAAACATTAAGACGCCAGCAACcgcgcc                 >  2:24510/1‑249 (MQ=255)                                                                                                                                                                                                                                            gcCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAAAAACCTTAAAACGCCCGCAACCGCGT‑AGGcccca         >  2:3774/1‑247 (MQ=255)                                                                                                                                                                                                                                                tGAGAATTTCATGACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGGATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACg       <  1:52619/249‑1 (MQ=255)                                                                                                                                                                                                                                                  agaATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCGCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAGACGCCTGCAAACGCGTCAGGCCGAACGTCAg   <  1:24510/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gaATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAgg  >  1:553090/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gaATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAgg  <  1:109341/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       TCCATTGAAGAAACCCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACGTCAGG  >  NC_000913/4410613‑4411105

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