Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,237	C→T	R25R (CGC→CGT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,237	0	C	T	97.6%	116.4 / ‑4.4	42	R25R (CGC→CGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTA  >  NC_000913/11993‑12478                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  aaaaaTTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTg                                                                                                                                                                                                                                             >  2:118848/1‑251 (MQ=255)          tGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCACCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACCGACTCACAACCACATGCTGACCGAATATATAGTGGAGACGTTTAAATGGGTAAAATAATTGGTATCGACCTGGGTCCTACCCACTCTTGTGAAGCGATTATGGATGGCACCCCTCATCGTGTGCTGGGGAACGcc                                                                                                                                                                                                                                    >  2:308310/1‑251 (MQ=255)             aTCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGa                                                                                                                                                                                                                                  >  2:26822/1‑250 (MQ=255)                    ccTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATc                                                                                                                                                                                                                          <  1:58575/251‑1 (MQ=255)                                               ttAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCGTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCtat                                                                                                                                                                                                <  2:237765/250‑1 (MQ=255)                                               ttAGGAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCtag                                                                                                                                                                                                <  2:248501/250‑2 (MQ=255)                                                tagtagTTAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATCGCCtata                                                                                                                                                                                               <  2:66905/250‑1 (MQ=255)                                                 agtagtCAACCGCAGTGAGTGAGTCTTCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCTCAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCTCCACGCCTTCTATCATTGCCTATAc                                                                                                                                                                                              <  2:349503/250‑1 (MQ=255)                                                           cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTg                                                                                                                                                                                   <  1:26822/251‑1 (MQ=255)                                                                 gagtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCCGGCTGGTGAAAcac                                                                                                                                                                             >  2:296631/1‑249 (MQ=255)                                                                 gagtgagtCTGCAAAAAAATGAAAGTGAGCAGTTGAAACCAGACGTTTCGCCCCTAGTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCCCCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAActc                                                                                                                                                                             <  2:335332/251‑1 (MQ=255)                                                                           gCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTc                                                                                                                                                                    <  1:373575/250‑1 (MQ=255)                                                                             aaaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCCCCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGc                                                                                                                                                                 >  1:82046/1‑251 (MQ=255)                                                                               aaaaaTGAAATTGGGCAGTTGAAATCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCg                                                                                                                                                               >  1:22319/1‑251 (MQ=255)                                                                                aaaaTGAACTTGGGCAGTTGAAACCAGCCGTTTCGCCCCTATTCAAGACTCACAACCACATGATGACCGAATATATAGTGGAGCCGTTTAGATGGTTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTCGCGCTTATGGATGGCAACACTCCTCGTGTGCTGGAGACCGCCGCATGCGCTCGCACCCCGCCTTCTATCATTGACTATACCCAGGATGGTGAAACTCTGGTTGGGCAACCgg                                                                                                                                                              >  2:38890/1‑251 (MQ=255)                                                                                             gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGc                                                                                                                                                 <  2:22319/251‑1 (MQ=255)                                                                                              cagatGAAACCAGACGTTTCGACCATATTCCAGACTCACAAACACATGATGCCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCAGAGTACTACCAACTCTTGTGTAGCGATTATGGGTGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCGTCTATACTTGCCTATACCCCGGATGGAGAAACTCTGGTTGGTCAGCCGGCTAACCGtacggc                                                                                                                                                 >  2:502715/5‑245 (MQ=255)                                                                                                gTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGt                                                                                                                                              <  1:323199/251‑1 (MQ=255)                                                                                                     aaCCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGCAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTCGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACg                                                                                                                                          >  2:234150/1‑250 (MQ=255)                                                                                                         aGACGTCTCGCCCTTATATCAGACTCACAACCACATGATGACCGAATAGATAGTGGAGACGTTTAGACGGGAAAAATAATGGGTATCGACCTGGGTACTGCCAACTCGTGTTTAGCGATTATGGATGGCACCACTCCTCGTGTGCTTGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTTGTTGGTCAGCCGGCTAAACGCCAGGCAGTGACGTAccc                                                                                                                                     <  2:128602/251‑1 (MQ=255)                                                                                                                tCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCaaaa                                                                                                                               <  1:318996/250‑1 (MQ=255)                                                                                                                        aTTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGtt                                                                                                                      >  1:491558/1‑251 (MQ=255)                                                                                                                                   acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGCTTAAAc                                                                                                           >  1:278604/1‑251 (MQ=255)                                                                                                                                   acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGGAGCGATTAGGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAAc                                                                                                           <  2:198251/251‑1 (MQ=255)                                                                                                                                         acaTGCTGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACACCGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGACTAACCGTCAGGCAGGGGCGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTg                                                                                                      >  2:151821/1‑250 (MQ=255)                                                                                                                                                              atgagaCGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCATCTCTTGTGTAGCGATTATGGATGGCACCACTCCGCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTGGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGAGTGGCCGCCGCGTCCAGGACg                                                                                <  2:282808/249‑1 (MQ=255)                                                                                                                                                                        tAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAc                                                                       >  1:92986/1‑250 (MQ=255)                                                                                                                                                                             gggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                 >  1:106687/1‑251 (MQ=255)                                                                                                                                                                              ggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTgg                                                                >  1:37624/1‑250 (MQ=255)                                                                                                                                                                              ggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGa                                                                <  1:25267/251‑1 (MQ=255)                                                                                                                                                                               gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                               <  2:92986/251‑1 (MQ=255)                                                                                                                                                                                   aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTGCCAGGACGAAGAAGTACAGCGTGATGttt                                                           <  1:156198/251‑1 (MQ=255)                                                                                                                                                                                          ggTATCGACCTGGGTACTACCAACTCTTGTTTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACGCCCTGTTGGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatcat                                                    <  2:278604/251‑1 (MQ=255)                                                                                                                                                                                               cGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGt                                               >  1:50247/1‑251 (MQ=255)                                                                                                                                                                                                gACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGtt                                              <  1:376013/251‑1 (MQ=255)                                                                                                                                                                                                                   ctTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCCATCATTGCCTATACCCAGGATGGTTAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGctgctg                             <  1:296631/249‑1 (MQ=255)                                                                                                                                                                                                                     tgtgTAGCGAGTATGGAGGGCACCACTCCTCCTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCGATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGGACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGa                            <  2:204541/248‑1 (MQ=255)                                                                                                                                                                                                                       tgtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGAACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACg                       >  1:37806/1‑251 (MQ=255)                                                                                                                                                                                                                          aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                                     >  1:37835/1‑122 (MQ=255)                                                                                                                                                                                                                          aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                                     <  2:37835/122‑1 (MQ=255)                                                                                                                                                                                                                                   ggatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGCTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATggg           >  2:153494/1‑251 (MQ=255)                                                                                                                                                                                                                                             acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa  >  1:506924/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  AAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTA  >  NC_000913/11993‑12478

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