Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,368,876	A→G	26.2%	E71G (GAG→GGG)	fxsA →	suppressor of F exclusion of phage T7

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,368,876	0	A	G	26.2%	108.1 / 8.5	66	E71G (GAG→GGG)	fxsA	suppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base A (12/34);  new base G (16/3);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 2.18e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  NC_000913/4368644‑4369113                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         ctgctATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGgatagt                                                                                                                                                                                                                                                    >  1:531673/1‑236 (MQ=255)              agaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGtt                                                                                                                                                                                                                                        <  2:269022/236‑1 (MQ=255)               gaAACCTTTGCGCTGGTTTCCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCg                                                                                                                                                                                                                                     <  1:740438/238‑1 (MQ=255)               gaAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTATCCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGC‑GGAGATG‑GTT‑TAAAAGGTTTctc                                                                                                                                                                                                                                    >  1:70285/1‑237 (MQ=255)                                ttccTTTTTTTGCCTTTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTtg                                                                                                                                                                                                                <  2:660815/237‑1 (MQ=255)                                  ccTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTtgct                                                                                                                                                                                                              <  1:500200/239‑1 (MQ=255)                                  ccTGTTATTGCCATTTTCCTTTATGGCTATATGGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTtgct                                                                                                                                                                                                              <  1:398266/239‑1 (MQ=255)                                        aTTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGGTG‑ATT‑AAAAATGTTTCGCC‑GA‑TCATTGCT‑GGTTTGCTGCtttt                                                                                                                                                                                                        >  1:1055254/1‑239 (MQ=255)                                                    cTTTAGTTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCg                                                                                                                                                                                                                                     <  1:233031/201‑1 (MQ=255)                                                            cTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑aaaa                                                                                                                                                                                                                                             <  1:835849/185‑1 (MQ=255)                                                              atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGGTG‑ATT‑TAAAGGGGTTTCCT‑GG‑TCATGGCT‑GGGTTGCTGCCTTTATTAC‑CGGGCTTTTTTACC‑Gaa                                                                                                                                                                                >  1:721054/1‑238 (MQ=255)                                                                  tGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCAGCGGC‑GGGGGTGCGTA‑AAAAGATTTTCG‑T‑GT‑TCTTTTCG‑GGTTTGGTGCGTTTTTTAC‑AGGGGGTTTTTTCC‑GACTTc                                                                                                                                                                             >  2:575333/1‑238 (MQ=255)                                                                     gaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCACGGA‑GGGGGTG‑GTT‑AAAAAGGTTTTCGT‑GA‑TAATTGTT‑GGGTTTTTGGTTTTTTTTA‑CGGGGTTTTTTTAC‑GGCTTCCCCg                                                                                                                                                                         >  2:1106456/1‑239 (MQ=255)                                                                       tttCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCCGCGGC‑GGGGATG‑ATT‑AAAAAGGTTTTCCT‑TG‑TGATTGGTGGGTTTGTGGCTTTTTTTAC‑CGGGTTTTTTTTAC‑GATTTCCCCgg                                                                                                                                                                        >  2:135992/1‑239 (MQ=255)                                                                          tactCTTTTTTTTGGTTGCCCATTTTTTTGGGGTTTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGCCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGt‑ctt                                                                                                                                                                   <  2:187831/236‑1 (MQ=255)                                                                          cAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGt‑ctt                                                                                                                                                                   <  1:698041/239‑1 (MQ=255)                                                                                ttATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAACCACCGGC‑GGGGGTGTGTT‑AAAAGGGTTT‑GCTGGA‑CAATTGTG‑GGTTTCTTGCTTTTTTTCC‑AGGCTTTTTTTTCC‑CACTTCCCCgcttttttctt                                                                                                                                                               >  2:10115/1‑229 (MQ=255)                                                                                        gTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTTCC‑GACTTCCTCGGG‑CTTCTACTtttttttc                                                                                                                                                      >  1:469523/1‑236 (MQ=255)                                                                                          ttgcctgtttttGGGGGGTTTGGTGACCCTTGTGCGGGTTATATTCCCGTCAGTTATCGGTATGGACCTGGTACGTAACCAGGGTTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTgccgc                                                                                                                                                   <  2:469523/229‑1 (MQ=255)                                                                                          ttcccTTGTTTTGGGGGTTTTGCTGACCCTCGTGCTGGTTATTTTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGTTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTgccgc                                                                                                                                                   <  2:730038/237‑1 (MQ=255)                                                                                           gCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGG‑GGGGAAGTATT‑AAAAG‑GTTTTCCT‑GATTCATGTCC‑GTTTTTGCGCTTTTTTTA‑‑CGGGCTTTTTTTCCTCTCTTCCTCGGT‑‑TTTTTCTTTTAttttcccc                                                                                                                                                  >  2:318201/1‑233 (MQ=255)                                                                                           gCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGG‑GGAGATG‑ATT‑TAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGGTTTCTGCTTTTTTTAC‑CGGGGTTTTTTAAC‑GACTTCCTCGGG‑CTTCTACTTTTATTgccgcc                                                                                                                                                  >  1:927138/1‑239 (MQ=255)                                                                                           gCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGGGTTTCGCT‑GA‑TCATTGCT‑GGTTTTCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGG‑CTTCTACTTTTATTtccgcc                                                                                                                                                  >  1:730038/1‑239 (MQ=255)                                                                                                  ggtagattgtgataagagaacgggtggtgGTTATTTTCAAGTCAGTTGTCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAACCCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                           <  2:583420/216‑1 (MQ=255)                                                                                                  agacaggggTATTGCTGACCCTCGTGCTGGTTAGATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                           <  1:708306/234‑1 (MQ=255)                                                                                                   ttttgggtttttgttgACCCTGGTGCTGGTTATATTCCAGTCATTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                           <  2:761232/228‑1 (MQ=255)                                                                                                   aTTGGGGGTATTGCTGACCCTCGTGCTGGGTATATGCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAAGTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCaa                                                                                                                                           <  1:263609/238‑1 (MQ=255)                                                                                                      gggggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAAACCAGCGGCGGGGGATG‑GTTAAAAAGTG‑TTCCTT‑GT‑GCTGTTCG‑GGTTTGTGGCTTTT‑TTTC‑CGGGCGTTTTTTCC‑CCCTTTCCCGGTGTTTTTTCTTTTTTTGCCCCCGGGGCaaaaa                                                                                                                                        >  2:655814/1‑238 (MQ=255)                                                                                                       ggggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGATG‑ATT‑AAAAGTGGTTCGCT‑GG‑TCATTGCT‑GGGTTTCTGCTTTTTTTAC‑CGGGGTTTTTTTCC‑GACTTCCTCGGT‑CTTTTTCTTTTATTTCCGCCGGTGCAAAAACa                                                                                                                                      >  1:262810/1‑239 (MQ=255)                                                                                                         ggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑aaaa                                                                                                                                                                                                                                             <  1:717039/140‑1 (MQ=255)                                                                                                           tATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGGGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTTCC‑GGATTCTCCGGG‑CTTCTCATTTTTTTGCCGCCGGGGCAAAAGAATCTg                                                                                                                                  >  1:487774/1‑239 (MQ=255)                                                                                                            aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGACAACCAAGCGGC‑GGGGTAG‑GAT‑AAAACAGGTTCGCG‑GA‑TCCTTATC‑GGGTTGCTGATTTTATTCA‑CGGGCTTTTTTTTC‑CACTTCCTCTGT‑CTTTTCTTTTTTTTTCCTCCGCTGTTAAACCAtatca                                                                                                                                 >  2:194731/1‑235 (MQ=255)                                                                                                              tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GG‑TCATTGCT‑GGGTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCCGTGCAAAAGCATCTGAcc                                                                                                                               >  1:178237/1‑238 (MQ=255)                                                                                                              tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACa                                                                                                                               <  1:575333/239‑1 (MQ=255)                                                                                                                   aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGATG‑ATT‑AAAAATGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGGCTTTTTACC‑GACTTCCCCGGG‑CTCTTACTTTTATTGCCGCCGGTGCAAAAGCATCTTACAGTGaa                                                                                                                          >  1:1068332/1‑239 (MQ=255)                                                                                                                   aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGC‑GGGGGTG‑ATT‑TAAAGGGTTTTCCT‑GG‑TCATTGGT‑GGGTTGTTGCTTTTATTTC‑CGGGGTTTTTTTCC‑GCCTTTCTCGGT‑CTTTTTCTTTTTTTTCCTCCgggggcaaaaaattttcaccgc                                                                                                                            >  2:708306/1‑217 (MQ=255)                                                                                                                     ccTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGGTTGCTGCTTTTATTCA‑CCGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACCTTTTTTTCCGCCGGTGCCAAAGCATCTGACAGTTAAgg                                                                                                                        >  1:93834/1‑238 (MQ=255)                                                                                                                                   tatTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGATG‑ATT‑AAAAGTGTTTCGCT‑GT‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGa                                                                                                                                 >  1:583420/1‑216 (MQ=255)                                                                                                                                     atCACTTTAGGTATCGGGATGTCCCTGGTACCTAACCAGGGCTTTAAGAATTTCGTGCGGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                        <  1:194731/238‑1 (MQ=255)                                                                                                                                         cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGC‑GGAGATG‑GTT‑AAAAATGTTTCGCT‑GA‑TCATTGCT‑GGGTTGAGTCTTTTATTACACGGGTGTTTTTTCC‑ACCTCCCTC‑GT‑CTTTTTCTTTTTTTTGCGCCGGGGCCAAAACATCTCGCACTGATGTTTATGCCGCATTTggcgttt                                                                                                    >  2:67266/1‑233 (MQ=255)                                                                                                                                          gTCAGTTTTCGTTTTGTGATTGGTACGTAACCAGGGGTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGG‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCtttt                                                                                                   <  2:531673/239‑1 (MQ=255)                                                                                                                                              gTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAgtg                                                                                                                                                                                                                                          <  1:157764/106‑1 (MQ=255)                                                                                                                                                        tGTCACTGGTAAGTAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑TGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggc                                                                                     <  1:1054562/239‑1 (MQ=255)                                                                                                                                                              tGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCCCCTTTGCGCTTTTTTTGCCGTCCTGGC‑‑GGCGGt                                                                             >  1:601593/1‑239 (MQ=255)                                                                                                                                                                gTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGGTG‑CTT‑AAAAAGGTTTTGCT‑GG‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTTCC‑GACCTACTCCGT‑CTTCTACTTTTTTTTCCGCCGGTGCAAAAACATCTGCCAGTGAAGTTGATGCCCCATTTGCGCTTTTTTCGCaaaccttgggggggggt                                                                             >  2:477814/1‑223 (MQ=255)                                                                                                                                                                           gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGGGATG‑ATT‑TAAAGAGTTTCGCC‑GA‑TCATTTCT‑GGTTTTGCGCGTTTTTTAA‑CGGGCTTTTTTACC‑CACTTTCCCGGG‑GTTCTTCTTTTTTTTTCCCCGGTGGCAAAAAATCTTTAAGTGGAGGTTATGTCCCATTTTCCCTTTTTTCCCATGTCCGGCGGGGGGGTTTTA‑GCGCggg                                                                 >  2:163194/1‑239 (MQ=255)                                                                                                                                                                           gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGG‑CTTCTACTTTTATTGCCGCCGGTTCAAAAGCATCTTACATTGAAGGTGATGCCGCCTTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTATGCGCCgg                                                                 >  1:842208/1‑238 (MQ=255)                                                                                                                                                                                ttAAGAATTTCTTTGGGATGATGAAAAAAATGGGGGCGGGTGAAACCCCAGCGGG‑GGGGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCCCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGt                                                          <  2:927138/239‑1 (MQ=255)                                                                                                                                                                                ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGt                                                          <  1:205447/239‑1 (MQ=255)                                                                                                                                                                                ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGt                                                          <  1:690341/239‑1 (MQ=255)                                                                                                                                                                                ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGGGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCCGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCCTTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTT‑TCGCCCGGTCCCGt                                                          >  1:197694/1‑239 (MQ=255)                                                                                                                                                                                ttAAGAATTTCGGTCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGt                                                          <  1:163194/239‑1 (MQ=255)                                                                                                                                                                                agaccgattgcggtcggctgcgacaAACAATGTCGGCGGGTAAAAACCCAGCGTC‑GGAGATG‑ATT‑AAAAGGGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGt                                                          <  1:582481/216‑1 (MQ=255)                                                                                                                                                                                       tttCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTa                                                   <  1:655814/239‑1 (MQ=255)                                                                                                                                                                                       tttCGTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCGGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTa                                                   <  1:930978/239‑1 (MQ=255)                                                                                                                                                                                           gTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATAc                                               <  1:135992/239‑1 (MQ=255)                                                                                                                                                                                               tGAGGCATCAAAAAATGGCCGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACttt                                            <  1:51598/238‑1 (MQ=255)                                                                                                                                                                                                     ataaaaaatTGGCGGGGGGGGAAAACCCAGCGGG‑GGGATTG‑ATT‑AAAGTTGTTTCTGT‑GA‑TCTTTGCT‑GGTTTGCTGTTTTTATTCC‑CGGGCTTTTTTACC‑GACTTCCCCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                     <  2:178237/236‑1 (MQ=255)                                                                                                                                                                                                        aaaaaaGGGGGGGGGGTTAAAACCCGGCGGG‑GGGGATG‑ATA‑AAAAGGGTTGGGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTGTAACC‑GAGTTCCTCGGT‑CTTTTACTTTGATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTGTGATGGt                                     <  2:197694/236‑1 (MQ=255)                                                                                                                                                                                                        aaaaaaGGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGGT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                     <  1:1106456/236‑1 (MQ=255)                                                                                                                                                                                                           aaaTTGCCGCGGGTTAAAACCCTGGGGG‑GGAAATG‑TTT‑AAAAGTTTTTCGCT‑TA‑TCATTGCT‑GGTTTTCTGCTTTTATTTC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                               <  2:487774/239‑1 (MQ=255)                                                                                                                                                                                                                gcggcgGGTTAAAACCCAGCGGG‑TGTGATT‑ATG‑AAAAGTGTTCCGCT‑GA‑TCATTTCT‑GGTTTGCTGCTTTTAGTAC‑CGTGCTTTTTTGCC‑GTCTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGGGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTTATGGTGAGTACCAGCg                          <  1:551453/239‑1 (MQ=255)                                                                                                                                                                                                                   gcggGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                        <  1:787237/238‑1 (MQ=255)                                                                                                                                                                                                                   gcggGTGAAAACCCAGCGGC‑GGAGATG‑ATA‑AAAAGTGTTACCGT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTCC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                        <  1:67266/238‑1 (MQ=255)                                                                                                                                                                                                                                 gggggc‑ggAGATG‑ATT‑AAAAGTTTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGcgcg         <  1:535133/236‑1 (MQ=255)                                                                                                                                                                                                                                 aGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTT‑GCGCCGGGACCGGTGGCGGTAAAACTTTTTTTGGGGGGGACCAACGCAAGGATGATGAgggc          >  1:446023/1‑238 (MQ=255)                                                                                                                                                                                                                                          agaTG‑ATT‑AAAAGTGTTTGGCT‑GA‑TCTTTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCt  <  1:318201/238‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         CTGCTATTACAGGAGAAACCTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGC‑GGAGATG‑ATT‑AAAAGTGTTTCGCT‑GA‑TCATTGCT‑GGTTTGCTGCTTTTATTAC‑CGGGCTTTTTTACC‑GACTTCCTCGGT‑CTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGC‑‑GGCGGTTTTA‑GCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCT  >  NC_000913/4368644‑4369113

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