Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,357	T→C	T65T (ACT→ACC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,357	0	T	C	93.9%	133.8 / ‑2.0	49	T65T (ACT→ACC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (26/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGC  >  NC_000913/12121‑12593                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              tCACAACCAAATGATGACCGAATATATAGTGGAGACGTTTCGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGCACGCCGCAGGCGATCGCACCACGACTTCTATCATTGCCTCTACCCAGGAGGGTGAAACTCTGGTTAGTCAGCCGGCTAAACGTAA‑GCCCGTGC‑CGAACCCGAAAACCACCATG‑TTTGCGTTTa                                                                                                                                                                                                                                 >  2:326692/1‑250 (MQ=255)   acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAAc                                                                                                                                                                                                                              >  1:102822/1‑251 (MQ=255)                 aCCGAATATATAGTGGAGACGTTTAGAGGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGAAGGCACCACTCC‑CGCGTGCTGGAGAACGCCGAAGGCGATAGCACCACGCCTTCTCTCATTGCCTATACCCAGGATGGTGAAAATCTAGTTGGTCAGCCGGCTAAACGTCC‑GGCAGTGA‑CGCACCCGCACAAAACTCTG‑TTTGCGATTAAACGCCTGATTGGTcgcc                                                                                                                                                                                                               >  2:25964/1‑251 (MQ=255)                            aGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTCTACCCAGGATGGTGAAACTCTGGTTGGTCAACCGGCTAAACGTCA‑GGCAGTGA‑CGAACACGCAAACAACCCAG‑TTTGCGTTTAAAAGCCTGATTGGCCGCCGCTTCCAggg                                                                                                                                                                                                     >  2:359242/1‑250 (MQ=255)                                  gaCGTTTAGATGGGCAAAATAATTGGTATCGACCTGGGTACTACCCACTATTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGCTGGTAAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGGGA‑CGAACACGCAAAACACCCTG‑TTTGCGATTAAAAGCCTGAGTggccgacaattacaggcaga                                                                                                                                                                                                  >  2:665322/1‑230 (MQ=255)                                     gTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGt                                                                                                                                                                                            >  1:211786/1‑251 (MQ=255)                                               gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                                                                                                                                                  <  1:154278/251‑1 (MQ=255)                                                   aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGttt                                                                                                                                                                              >  1:635786/1‑251 (MQ=255)                                                         tGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatc                                                                                                                                                                         <  1:160228/250‑1 (MQ=255)                                                          ggTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatca                                                                                                                                                                        <  1:503846/250‑1 (MQ=255)                                                           gTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATg                                                                                                                                                                      >  1:271009/1‑251 (MQ=255)                                                                  ccTGGGTACTAACAACTCTTGTGTAGCGATTATGGATGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTAA‑GGCAGTGA‑CGAACACGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCAAGGACGAAGAAGAACAGCGTGAAGGTTCCATCATGCCGTTc                                                                                                                                                                >  1:450239/1‑250 (MQ=255)                                                                   cTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGCCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCa                                                                                                                                                               >  1:314579/1‑250 (MQ=255)                                                                     gggTACTACAAACTCTTGTGTAGCGTGTATTGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCTATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCaaaa                                                                                                                                                            <  2:113325/251‑1 (MQ=255)                                                                            acCAACTCTGGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCGGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgc                                                                                                                                                     <  1:665322/251‑1 (MQ=255)                                                                              cAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctg                                                                                                                                                   >  1:497239/1‑251 (MQ=255)                                                                               aaCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATtgctgc                                                                                                                                                  >  1:164416/1‑251 (MQ=255)                                                                                   cctgtgTAGAGATTATGGATGGCACCACTCCTCTTGTGCTGGAGAACTCCGAAGTCGAGCGCACCACGCCTGCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGACGGCTAACCGTCA‑TGCAGTGA‑CGAACCTGCAAAACTCCCTG‑TTTGCGATTAAACGCCTGAGTAGCCGCCGATTCCAGGACGAAGAAGTACAGCGGGATGTTTCCTTCATGCCGTCCAAAATTATTGCTGCTGAt                                                                                                                                              <  2:118647/249‑1 (MQ=255)                                                                                     tgtgTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGa                                                                                                                                               <  2:102822/248‑1 (MQ=255)                                                                                      gtgtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTACGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAAAGCGTGATGTTTCCATCCTGCCGTTCAAAATTATTGCTGCTGATAAc                                                                                                                                           >  2:557820/1‑251 (MQ=255)                                                                                        gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                                                                                                                                         <  1:183598/203‑1 (MQ=255)                                                                                        gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                                                                                                                                         >  2:183598/1‑203 (MQ=255)                                                                                         tAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                                                                                                                                        >  2:388563/1‑251 (MQ=255)                                                                                          aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGc                                                                                                                                        <  1:321035/250‑1 (MQ=255)                                                                                                    gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                                                                                                                                      <  2:218456/178‑1 (MQ=255)                                                                                                    gatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAgg                                                                                                                                                                                                      >  1:218456/1‑178 (MQ=255)                                                                                                             acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa                                                                                                                     >  1:176007/1‑250 (MQ=255)                                                                                                                ccgcttgtgcttgagtccGCCGAAGGCGGTCGCTCCACGCCTTCTATCATTGCCTATCCCCAGGATGGTGAAACTCTGGTTGTTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CTAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAAGGATTGCTGCGGATAACGGCGACGCATGGGTCGAAGTTAAAgg                                                                                                                 <  2:358792/234‑1 (MQ=255)                                                                                                                      gTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTCAACGCCTGATTGGCCGCCGCTTCCCGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGAAGTTAAGGGCCAGa                                                                                                            >  2:150574/1‑250 (MQ=255)                                                                                                                          tGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAACATg                                                                                                       >  2:138587/1‑251 (MQ=255)                                                                                                                           ggAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATg                                                                                                       <  1:557820/250‑1 (MQ=255)                                                                                                                            gagaACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTa                                                                                                                                                                                           >  2:459871/1‑165 (MQ=255)                                                                                                                            gagaACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTa                                                                                                                                                                                           <  1:459871/165‑1 (MQ=255)                                                                                                                               aaCGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAcc                                                                                                  <  1:308964/251‑1 (MQ=255)                                                                                                                                aCGCCGAAGGCAAACGCACCAAGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACTCCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCCccg                                                                                                 <  2:169945/251‑1 (MQ=255)                                                                                                                                            aTCGCACCACGACTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGAAGAGA‑CGAAACCGCAAAACACCCTG‑CTTGCGATTAAACGCCTGATTGGCCGACGCTTCCAGGACGAAGAAGTACAGCGTGATGTGTCAATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGAAGCATGGGTCGAAGTTAAAGGCCAGACAATGGCAACGCCGCCGAtattt                                                                                     >  2:28539/1‑247 (MQ=255)                                                                                                                                                    acgctttcttccatgcccTATACCCAGGAGGGTGAACCGCTGGTTGGCCAGCGGGATAAACG‑CACGGCAGGGA‑AGAACCCCAAAAACACCCTG‑TTTGCGATTAAACGCCTGATGGGCCGCCGCTCCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTCAAAAATTTTTGCTGCGGATACCGGCGCCGCATGGGGCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGt                                                                             <  1:150574/235‑1 (MQ=255)                                                                                                                                                     cGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCATCGCATGGGTCGAATTTAAAGGCCAGAAAATGGAACCGCCGCAGATTTCTGCTGAAggg                                                                            >  2:270393/1‑249 (MQ=255)                                                                                                                                                       ccTTCTATCATTGCCTAGACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGAAAAACCCCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCt                                                                          <  1:388563/251‑1 (MQ=255)                                                                                                                                                                           ccAGGATGGTGACACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGAGCTGAAAAAAATGAAGAAACCCg                                                      >  2:443157/1‑251 (MQ=255)                                                                                                                                                                                 gggaaaaaCTCTGGTTGGTCGCCCGGCTAAACGTCG‑GACAGGAA‑CGAACCCGCAAAACCCCCTG‑TTTGCGATTAAACCCCGTATTGGCCGCCCCTTCAGGGACAAAAAAGTACAGCGTTATTTTTCCATCATGCCGTTAAAAATGATTGCTGTGGAAAACGGCGACGAAGGGGTCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAg                                                <  1:227508/246‑1 (MQ=255)                                                                                                                                                                                            tggTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGt                                     >  1:129539/3‑251 (MQ=255)                                                                                                                                                                                             gtttGGCCAGCCGGCTAAACGTCA‑GGCAG‑GAGCCACCCCTAAAAACCCCCTG‑TTTGCGATTAACCGCCTGATGGGCCGCCGCTTCCAGAACGAAAAAGTACAGCGGGATGTTTCCATCATGCCGGTAAAATTTATTGCTGCTGATACCGGCGCCGCTTGGGCCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTg                                    <  1:326692/249‑1 (MQ=255)                                                                                                                                                                                                        cGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACt                         >  1:88220/1‑251 (MQ=255)                                                                                                                                                                                                            aaaaCGTGA‑GGCAGGGA‑CGACCCCGAAAAACACCCTG‑TTTGCGATAAAACGCCTGATTGGCCGCCTCTTCCAGGACGAAAAAGTACAGCGTGATTTTGCCATCATGCCTTTCAAAATTATTGCTGCTGATAACGGCTACGCATGGGTCGAAGTAAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAg                     <  1:627679/250‑1 (MQ=255)                                                                                                                                                                                                             aaaCGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGc                    >  1:678034/1‑251 (MQ=255)                                                                                                                                                                                                                 gTCA‑GGCAGTGA‑CGAACCCG‑AAAACCCTTTGGTTGGCGATAAACCGCCTGTTTGGCCGCCGCTTCAAGGACGAAGAAGTACAGGGTTATGTTTCAAACATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGAATGGTTCAAAGTAAAAGGCCAAAAAATGGCTCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGtt                <  1:270393/251‑1 (MQ=255)                                                                                                                                                                                                                            gA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTAc      >  2:314256/1‑251 (MQ=255)                                                                                                                                                                                                                                  aCCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATAACCGTACCGGc  >  2:280995/1‑250 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              TCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGC  >  NC_000913/12121‑12593

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