Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,137,212	C→T	R313R (CGC→CGT)	flgI →	flagellar P‑ring protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,137,212	0	C	T	95.7%	131.7 / ‑3.3	47	R313R (CGC→CGT)	flgI	flagellar P‑ring protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (12/33);  total (12/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGCCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAGCGGGCGAAGATC  >  NC_000913/1136981‑1137455                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         tgtcGATAGTCAGAAAAGGCAGGTTAATGTCGCCCCGCAGGACGCTAAAGTAGTGATTAGCTGGCGCACCGGTTCGGTGGTGAGGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAAGCTCTCAGTAACAGTTAGTCGTCAGGCCAATGTCAGCCAGCCAGATGCACCGTTTGGTGGTGGGCAGACGGTGGTTACGCCACAAACGCAGATCGATTTACGTCAGAT‑CGGCGGTG‑CGCTg                                                                                                                                                                                                                                      <  2:233362/247‑1 (MQ=255)  gCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCGAATCGATTTACGTCAGAG‑cggcgggtacccgcc                                                                                                                                                                                                                                     >  2:147221/1‑242 (MQ=255)     gATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTCCGTCAGAG‑cggcggggtagcggcac                                                                                                                                                                                                                                   >  2:431362/1‑239 (MQ=255)     gATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTCCGTCAGAG‑cggcggcttgccgccca                                                                                                                                                                                                                                   >  1:401439/1‑239 (MQ=255)     gATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCaaa                                                                                                                                                                                                                                  >  1:204150/1‑250 (MQ=255)           cAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTggcaagcggacg                                                                                                                                                                                                                           >  1:567679/1‑240 (MQ=255)              aaTATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTc                                                                                                                                                                                                                        <  1:42897/251‑1 (MQ=255)                 atGCAGGTTAAGGTCACCCCTCAGGACGCTAAAGTAGGGATTAACGCGCGTACCGGTGCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCTTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAg                                                                                                                                                                                                                     <  2:398670/251‑1 (MQ=255)                       gTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCGCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAg                                                                                                                                                                                                               <  2:261286/251‑1 (MQ=255)                             gTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑GTGTGCAAAGAGTGCGATCTCGCGCCAGCCTCa                                                                                                                                                                                                          >  1:290973/1‑250 (MQ=255)                                                          tAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTTTCG‑TGCCAAGAGTACGGTCCCACGCCAGCCTCAAGAACGCGGTgggcgggcccaaagccccgg                                                                                                                                                                            >  1:425182/1‑236 (MQ=255)                                                          tAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTg                                                                                                                                                                             <  1:190959/250‑1 (MQ=255)                                                                gcgcACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTa                                                                                                                                                                      <  1:152851/251‑1 (MQ=255)                                                                 cgcACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCCCGCCAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGGTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTAc                                                                                                                                                                     <  2:387595/251‑1 (MQ=255)                                                                    aacggttcggttGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGcc                                                                                                                                                                  <  2:186092/249‑1 (MQ=255)                                                                          tCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATgg                                                                                                                                                             <  1:431362/250‑1 (MQ=255)                                                                          tCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATgg                                                                                                                                                             <  1:147221/250‑1 (MQ=255)                                                                           gggtgttGATTAATCGCGAAGTGACCCTCGACAGCGGCGCTGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGGTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACTTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGTCTCTACGCCGATgtat                                                                                                                                                           <  2:387006/250‑4 (MQ=255)                                                                            ggtggtGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTAc                                                                                                                                                                     <  1:509112/240‑1 (MQ=255)                                                                                   atgaATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTc                                                                                                                                                  <  1:420347/251‑1 (MQ=255)                                                                                       aTCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATa                                                                                                                                             <  1:419105/251‑1 (MQ=255)                                                                                                 gacctccGACTGCTGCTCGGTAGCGCAGGGGAATCGCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGTTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTT‑ATTTCC‑ATACTgtaagcaa                                                                                                                                   <  2:354032/245‑8 (MQ=255)                                                                                                 gACCATCGACTGCCGCGCGGTCGAGCAGGGGAATCCCTCAGTAAGAGGTAATCGTCAGGCCAATGTCATCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGAGCTATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTGCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGAGCTG‑ATGTCC‑ATACTGcaatca                                                                                                                                    <  2:425346/250‑1 (MQ=255)                                                                                                    ccTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGc                                                                                                                                <  1:142934/251‑1 (MQ=255)                                                                                                      ccGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCaa                                                                                                                              <  1:260137/250‑1 (MQ=255)                                                                                                               gcgcGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGGCAGGCCAATGTCAGCCAGCCAGATACACCGGTTGGTGGTGGACAGACTGTGGTTACACCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGGCC‑ATACTGCAAGCAATGCAAA‑GTGCggg                                                                                                                     <  2:401780/250‑1 (MQ=255)                                                                                                                      aGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGc                                                                                                             <  1:250920/251‑1 (MQ=255)                                                                                                                            ggggAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAAAGTAAGTTCCAGCGCCCACCCCCATAACGTGGGGCGCGCGCTCCACGCGCTGGGCGCTACGCCCGTGGGTCTGTATGTCC‑‑TAATTCAATCAATGAAAA‑AAGCGGGGAGTTTGCGGGCa                                                                                                        >  1:433046/1‑250 (MQ=255)                                                                                                                                  tctctcAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑TGCTGCCAAACGTACGTTCCCGCGCCAGCCTCAATAACGTGGTGCGGGCGCCCAATGCGCTGGGGGCTACACCGCTGGGTCTG‑ATGTCC‑ATAATGCAATCAATGCCAA‑GGGCGGGGTGGCTGCGGGGAAAACTgg                                                                                                 >  1:110166/1‑251 (MQ=255)                                                                                                                                             aaGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGCGGCCAGACTGTGGTTACTCCACAAACGCAGATCGATTTAAGTCAGAG‑TGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCGCAATAACGTGGTGCGCGCGCTCAATGCGCTGGTCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCtga                                                                                      <  1:380966/250‑1 (MQ=255)                                                                                                                                                   aTCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAACGCACGTTTCAGCGCCAGCCTCCATAAAGTGGGGGGCGCGCTCAATGCGCGGGGGGCTACGCCCGTGGGGCTG‑ATGTGC‑CTACTGCAATCAATTCAAAAGGACGGG‑TGTGTGCGGGCGAAACCGGAAAACATCttgtggttc                                                                                >  1:36862/1‑243 (MQ=255)                                                                                                                                                      gTCATGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCGGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGAGCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGGCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCg                                                                             <  2:422962/251‑1 (MQ=255)                                                                                                                                                      gTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCGCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGa                                                                                                <  2:286883/232‑1 (MQ=255)                                                                                                                                                      gTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑AGCTGCAAAGCGGACGTTTCCGCGCCCCCCTCCATAAAGTGGGGCGCGCGCCCAATGCGCTGGGCGCCACCCCGATGGGTCTT‑ATGTCC‑ATACTGCAATCAATGCAAA‑GGGCGGGAGGTCTGCGGGCAAAAATGGaagaacggaagagcacacac                                                                             >  1:286883/1‑233 (MQ=255)                                                                                                                                                         aGGCCAATGTCAGCCAGCCAGATACCCCGTTTGGTGGTGGACAGACTGGGGTTACCCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTG‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAatcat                                                                                          <  2:252414/235‑1 (MQ=255)                                                                                                                                                         aGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCCAAACGTACGTTCCAGCGCCAGCCCCACAAACGTGGGGCGCGCGCTCAACGCGCTGGGCGCTACCCCGCTGGATCTG‑ATGTCCCATA‑TTCAATCAAATCAAA‑GTGCGGGGTGTCTTCTGGCGAAAATGGAAaacatagatcggaaaagcaac                                                                          >  1:252414/1‑235 (MQ=255)                                                                                                                                                                 gTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACGGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAacta                                                                  <  1:33359/251‑1 (MQ=255)                                                                                                                                                                                gcacCGGTTGGTGGGTGACAGACTGTGGTTACCCCACAAACGCAGATCGATGTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑AGACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCgg                                                     <  2:425182/248‑1 (MQ=255)                                                                                                                                                                                           gtgatgACAGACTGTGGTTACTCCACAAACGCAGATCTAGTTACGTCAGAG‑CGGCGGGT‑CGCTGCAAAGCGTTCGTTCCAGCGTCAGCCTCAATAACGTGGGGCGCGCGCCCAATGCGCTGGGCGCTACGCTGATGGATCTG‑ATGGCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCa                                        <  2:290973/246‑1 (MQ=255)                                                                                                                                                                                            ggtggACAGACTGTGGTGACTCCACAAACGCAGATCGAGTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCa                                        <  2:567679/250‑1 (MQ=255)                                                                                                                                                                                                  cAGACTGTGGGTTCGCCACAACCGCAGATCGACTTACGTC‑GAGGCGGTGGTT‑CGCTTCCACGCGTACGTTCCAGCGCCTGCCTCAATAACGTGGTGCGCGCGCACAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAACCAAAGTAAA‑GTGCGGGATGTCTGCTGGCAAAACTGGAAATCATCTGATGATCAGCGAAAGTAAACTACTCGCAAGTGCGGCCTGGGATGCTCAAtctctc                                 <  2:166749/251‑1 (MQ=255)                                                                                                                                                                                                  aagattgttgTTACTCCACAAACGAAGAGCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGGACGTTCCAGCGCCAGCCGCAATAACGTGGTGCGCGCGCTCAATGCGCGGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTc                                 <  2:542757/242‑1 (MQ=255)                                                                                                                                                                                                        gtgGTTACTCCACAAACGCAGATCGATTTACGCCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGa                            <  1:283884/250‑1 (MQ=255)                                                                                                                                                                                                                 ccACAATCGCAGATCGATTTACGCCAGAG‑CGGCGTTG‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAGCGTGGTGGGCGCTCTAAATGCGCTGGGCGCTACGCTGATGGATCTG‑ATGTCC‑ATACTGCAATCATTGCAAA‑GTGCGGGATGTCTGCGGGCAAAACGGGAAATCATCTGATGATCAGCGACATCAAACTGCTGTCAAGTGCGGCCTGTGTGGCGCAGTCACTGAACGAACGAAAgg                    <  2:518244/249‑1 (MQ=255)                                                                                                                                                                                                                     aaaCGCAGAGCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAg              <  1:63714/251‑1 (MQ=255)                                                                                                                                                                                                                            gATCGATTTACGTCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACTTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGcaa                                                                                                                                       <  1:241460/124‑1 (MQ=255)                                                                                                                                                                                                                                  tttACGTCAGTG‑CGGAGGTG‑CGCGGCAAAGCGTACGTTCCAGCGCCAGGCTCAACAACGTGGTGCGGGCGCTCAATGCGCTGGGCGCGTCGCCGATGGATCTG‑ATGTCC‑ATACTGCAAACAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACGCAACGAACTAAAGGCGAAAGCGGGCGAAGATc  <  2:204150/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         TGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGCCAGAG‑CGGCGGTT‑CGCTGCAAAGCGTACGTTCCAGCGCCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTG‑ATGTCC‑ATACTGCAATCAATGCAAA‑GTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAGCGGGCGAAGATC  >  NC_000913/1136981‑1137455

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