Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,309	A→G	L49L (CTA→CTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,309	0	A	G	94.7%	171.0 / ‑4.0	58	L49L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (32/22);  total (35/23)
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.92e-01

AGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAG  >  NC_000913/12062‑12547                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 agtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTCCCTATACCCAGGATGGTGAAACTCTGGtt                                                                                                                                                                                                                                               <  1:533352/251‑1 (MQ=255)      gCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAAGTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCGCGTGTGCTGTAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCa                                                                                                                                                                                                                                         <  1:83223/251‑1 (MQ=255)       cAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTTGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATTCCCAGGATGGGGAAACTCTAGTTGG‑ttag                                                                                                                                                                                                                                        <  1:192684/251‑4 (MQ=255)         aaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGcc                                                                                                                                                                                                                                      >  1:652953/1‑251 (MQ=255)              aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCAAATGATAACCGAATATCTAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACATGGGTACTACCAACTCTTGTGTAGCGAATATGGACGGCACCACTCCTCGTGTGCTGGAGAACGCCGCAGTCGATCGCTCCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCT‑GTTGGGTCAACCGGCt                                                                                                                                                                                                                                  >  2:485011/1‑250 (MQ=255)               tGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGACACTCTGGTTGG‑TCAGCCGGCTa                                                                                                                                                                                                                                 >  1:129902/1‑250 (MQ=255)                           gTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGGGAAGACGTTTAAATGAATAAAATCATTGGTAAAGAACTAGGGACTACCGACTCTTGTGTAGCGATTATGGATGGCACAACTCAGCGTGTGCTGGAGCACACCGAAGGCGATAGCACCACACCTTCTATCATTGCCTATAACCAAGATAGTGAAACTCTGGTTGG‑TCAGCAGGATAAACGtaagacagg                                                                                                                                                                                                                    >  2:150810/1‑243 (MQ=255)                            ttGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTg                                                                                                                                                                                                                  >  2:321035/1‑251 (MQ=255)                            ttGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTATGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTg                                                                                                                                                                                                                  >  1:451867/1‑251 (MQ=255)                                   cAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAAcc                                                                                                                                                                                                           >  1:389693/1‑251 (MQ=255)                                        gTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCa                                                                                                                                                                                                       >  1:430425/1‑250 (MQ=255)                                        gTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATACAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCa                                                                                                                                                                                                       >  1:537419/1‑250 (MQ=255)                                           tCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTTGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCaaaa                                                                                                                                                                                                    <  2:251490/250‑1 (MQ=255)                                            cGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAac                                                                                                                                                                                                   >  1:334932/1‑250 (MQ=255)                                               cccTATTAAAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCATCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAacaccc                                                                                                                                                                                               >  1:307285/1‑249 (MQ=255)                                                            tCACAACCAAATGATGACCGAATATATAGTGGAGACGTTTCGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGCACGCCGCAGGCGATCGCACCACGACTTCTATCATTGCCTCTACCCAGGAGGGTGAAACTCTGGTTAG‑TCAGCCGGCTAAACGTAA‑GCCCGTGCCGAACCCGAAAACCACCATGTTTGCGTTTa                                                                                                                                                                                   >  2:326692/1‑250 (MQ=255)                                                              acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAAc                                                                                                                                                                                >  1:102822/1‑251 (MQ=255)                                                                        gatgaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCaa                                                                                                                                                                                                      >  1:232268/1‑219 (MQ=255)                                                                        gatgaCCGAATATATAGTGGAGACGGTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑GCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCaa                                                                                                                                                                                                      <  2:232268/219‑1 (MQ=255)                                                                            aCCGAATATATAGTGGAGACGTTTAGAGGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGAAGGCACCACTCC‑CGCGTGCTGGAGAACGCCGAAGGCGATAGCACCACGCCTTCTCTCATTGCCTATACCCAGGATGGTGAAAATCTAGTTGG‑TCAGCCGGCTAAACGTCC‑GGCAGTGACGCACCCGCACAAAACTCTGTTTGCGATTAAACGCCTGATTGGTcgcc                                                                                                                                                                 >  2:25964/1‑251 (MQ=255)                                                                                       aGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTCTACCCAGGATGGTGAAACTCTGGTTGG‑TCAACCGGCTAAACGTCA‑GGCAGTGACGAACACGCAAACAACCCAGTTTGCGTTTAAAAGCCTGATTGGCCGCCGCTTCCAggg                                                                                                                                                       >  2:359242/1‑250 (MQ=255)                                                                                             gaCGTTTAGATGGGCAAAATAATTGGTATCGACCTGGGTACTACCCACTATTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGCTGGTAAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGGGACGAACACGCAAAACACCCTGTTTGCGATTAAAAGCCTGAGTggccgacaattacaggcaga                                                                                                                                                    >  2:665322/1‑230 (MQ=255)                                                                                                gTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGt                                                                                                                                              >  1:211786/1‑251 (MQ=255)                                                                                                          gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                                                                                                    <  1:154278/251‑1 (MQ=255)                                                                                                              aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGttt                                                                                                                                >  1:635786/1‑251 (MQ=255)                                                                                                                    tGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatc                                                                                                                           <  1:160228/250‑1 (MQ=255)                                                                                                                     ggTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatca                                                                                                                          <  1:503846/250‑1 (MQ=255)                                                                                                                      gTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATg                                                                                                                        >  1:271009/1‑251 (MQ=255)                                                                                                                             ccTGGGTACTAACAACTCTTGTGTAGCGATTATGGATGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTAA‑GGCAGTGACGAACACGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCAAGGACGAAGAAGAACAGCGTGAAGGTTCCATCATGCCGTTc                                                                                                                  >  1:450239/1‑250 (MQ=255)                                                                                                                              cTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGCCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCa                                                                                                                 >  1:314579/1‑250 (MQ=255)                                                                                                                                gggTACTACAAACTCTTGTGTAGCGTGTATTGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCTATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCaaaa                                                                                                              <  2:113325/251‑1 (MQ=255)                                                                                                                                       acCAACTCTGGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCGGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgc                                                                                                       <  1:665322/251‑1 (MQ=255)                                                                                                                                         cAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctg                                                                                                     >  1:497239/1‑251 (MQ=255)                                                                                                                                          aaCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctgc                                                                                                    >  1:164416/1‑251 (MQ=255)                                                                                                                                              cctgtgTAGAGATTATGGATGGCACCACTCCTCTTGTGCTGGAGAACTCCGAAGTCGAGCGCACCACGCCTGCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGACGGCTAACCGTCA‑TGCAGTGACGAACCTGCAAAACTCCCTGTTTGCGATTAAACGCCTGAGTAGCCGCCGATTCCAGGACGAAGAAGTACAGCGGGATGTTTCCTTCATGCCGTCCAAAATTATTGCTGCTGAt                                                                                                <  2:118647/249‑1 (MQ=255)                                                                                                                                                tgtgTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGa                                                                                                 <  2:102822/248‑1 (MQ=255)                                                                                                                                                 gtgtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTACGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAAAGCGTGATGTTTCCATCCTGCCGTTCAAAATTATTGCTGCTGATAAc                                                                                             >  2:557820/1‑251 (MQ=255)                                                                                                                                                   gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                                                                                           >  2:183598/1‑203 (MQ=255)                                                                                                                                                   gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                                                                                           <  1:183598/203‑1 (MQ=255)                                                                                                                                                    tAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                                                                                          >  2:388563/1‑251 (MQ=255)                                                                                                                                                     aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                                                                                          <  1:321035/250‑1 (MQ=255)                                                                                                                                                               gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                                                                                        <  2:218456/178‑1 (MQ=255)                                                                                                                                                               gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                                                                                        >  1:218456/1‑178 (MQ=255)                                                                                                                                                                        acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa                                                                       >  1:176007/1‑250 (MQ=255)                                                                                                                                                                           ccgcttgtgcttgagtccGCCGAAGGCGGTCGCTCCACGCCTTCTATCATTGCCTATCCCCAGGATGGTGAAACTCTGGTTGT‑TCAGCCGGCTAAACGTCA‑GGCAGTGACTAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAAGGATTGCTGCGGATAACGGCGACGCATGGGTCGAAGTTAAAgg                                                                   <  2:358792/234‑1 (MQ=255)                                                                                                                                                                                 gTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTCAACGCCTGATTGGCCGCCGCTTCCCGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGAAGTTAAGGGCCAGa                                                              >  2:150574/1‑250 (MQ=255)                                                                                                                                                                                     tGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAACATg                                                         >  2:138587/1‑251 (MQ=255)                                                                                                                                                                                      ggAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATg                                                         <  1:557820/250‑1 (MQ=255)                                                                                                                                                                                       gagaACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTa                                                                                                                                             <  1:459871/165‑1 (MQ=255)                                                                                                                                                                                       gagaACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTa                                                                                                                                             >  2:459871/1‑165 (MQ=255)                                                                                                                                                                                          aaCGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAcc                                                    <  1:308964/251‑1 (MQ=255)                                                                                                                                                                                           aCGCCGAAGGCAAACGCACCAAGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACTCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCCccg                                                   <  2:169945/251‑1 (MQ=255)                                                                                                                                                                                                       aTCGCACCACGACTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGAAGAGACGAAACCGCAAAACACCCTGCTTGCGATTAAACGCCTGATTGGCCGACGCTTCCAGGACGAAGAAGTACAGCGTGATGTGTCAATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGAAGCATGGGTCGAAGTTAAAGGCCAGACAATGGCAACGCCGCCGAtattt                                       >  2:28539/1‑247 (MQ=255)                                                                                                                                                                                                               acgctttcttccatgcccTATACCCAGGAGGGTGAACCGCTGGTTGG‑CCAGCGGGATAAACG‑CACGGCAGGGAAGAACCCCAAAAACACCCTGTTTGCGATTAAACGCCTGATGGGCCGCCGCTCCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTCAAAAATTTTTGCTGCGGATACCGGCGCCGCATGGGGCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGt                               <  1:150574/235‑1 (MQ=255)                                                                                                                                                                                                                cGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCATCGCATGGGTCGAATTTAAAGGCCAGAAAATGGAACCGCCGCAGATTTCTGCTGAAggg                              >  2:270393/1‑249 (MQ=255)                                                                                                                                                                                                                  ccTTCTATCATTGCCTAGACCCAGGATGGTGAAACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGAAAAACCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCt                            <  1:388563/251‑1 (MQ=255)                                                                                                                                                                                                                                      ccAGGATGGTGACACTCTGGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGAGCTGAAAAAAATGAAGAAACCCg        >  2:443157/1‑251 (MQ=255)                                                                                                                                                                                                                                            gggaaaaaCTCTGGTTGG‑TCGCCCGGCTAAACGTCG‑GACAGGAACGAACCCGCAAAACCCCCTGTTTGCGATTAAACCCCGTATTGGCCGCCCCTTCAGGGACAAAAAAGTACAGCGTTATTTTTCCATCATGCCGTTAAAAATGATTGCTGTGGAAAACGGCGACGAAGGGGTCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAg  <  1:227508/246‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 AGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGG‑TCAGCCGGCTAAACGTCA‑GGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAG  >  NC_000913/12062‑12547

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