Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,368,930	C→T	11.5%	G89G (GGC→GGT)	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,930	0	C	T	11.5%	152.8 / 2.1	68	G89G (GGC→GGT)	fxsA	suppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base C (24/35);  new base T (2/6);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTA  >  NC_000913/4368704‑4369166                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCCGGGGT‑GATTAAAAATGTTTTGCCGTATCTT‑TCCGGGTTT‑CTGC‑CTTTTTTACCCGGGTTTTTTTAcca                                                                                                                                                                                                                                    >  1:289546/1‑237 (MQ=255)        attttCAATCCTTATTCAGGTTGCCCATGGCTTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCTGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCg                                                                                                                                                                                                                           <  1:612092/237‑1 (MQ=255)         ttttCAATCTTTTTTCAGTTTGCCCATTTATTGGGGGTTTTGCTGCCCCTCGTGCTGGTTTATTTCACGTCAGTTACCGGTATGTAACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCg                                                                                                                                                                                                                           <  2:982989/237‑1 (MQ=255)            tCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAACAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGtct                                                                                                                                                                                                                       <  1:958221/239‑1 (MQ=255)                     aTTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAAGCCAGGGGGGGGGGT‑TGTTCAAAAAGTTTCCCTG‑ATCGT‑ATCTGGGCTGCTGCTTTTTATTT‑CCGGCTTTTTTTACGCAC‑TCCTTGGTGTTTTGCtttt                                                                                                                                                                                                              >  2:612092/1‑239 (MQ=255)                             tggccATGTATGGGGGGTATTGCTGGCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                                                                      <  1:580658/236‑1 (MQ=255)                              gCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGAT‑GATTAAAAGTGTTTCGCTG‑GTCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCCCGGGCTTTTACTTTTATTtccgcc                                                                                                                                                                                                     >  1:506209/1‑239 (MQ=255)                                     ttttGGGGGTTTTGCTTACCCTCGTTCTTGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTgccgc                                                                                                                                                                                                      <  1:401503/229‑1 (MQ=255)                                            ggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCg                                                                                                                                                                                                                           <  1:114421/203‑1 (MQ=255)                                               tttgcttgcCCTCGTTGTGGGTTTATTCACGTTAATTATCGGTAAGTCACTGGTAGGGAACCAGGGTTTTAAGAATTTCGTGCTGATGAAGAAAAAAAGTGCGGCGGGTGAATACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TCTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGa                                                                                                                                                                                    <  2:386996/238‑1 (MQ=255)                                               aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTGAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGGTGGATTAAAAG‑GGTTCG‑TG‑TGCCTATGCTGGGTTTTTGC‑TTTTTTTA‑CCGGGGTTTTTTTCCGCCTTCCCTGGGCTTTCACTTTTTTTTGCCCCGGGGTaaaaaagagtat                                                                                                                                                                                    >  2:401503/1‑231 (MQ=255)                                               aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATATCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGGTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTTTTA‑CCGGGCTTTTTTACCGACTTCCTCGGGCTTCTACTTTTTTTTCCGCCGGTGCAAAAGCATCtta                                                                                                                                                                                    >  1:128236/1‑237 (MQ=255)                                                      aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCACCGGCGGAGAT‑GAGTAAAAGTGGTTCACTG‑ATCAT‑TGCTGGGTTTCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGGGCAAAAGCATCttgcagggca                                                                                                                                                                             >  1:821052/1‑230 (MQ=255)                                                        ccTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGGTTTCTCCCTTTTTTACCGCCGGGGCAAAAGAATCTGGCAGTGAAGt                                                                                                                                                                           >  1:149061/1‑239 (MQ=255)                                                           cGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGGTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTTCCGGCTTCCTCCGGCTTCTAATTTTATTGCCGCCCGGGCCAAAGAACTTTACAAGGAAGTTgt                                                                                                                                                                        >  1:689488/1‑238 (MQ=255)                                                              gggggTTATTTTCCCGTAATTTTTCGGTAGTTCACTGGTACGGAACCAGGGTTTTAAGAATGTCGTGGTGATGCAGCAAAAAATGGCGGCGGGTGAAATCCCAGCGGTGGAGAT‑GATTAAAAGTTTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACGGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGc                                                                                                                                                                     <  2:816581/236‑1 (MQ=255)                                                               cTGGTTATATGCTCGGCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGAGGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGc                                                                                                                                                                     <  1:295013/238‑1 (MQ=255)                                                                            cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGGGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGGCTTCTACTTTTTTTGCCCCCGGGGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCttt                                                                                                                                                       >  1:874210/1‑239 (MQ=255)                                                                             gTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGGGTTTCGCTG‑ATCAT‑TGCTGGGTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTTCCGACTTCCTCGGGCTTCTACTTTTTTTGCCGCCGGTGCAAAAGAATCTGACCTGGAAGTTGATGCCCCATTTGCGCtttt                                                                                                                                                      >  1:104596/1‑239 (MQ=255)                                                                                aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGAGAT‑GATTAAAAGGGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAACATCTGACAGTGAAGTTGATGCCGCATTTGCGCtttttt                                                                                                                                                    >  1:743898/1‑236 (MQ=255)                                                                                  ttAGCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTgcg                                                                                                                                                           <  1:330361/229‑1 (MQ=255)                                                                                   tagcGGTATGTCACTGGTACGTAGCCAGGGCTTTAAGAATTTCGTTCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCg                                                                                                                                                  <  1:1025050/234‑1 (MQ=255)                                                                                      cGGTATGGCACTGGTACGTAACCAGGGCGTTAAGAATTTCGTGCGGAGGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGATG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGc                                                                                                                                             <  1:184103/239‑1 (MQ=255)                                                                                            gTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTtgcg                                                                                                                                       >  1:409661/1‑239 (MQ=255)                                                                                              cccTGGTCCGTAACCAGGGCTTTAAAAATTTCGTGCTGATGCAGCAAAAAAGGGGGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCCCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTggcggc                                                                                                                                     <  2:689488/237‑1 (MQ=255)                                                                                               aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGGTTCGCTG‑ATCAT‑TGCTGGTTTTCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGGCTTTTACTTTTTTTTCCGCCGGTTCAAAAGCATTTGACAGGGAAGTTGATTCCGCATTTTCGCTTTTTTCGCATGCCTGggggcg                                                                                                                                    >  1:814595/1‑239 (MQ=255)                                                                                                              gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTTTTGCCGCCGGGGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGGGGTTTTAGGCCCggg                                                                                                                      >  1:354980/1‑238 (MQ=255)                                                                                                              gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTTCGCTTTTTTCGCCTGCCTGGCGGCGGGTTTatggcccgg                                                                                                                      >  1:620694/1‑230 (MQ=255)                                                                                                               ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGGGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGGGAAGTTGATGCCGCATTTTGGCTTTTTTCCCAGTCCTTGGGGCGGGTTTTACGCCgggg                                                                                                                     >  1:681531/1‑237 (MQ=255)                                                                                                                gCTTTAAGAATTTCGGGCTGAGACAGCAAAAAATGGCGGCGGGTAAAACCCCAGGGGGGGAAAT‑GTTTAAAAGTTTTTCGCTT‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCggg                                                                                                                      <  2:681531/236‑1 (MQ=255)                                                                                                                   tttaaattttCGTGCTGATGCAGCAAAAAATGTGGGGGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGGTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCCGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                                                                                <  2:409661/232‑1 (MQ=255)                                                                                                                   ttAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                                                                                <  1:69230/239‑1 (MQ=255)                                                                                                                   taaaaattttCGTGGTGATGCAACAAAAAAGGGCGGGGGGTGAAAACCCAGCGGGGGAGAT‑GAATAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                                                                                <  2:104596/232‑1 (MQ=255)                                                                                                                        aaTTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCgg                                                                                                           <  1:778415/239‑1 (MQ=255)                                                                                                                          ttttggccTGGTGCAGAAAAAAATGGGGGGGGGTTAAAACCCGGGGGCGGAGAT‑GATTAAAAGTTTTTCGCTT‑GTCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTa                                                                                                         <  2:620694/232‑1 (MQ=255)                                                                                                                                   gATGCAGAAAAAAATGGCGGCGGGTGAAAACCCAGCGGGGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTg                                                                                                <  2:743898/239‑1 (MQ=255)                                                                                                                                      gcagcaAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATg                                                                                             <  1:470706/239‑1 (MQ=255)                                                                                                                                      gcagaaAAAAATGGGGGGGGGTGAAAACCCAGGGGGGGAGAT‑GATTAAAAGTTTTTCGCTG‑ATCAT‑TGCGGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATg                                                                                             <  2:874210/239‑1 (MQ=255)                                                                                                                                        agcaAAAAAGGGCGGCGGGTGAAAGCCGGGCGGCGGAGAT‑GATTAAAAGTTTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACGTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATg                                                                                             <  1:667797/237‑1 (MQ=255)                                                                                                                                         gtaaaaaaGGGCGGCGGTTGAAAACCCAGGGGCGGGGAT‑GATTAAAACTGTTGCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                           <  1:433411/236‑1 (MQ=255)                                                                                                                                             aaaaTGGCGGCGGGTGAAAACCCAGGGGCGGAGTT‑GATTAAAAGGTTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGt                                                                                                                <  1:337968/213‑1 (MQ=255)                                                                                                                                                                   cAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCt                                                                                                                                                                                      >  1:369029/1‑121 (MQ=255)                                                                                                                                                                   cAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCt                                                                                                                                                                                      <  2:369029/121‑1 (MQ=255)                                                                                                                                                                    aGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGTGGCGGTTTTAGCGCCGGGACCCGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGTTGATGAacgc                                                                >  1:666102/1‑238 (MQ=255)                                                                                                                                                                     gggggggAAAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                              <  1:10478/234‑1 (MQ=255)                                                                                                                                                                           gagaT‑GATTAAAAGTGTTTCGCTT‑ATCAT‑TGTTGGTTTGCGTC‑TTTTATTA‑CCGGGCTTTTTTACCGTTTTCCTCGGGTTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGcc                                                         <  1:769087/238‑1 (MQ=255)                                                                                                                                                                           gagaT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                                                                              <  1:611065/217‑1 (MQ=255)                                                                                                                                                                            agaT‑GATTAAAAGTGGTTCGCTG‑ATCAT‑TGCTGGTTTGCTTC‑TTTTATTA‑CCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGAATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCt                                                        <  1:819376/238‑1 (MQ=255)                                                                                                                                                                               t‑gatTAAAAGTTTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTC‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGa                                                     <  1:88871/238‑1 (MQ=255)                                                                                                                                                                                       aaGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGCCGGCGGTTTTAGGGCCGGGACCGGGGGCGGTAATACTTTTGATGGGTAGTACCAGCGAAAGGAGGATGTGCCCGACCCCCTTGTTCATAAAg                                             >  1:875635/1‑239 (MQ=255)                                                                                                                                                                                                gCTG‑ATCGT‑TTGTGGTTTTTTGT‑TTTTATTA‑CCGGGTTTTTTTCCCGCCTCCCTCGGTTTTCTTCTTTTATTGCCGCCGGTGCAAAAGCATCTGCCAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                    <  2:149061/239‑1 (MQ=255)                                                                                                                                                                                                gCTG‑ATCAT‑TGTTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGcc                                    <  1:437752/239‑1 (MQ=255)                                                                                                                                                                                                     aTCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGGTCCTAAAGACGATCGCCCGGa                                >  1:220769/1‑239 (MQ=255)                                                                                                                                                                                                     aTCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAAACCATCGCCCGGa                                >  1:12329/1‑239 (MQ=255)                                                                                                                                                                                                                gTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTACTACTTTTGATGTTGAGTACCAACGCAAGGATGATGAGCGCGAACGCCTTGATCATAAAGACGATTGCCAGGATTAATGTCGa                      >  1:108959/1‑239 (MQ=255)                                                                                                                                                                                                                gTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                                                       <  1:507039/206‑1 (MQ=255)                                                                                                                                                                                                                gTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                                                       >  2:507039/1‑206 (MQ=255)                                                                                                                                                                                                                  ttttttgtttttttTA‑CCGGGCTTTTTTACCGACTTCCCCGGTCTTCTACTTTTTTTCGCCCCGGTGAAAAAGCATCGGACAGTGAAGTTGATGCCGCAGTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGa                      <  2:814595/230‑1 (MQ=255)                                                                                                                                                                                                                   tggtgc‑tTTTATTA‑CCGGGTTTTTTTACCGATTTCCTCGGTCTTCTACTTTTATGGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAAc                   <  1:102697/236‑1 (MQ=255)                                                                                                                                                                                                                   tgctgc‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATTGTGAGTACCAGCGCAAGGATGATGAGCGCGACCGCCTTGTTAATAAAGACGACGGCCCGGATTAATTTCCaaa                    >  1:820086/1‑238 (MQ=255)                                                                                                                                                                                                                     ctgc‑tTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGGAGTACCAGCGAAAGGCTGTTGAGCGCGACCGCCTTGATCTTAAAGCCGATCGCCAGGATTAATGTCGAAACg                  >  1:585077/1‑238 (MQ=255)                                                                                                                                                                                                                     ctgc‑tTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTAGGGCCGGGGCCGGTGGCGGTAATACTTTTGATGGGGGGTACCAGCGAAAGGGTGGTGGGCGCGAACGCCTTGATTATAAAGACGATCGCCAGGATTAATgggcaaaacc                 >  2:272673/1‑230 (MQ=255)                                                                                                                                                                                                                        ttttttttTC‑CCGGGTTTTTTTCCCGCTTTCCTCGGTTTTTTTCTTTTTTTGCCCGCGGTGCAAAGGCATCTGACGGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCg               <  2:12329/237‑1 (MQ=255)                                                                                                                                                                                                                        c‑tTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGt                                                                                       <  1:949460/166‑1 (MQ=255)                                                                                                                                                                                                                           tttATTC‑CCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGAAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACTAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGa             <  1:803747/238‑1 (MQ=255)                                                                                                                                                                                                                             tattt‑CCGGGCTTTTTTACCGATTTCCTCGGTTTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTa          <  1:497296/239‑1 (MQ=255)                                                                                                                                                                                                                                ttcccGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTAtg        <  1:272673/236‑1 (MQ=255)                                                                                                                                                                                                                                      ggCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACACCGGATTATTTGGTTa  >  1:995825/1‑239 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             ATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGAT‑GATTAAAAGTGTTTCGCTG‑ATCAT‑TGCTGGTTTGCTGC‑TTTTATTA‑CCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCGCCAGGATTAATGTCGAAACGCCGGATTATGTGGTTA  >  NC_000913/4368704‑4369166

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