Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,014,217	C→A	21.8%	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	21.8%	91.1 / 3.5	47	A330E (GCA→GAA)	fliF	flagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base C (34/2);  new base A (2/9);  total (36/11)
Fisher's exact test for biased strand distribution p-value = 2.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCT  >  NC_000913/2013992‑2014448                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          gTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATaa                                                                                                                                                                                                                              >  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‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCg                                                                                                                                                                                                              >  1:640511/1‑239 (MQ=255)                  gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCAAGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGa                                                                                                                                                                                                             >  1:1032231/1‑238 (MQ=255)                  gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGGTGTCGACTCAACCGGACCCTCCGAAGAACCCCCCACTACGCACGCCTCCGGCAACC‑CAAAACAACCCCCAGACGCAGGcaaa                                                                                                                                                                                                            >  2:81243/1‑236 (MQ=255)                       ggTTACGGCGCAGCTGGACTTCGCCAGTCAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGTCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGTGCAAAGCGGTTCCGGTTATCCGGGC‑GCCGTACCGGGGGCGGTGTCGAATCAACCGGAACCTGGGGAAAACGCCCCCATCGGCACGCCCCCGGCAAATACAAAAAAAccccaaaaaaaaacgacaaaaa                                                                                                                                                                                                        >  2:1047227/1‑219 (MQ=255)                                 cAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAACCGGTGCCGGTTATCCGGGC‑GGGGCACCGGGGGCGTTGTCGAATCAAACGGCACAGCCGAAAAACGCGCCAATCGCTACGCCCCCGCAACAA‑CAAAACAAACCCCCGCACCAGCCGAACAccaccaaccaaa                                                                                                                                                                                              >  2:403543/1‑232 (MQ=255)                                       ggcTTCCCCAGTAAAAAAAAAAAGAAGAAACAGTTGCGCCCAAAGGGTGATGATTCCTGTGGCGCGT‑‑TTGGTTCACGCCAGCTTAATGAGAGCGAGAAAAGTGGTCCGGGTTA‑CCGGGCGGGCGTACCGGGGGCGTTGTCGAACCAACCGGCACCTGCAAATAACGCGCCAATAAGCACGCCCCCGGAAAAT‑AAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGcc                                                                                                                                                                                       <  2:1010093/237‑1 (MQ=255)                                         cgtCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAACCAGCACGCCCCCGGAAAAC‑CAAAATAACCGCCAGCAGCAGGCGGCGCCCACCAGCAAAAACGGGCcgc                                                                                                                                                                                     >  2:987759/3‑239 (MQ=255)                                            cGCCAATAAAGAAAAAAGGGAAGAACTGTGTCGCCATAGGGGTG‑TGAA‑CTCATGTGGCGCCTTTCGTCCCCGCCTGTTAAGTGGGGGAGAGAAACGCGTTCGCGGTTTCCGGGGG‑GGCGCACCGGGGGGGTTTCCAAACAAACCGGCCCCTGCAAAAACCCCGCAAATCACCCCCCCCCCGAAAAAC‑AAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTa                                                                                                                                                                                  <  2:376305/239‑1 (MQ=255)                                              ccAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCCGCAAAACTggctgcacctta                                                                                                                                                                                >  1:397534/1‑229 (MQ=255)                                                            aCGGACGCACAGTATCGATCTAACGGTGATGAATCTAATGCGGCGC‑‑CTCGTTCACGCCAGCTTAATGAGAGCGACCAAAGCGGATCCGGTTATCCGGGC‑GGCGCACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAAGAACGCGCCCATAAGCACCCCTCCGGCCAAG‑CAAAATACCCGCCATCAGCAGGCAAGCACCCCCCGCAATAGAGGGcagaggcgcacccacaggaacta                                                                                                                                                                  >  2:706952/1‑217 (MQ=255)                                                            aCGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGcctgtctcttataaacatatccg                                                                                                                                                                  >  1:288341/1‑218 (MQ=255)                                                                       gTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGGAAcg                                                                                                                                                       >  1:808968/1‑238 (MQ=255)                                                                        tATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCCGGCGAGCCCCACCAGCAATAATGGGCCGCGTACCACACAGCGGAGAGAAACCAGCaaaca                                                                                                                                                      >  2:1075559/1‑236 (MQ=255)                                                                             cccTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGAAATAGTGGGCGGAGTCGCACACAGCGGAATGAAACCAGTAACTaagaag                                                                                                                                                 >  1:662982/1‑239 (MQ=255)                                                                             cccTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCACCAATAGTGGGCCGCGTAGCACACACCGGAATGAAACCAGTAACACCGaa                                                                                                                                                  >  1:401092/1‑238 (MQ=255)                                                                                  aCGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCCGCAATAGTGGGCCGCGTAGCAACCAGCGGAATGAAACCAGAAACCACGaaagcaa                                                                                                                                             >  2:760772/1‑233 (MQ=255)                                                                                  aCGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGCAACTACGaaagcga                                                                                                                                             >  1:377910/1‑233 (MQ=255)                                                                                        atgaATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGATCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCAcc                                                                                                                                     >  1:574204/1‑239 (MQ=255)                                                                                                                   aCGCCAGCTTAATGTGAGCGAGCAAAGCGGATCCAGTTATCCGGGC‑GGCGCACCGGGGGCGTTGTCGAATCAACCAGCACCTGCGAATAACGCGCCAATCAGCACGCCTACGGCACAT‑CAAAATAACCGCCAGCAGCAGGGGAGCCCCACCCGCAAAAGGGGGCCGCCTAACACCCACCGGGATGGAAACACAAGCTCCGAATACGAAGCGCCCCA‑TCGTCATACCAAAATGAACGTgg                                                                                                            >  2:155740/1‑239 (MQ=255)                                                                                                                                   agcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACATGGGCGATGTGCAACGTCt                                                                                            >  1:211836/1‑239 (MQ=255)                                                                                                                                   agcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCAAACAGCGGAAGGAAACCAGTCACAACGAAGTCG‑AGCGCACCATTCGCCATACCAAAAAGAACGGGGGCGATGTGAAACg                                                                                               >  2:111464/1‑236 (MQ=255)                                                                                                                                    gcgagcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGCAATGAAACCAGTAACTACGAAGTCG‑ATCGAACCATTCGTCATACCAAAATGAACGTGGACGATGTGAAActact                                                                                            >  1:133093/1‑234 (MQ=255)                                                                                                                                       agcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAAGGAACGTGGGCGAGTGGCAACGCCttcca                                                                                        >  1:645204/1‑235 (MQ=255)                                                                                                                                       agcAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCAGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAATCG‑ATCGCACCATTCGTCATACCAAAATGAAAGAGGGCAAAGGGCAACggcgggca                                                                                        >  2:283728/1‑232 (MQ=255)                                                                                                                                                   cgggTTACCGGGGC‑GGCGTCCGGGGGGTGTTGTCAAATCAACCGGCACCTGCGAATAACGCGCAAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCgtg                                                                            <  1:760772/237‑1 (MQ=255)                                                                                                                                                              ggc‑gccGAACCGGGGGCTTGGTGGAATCACACAGACTCTGTAAATGCCCCGCCAACCTGCACGCCTCGGAAAAAT‑AAAAAAAACCGCAAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACaa                                                                 <  2:574204/239‑1 (MQ=255)                                                                                                                                                                ggggcttgccGGGGGCGTTGACAAACCAACCGGCCCCGGAAAAAAACGCCCCAACCGGCCCCCCCCGGGAAAAT‑AAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACaa                                                                 <  1:545564/230‑1 (MQ=255)                                                                                                                                                                      tACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGCCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGCCTGTCAGTCGCGGCCGGGGTGAATTACAAAACCtt                                                           >  1:1039147/1‑238 (MQ=255)                                                                                                                                                                          gggggCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGGCATACCAAAATGAACGTGGGCGACTTGCAACGTCGGTCAGTGCCGGACGTGGGGAATAAAAAAAACTAGCCAg                                                      >  2:111710/1‑239 (MQ=255)                                                                                                                                                                               gttgttCGAACCACCCGTCACCGGCAAATAACGCGCAAATCAGCCCGCCGCCGGAAAAT‑CAAAATAACCGCCAGCAGCAGGCGGGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTTTGTCAGTCGCGGTTGTGGTGAATTACAAAACCTTGCCAGATgg                                                  <  1:155740/233‑1 (MQ=255)                                                                                                                                                                                gtgtgcGAATCAGCCGGCACAGGCAAAAAACGCCCCAACAAGCACGCCCCCGGAAAAC‑AAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGc                                                 <  1:403543/233‑1 (MQ=255)                                                                                                                                                                                 cggcccacccaacccaccccTGCAAACAACGCGCAAATCACCCCGCCTCGGAAAAAA‑CAAAATACCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCaa                                               <  2:889618/221‑1 (MQ=255)                                                                                                                                                                                              cGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACCCAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGGCCGGGGGAATTAACAAACCTTGCCaacaggaaacccgtggccccca                                  >  2:1031351/1‑218 (MQ=255)                                                                                                                                                                                              cGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGGGGTGAATTCCAAACACTTGCCAGATGGCAAACCGTTGcccccg                                  >  1:889618/1‑235 (MQ=255)                                                                                                                                                                                                  aCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCaaa                                                                                                                      >  1:759019/1‑151 (MQ=255)                                                                                                                                                                                                      ggaaaaaaCGCGCCAACCAGCCCGCCACCGGAAAAA‑AAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGa                          <  1:1010065/233‑1 (MQ=255)                                                                                                                                                                                                       cGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGc                                       >  1:728261/1‑225 (MQ=255)                                                                                                                                                                                                        aaaTACCGCCCAAAACACCACGCCTCCGGAAAAA‑GAAAAAACCCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGaa                         <  1:1031351/237‑1 (MQ=255)                                                                                                                                                                                                              cgcgCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGAACACAGCGGAATGAAACCAGTAACTACGACGCCG‑CTCGCACCATTCGTCATACCAAAATGAACGTCAGCGATGTGCAACGTCGGCCACTCGCCGTCGTGGTGACTTTCAAAACCTAGCCAATCGGAAAACCGGTTCCCAGCAGCAAcacgaaaaaga                  >  2:898189/1‑229 (MQ=255)                                                                                                                                                                                                                 gcCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAAATGGCAAACCGTTGCCTCCCAGCAACGAACAGATGAAGc               >  1:505725/1‑239 (MQ=255)                                                                                                                                                                                                                 gcCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGCCGGGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGc               >  1:1090022/1‑239 (MQ=255)                                                                                                                                                                                                                             gCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGACTGTCAGTCGCGGTCGTGGTGAATTACaaaa                                                               >  1:549501/1‑179 (MQ=255)                                                                                                                                                                                                                               cTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGCAGCAAATTGAAGATCt  >  1:161195/1‑238 (MQ=255)                                                                                                                                                                                                                               cTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑ATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCCCGGTCGTGGGGAATTACAAAACCTTGCCAGATTGCAAACCGTGGCCTCTCAGCAACGACAAGATCAAGCAAAtttacaatct  >  1:1013067/1‑230 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          GTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGC‑‑TTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGC‑GGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAT‑CAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCG‑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 < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑