Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,014,217	C→A	22.0%	A330E (GCA→GAA)	fliF →	flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,014,217	0	C	A	22.0%	84.5 / 3.0	42	A330E (GCA→GAA)	fliF	flagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base C (30/2);  new base A (2/8);  total (32/10)
Fisher's exact test for biased strand distribution p-value = 1.54e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCT  >  NC_000913/2013992‑2014448                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         gTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATaa                                                                                                                                                                                                                              >  1:637183/1‑238 (MQ=255)                aCGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCaa                                                                                                                                                                                                                                        >  2:672949/1‑213 (MQ=255)                aCGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCaa                                                                                                                                                                                                                                        <  1:672949/213‑1 (MQ=255)                aCGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCg                                                                                                                                                                                                              >  1:640511/1‑239 (MQ=255)                  gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCAAGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGa                                                                                                                                                                                                             >  1:1032231/1‑238 (MQ=255)                  gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGGTGTCGACTCAACCGGACCCTCCGAAGAACCCCCCACTACGCACGCCTCCGGCAACCCAAAACAACCCCCAGACGCAGGcaaa                                                                                                                                                                                                            >  2:81243/1‑236 (MQ=255)                                 cAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAACCGGTGCCGGTTATCCGGGC‑GGGGCACCGGGGGCGTTGTCGAATCAAACGGCACAGCCGAAAAACGCGCCAATCGCTACGCCCCCGCAACAACAAAACAAACCCCCGCACCAGCCGAACAccaccaaccaaa                                                                                                                                                                                              >  2:403543/1‑232 (MQ=255)                                       ggcTTCCCCAGTAAAAAAAAAAAGAAGAAACAGTTGCGCCCAAAGGGTGATGATTCCTGTGGCGCGT‑‑TTGGTTCACGCCAGCTTAATGAGAGCGAGAAAAGTGGTCCGGGTTA‑CCGGGCGGGCGTACCGGGGGCGTTGTCGAACCAACCGGCACCTGCAAATAACGCGCCAATAAGCACGCCCCCGGAAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGcc                                                                                                                                                                                       <  2:1010093/237‑1 (MQ=255)                                         cgtCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAACCAGCACGCCCCCGGAAAACCAAAATAACCGCCAGCAGCAGGCGGCGCCCACCAGCAAAAACGGGCcgc                                                                                                                                                                                     >  2:987759/3‑239 (MQ=255)                                            cGCCAATAAAGAAAAAAGGGAAGAACTGTGTCGCCATAGGGGTG‑TGAA‑CTCATGTGGCGCCTTTCGTCCCCGCCTGTTAAGTGGGGGAGAGAAACGCGTTCGCGGTTTCCGGGGG‑GGCGCACCGGGGGGGTTTCCAAACAAACCGGCCCCTGCAAAAACCCCGCAAATCACCCCCCCCCCGAAAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTa                                                                                                                                                                                  <  2:376305/239‑1 (MQ=255)                                              ccAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCCGCAAAACTggctgcacctta                                                                                                                                                                                >  1:397534/1‑229 (MQ=255)                                                                       gTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGGAAcg                                                                                                                                                       >  1:808968/1‑238 (MQ=255)                                                                        tATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCCGGCGAGCCCCACCAGCAATAATGGGCCGCGTACCACACAGCGGAGAGAAACCAGCaaaca                                                                                                                                                      >  2:1075559/1‑236 (MQ=255)                                                                             cccTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGAAATAGTGGGCGGAGTCGCACACAGCGGAATGAAACCAGTAACTaagaag                                                                                                                                                 >  1:662982/1‑239 (MQ=255)                                                                             cccTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCACCAATAGTGGGCCGCGTAGCACACACCGGAATGAAACCAGTAACACCGaa                                                                                                                                                  >  1:401092/1‑238 (MQ=255)                                                                                  aCGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCCGCAATAGTGGGCCGCGTAGCAACCAGCGGAATGAAACCAGAAACCACGaaagcaa                                                                                                                                             >  2:760772/1‑233 (MQ=255)                                                                                  aCGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGCAACTACGaaagcga                                                                                                                                             >  1:377910/1‑233 (MQ=255)                                                                                        atgaATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGATCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCAcc                                                                                                                                     >  1:574204/1‑239 (MQ=255)                                                                                                                   aCGCCAGCTTAATGTGAGCGAGCAAAGCGGATCCAGTTATCCGGGC‑GGCGCACCGGGGGCGTTGTCGAATCAACCAGCACCTGCGAATAACGCGCCAATCAGCACGCCTACGGCACATCAAAATAACCGCCAGCAGCAGGGGAGCCCCACCCGCAAAAGGGGGCCGCCTAACACCCACCGGGATGGAAACACAAGCTCCGAATACGAAGCGCCCCA‑TCGTCATACCAAAATGAACGTgg                                                                                                            >  2:155740/1‑239 (MQ=255)                                                                                                                                   agcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACATGGGCGATGTGCAACGTCt                                                                                            >  1:211836/1‑239 (MQ=255)                                                                                                                                   agcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCAAACAGCGGAAGGAAACCAGTCACAACGAAGTCG‑AGCGCACCATTCGCCATACCAAAAAGAACGGGGGCGATGTGAAACg                                                                                               >  2:111464/1‑236 (MQ=255)                                                                                                                                    gcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGCAATGAAACCAGTAACTACGAAGTCG‑ATCGAACCATTCGTCATACCAAAATGAACGTGGACGATGTGAAActact                                                                                            >  1:133093/1‑234 (MQ=255)                                                                                                                                       agcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAAGGAACGTGGGCGAGTGGCAACGCCttcca                                                                                        >  1:645204/1‑235 (MQ=255)                                                                                                                                       agcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCAGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAATCG‑ATCGCACCATTCGTCATACCAAAATGAAAGAGGGCAAAGGGCAACggcgggca                                                                                        >  2:283728/1‑232 (MQ=255)                                                                                                                                                   cgggTTACCGGGGC‑GGCGTCCGGGGGGTGTTGTCAAATCAACCGGCACCTGCGAATAACGCGCAAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCgtg                                                                            <  1:760772/237‑1 (MQ=255)                                                                                                                                                              ggc‑gccGAACCGGGGGCTTGGTGGAATCACACAGACTCTGTAAATGCCCCGCCAACCTGCACGCCTCGGAAAAATAAAAAAAACCGCAAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACaa                                                                 <  2:574204/239‑1 (MQ=255)                                                                                                                                                                ggggcttgccGGGGGCGTTGACAAACCAACCGGCCCCGGAAAAAAACGCCCCAACCGGCCCCCCCCGGGAAAATAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACaa                                                                 <  1:545564/230‑1 (MQ=255)                                                                                                                                                                      tACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGCCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGCCTGTCAGTCGCGGCCGGGGTGAATTACAAAACCtt                                                           >  1:1039147/1‑238 (MQ=255)                                                                                                                                                                          gggggCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGGCATACCAAAATGAACGTGGGCGACTTGCAACGTCGGTCAGTGCCGGACGTGGGGAATAAAAAAAACTAGCCAg                                                      >  2:111710/1‑239 (MQ=255)                                                                                                                                                                               gttgttCGAACCACCCGTCACCGGCAAATAACGCGCAAATCAGCCCGCCGCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGGGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTTTGTCAGTCGCGGTTGTGGTGAATTACAAAACCTTGCCAGATgg                                                  <  1:155740/233‑1 (MQ=255)                                                                                                                                                                                gtgtgcGAATCAGCCGGCACAGGCAAAAAACGCCCCAACAAGCACGCCCCCGGAAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGc                                                 <  1:403543/233‑1 (MQ=255)                                                                                                                                                                                              cGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGGGGTGAATTCCAAACACTTGCCAGATGGCAAACCGTTGcccccg                                  >  1:889618/1‑235 (MQ=255)                                                                                                                                                                                                  aCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCaaa                                                                                                                      >  1:759019/1‑151 (MQ=255)                                                                                                                                                                                                      ggaaaaaaCGCGCCAACCAGCCCGCCACCGGAAAAAAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGa                          <  1:1010065/233‑1 (MQ=255)                                                                                                                                                                                                       cGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGc                                       >  1:728261/1‑225 (MQ=255)                                                                                                                                                                                                        aaaTACCGCCCAAAACACCACGCCTCCGGAAAAAGAAAAAACCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGaa                         <  1:1031351/237‑1 (MQ=255)                                                                                                                                                                                                              cgcgCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGAACACAGCGGAATGAAACCAGTAACTACGACGCCG‑CTCGCACCATTCGTCATACCAAAATGAACGTCAGCGATGTGCAACGTCGGCCACTCGCCGTCGTGGTGACTTTCAAAACCTAGCCAATCGGAAAACCGGTTCCCAGCAGCAAcacgaaaaaga                  >  2:898189/1‑229 (MQ=255)                                                                                                                                                                                                                 gcCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAAATGGCAAACCGTTGCCTCCCAGCAACGAACAGATGAAGc               >  1:505725/1‑239 (MQ=255)                                                                                                                                                                                                                 gcCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGCCGGGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGc               >  1:1090022/1‑239 (MQ=255)                                                                                                                                                                                                                             gCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGACTGTCAGTCGCGGTCGTGGTGAATTACaaaa                                                               >  1:549501/1‑179 (MQ=255)                                                                                                                                                                                                                               cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGCAGCAAATTGAAGATCt  >  1:161195/1‑238 (MQ=255)                                                                                                                                                                                                                               cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCCCGGTCGTGGGGAATTACAAAACCTTGCCAGATTGCAAACCGTGGCCTCTCAGCAACGACAAGATCAAGCAAAtttacaatct  >  1:1013067/1‑230 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         GTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCT  >  NC_000913/2013992‑2014448

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