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	96.0%	138.0 / ‑2.4	50	G74G (GGT→GGC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (27/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

ACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAA  >  NC_000913/12137‑12626                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   aCCGAATATATAGTGGAGACGTTTAGAGGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGAAGGCACCACTCC‑CGCGTGCTGGAGAACGCCGAAGGCGATAGCACCACGCCTTCTCTCATTGCCTATACCCAGGATGGTGAAAATCTAGTTGGTCAGCCGGCTAAACGTCC‑GGCAGTGA‑CGCACCCGCACAAAACTCTG‑TTTGCGATTAAACGCCTGATTGGTcgcc                                                                                                                                                                                                                                                >  2:25964/1‑251 (MQ=255)            aGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTCTACCCAGGATGGTGAAACTCTGGTTGGTCAACCGGCTAAACGTCA‑GGCAGTGA‑CGAACACGCAAACAACCCAG‑TTTGCGTTTAAAAGCCTGATTGGCCGCCGCTTCCAggg                                                                                                                                                                                                                                      >  2:359242/1‑250 (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                                                                                                                                                                                                                          >  2:183598/1‑203 (MQ=255)                                                                        gtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACa                                                                                                                                                                                                                          <  1:183598/203‑1 (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)                                                                                                                                                                                                                            caccctG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACg                        >  1:130886/5‑251 (MQ=255)                                                                                                                                                                                                                                      tGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAAGATGCTCAGc               >  2:180208/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                              >  1:536516/1‑177 (MQ=255)                                                                                                                                                                                                                                                    tGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCCGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaa  >  1:472739/1‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   ACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAA  >  NC_000913/12137‑12626

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