Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,384	T→C	G74G (GGT→GGC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,384	0	T	C	89.7%	93.1 / 2.4	40	G74G (GGT→GGC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (19/16);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

TGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAA  >  NC_000913/12150‑12627                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    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)                                                                      ggatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGCTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATggg                                                                                                                                                                >  2:153494/1‑251 (MQ=255)                                                                                acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa                                                                                                                                                       >  1:506924/1‑250 (MQ=255)                                                                                      cgtgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGcc                                                                                                                                                 <  1:217678/247‑1 (MQ=255)                                                                                                        gAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGCGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAACCACCCTGTTTGCGATTAAACGCCTGATTGGACGCCGCTTCCAGGACGAAGAAGAACAGCGTGATGTTTCCATCATGCCGTTCAACATTATTGCTGCTGATAACGGCGACGCATAGGTCGAAGTTAAAGGCCAGGAAATggacacgccgc                                                                                                                               >  2:48052/1‑241 (MQ=255)                                                                                                               cacgacccacCCCTTCTAACATTGCCTATACCCGGGATGGTGAACCCCTGGTGGGTCAGCCGGCTAAACGTCAGGCGGTGACGAACCCGCAAACCACCCTGTTTGCGATAAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCAAGAAAATGGCACCGCCGCAGATctct                                                                                                                       <  2:106687/245‑1 (MQ=255)                                                                                                                    accacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTctgctg                                                                                                                   >  1:251452/1‑250 (MQ=255)                                                                                                                             tCTATCATTCCCTATCCCCAGGATGGTGAACCTCTAGTGGGTCAGCCGGCAAAACGTCAGGCAGTGACGACCCCGCAAAACCCTCTGTTTGCGATTAACCGCCTGTTTGGCCGCCGTTTCCAGGACGAAGAAGTACAGCGTGATTTTTCCATCATGCCGTTAAAAATTATTTCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGa                                                                                                          <  1:234150/250‑1 (MQ=255)                                                                                                                                   agtGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGaagaag                                                                                                                                                                                                                               <  2:242261/125‑1 (MQ=255)                                                                                                                                   aTTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGaagaag                                                                                                                                                                                                                               >  1:242261/1‑127 (MQ=255)                                                                                                                                        cTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATga                                                                                                 >  1:187465/1‑248 (MQ=255)                                                                                                                                                     ggTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGCTCTCTGCTGGAGTGCTGAAAAAAATGAAGAAAACCGCTGa                                                                                    >  2:117177/1‑248 (MQ=255)                                                                                                                                                         aaaCTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTTATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGACGTTAAAGGCCAGAAACTGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAc                                                                             >  1:185339/1‑251 (MQ=255)                                                                                                                                                              ctggTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATgg                                                                                                                                                                 >  2:232390/4‑162 (MQ=255)                                                                                                                                                              ctggTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATgg                                                                                                                                                                 <  1:232390/159‑1 (MQ=255)                                                                                                                                                                ggtgggTCAGCCGGCTAACCGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGCCCGCCGCTTCCAGGACGAAAAATTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGt                                                                       <  1:218885/246‑1 (MQ=255)                                                                                                                                                                ggTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGAAGTTAAAGGCCAGAAAATGGAACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGACGATTACCTGGGTg                                                                      >  2:436727/2‑251 (MQ=255)                                                                                                                                                                   tGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                   >  2:79492/1‑251 (MQ=255)                                                                                                                                                                                           aGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACATGGGTGAACCGGTAACTGAAGCTGTTATCACCg                                           >  2:22173/1‑251 (MQ=255)                                                                                                                                                                                                    cccGAAAAACCCCCTGTTTGCGATTAACCGCCTTATTGGCCGCTGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCGTTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAATTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAt                                  <  1:308310/251‑1 (MQ=255)                                                                                                                                                                                                     ccGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTCCAGCGAGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGACAAACATGACGAAAACCGCTGAAGATTACCTGGGTAAACCGGTAACTGAAGATGTTATCCCCGTACCgacca                                  >  2:159682/1‑246 (MQ=255)                                                                                                                                                                                                            accccctGTTTGCGATTAAACGCCTGATTGCCCGCCGCTCCCAGGACGAAGAAGTACAGCGTGATTTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAATTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAc                          <  1:502702/246‑1 (MQ=255)                                                                                                                                                                                                               ccctGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCGGGACGAAGCAGTACAGCGTGATGTTTCCATCATGCCGTTCACAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGCAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATCCTTTAACGAt                       >  2:427324/3‑251 (MQ=255)                                                                                                                                                                                                                   gTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTc                   >  1:497469/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaCGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGc         >  1:488404/1‑251 (MQ=255)                                                                                                                                                                                                                                   cTGATTGGCCGCCGCTTCCAGGACGAAGAAGTCCAGCGAGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTCCTGAAAAAAATGACGAAACCCGCTGAAAATTACCTGGGGGAACCGGTAACTGAAGCTGTTATCACCATACCGGCATACTTTCACGATGCTCAGCGGCAGGCAACCa    >  2:321609/1‑250 (MQ=255)                                                                                                                                                                                                                                    tGATTGGCCGCAGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCAATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAAcaaa   >  2:443461/1‑247 (MQ=255)                                                                                                                                                                                                                                    gtaTTGGCCGCCGCTTCCAGGACGAAAAAGTACGGCTGGATGTTTCAATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGCCGCATGGGTCAAATTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaaa  <  1:253051/249‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    TGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAA  >  NC_000913/12150‑12627

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