Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,977,062	T→C	15.0%	G26G (GGA→GGG)	motA ←	proton conductor component of flagella motor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,977,062	0	T	C	15.0%	160.1 / 3.0	68	G26G (GGA→GGG)	motA	proton conductor component of flagella motor
Reads supporting (aligned to +/- strand):  ref base T (39/17);  new base C (2/9);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGT  >  NC_000913/1976825‑1977299                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              aCGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCCGCGGGTTGATAGA‑GGGCTc                                                                                                                                                                                                                                              >  1:393725/1‑239 (MQ=255)         tttCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAAA‑GTGCTCCAAGGCt                                                                                                                                                                                                                                       >  1:872010/1‑238 (MQ=255)             aaactccccttCCAGCGAAAAACTCCCATACTGCCGAGATTTCGCCATCACCCGATAAAGCAGAGCCAGCAGAGCCATATACCTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGAACCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCAc                                                                                                                                                                                                                                  <  2:143587/230‑1 (MQ=255)                 tctcGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGGAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGt                                                                                                                                                                                                                              <  1:306267/236‑1 (MQ=255)                 tcccGTTCAAACAAAAACACCCCCATCTGCGGCGCGTTCGCCTGCACCCGATAAAGCCAGGCCAGAAGTACCATATACATTGCATTGGTGTATTTGGGGCGACGAAAAAGCAACGGCA‑GCGCCTACAGCGT‑GC‑CTTTAACCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAACCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGt                                                                                                                                                                                                                              <  2:741396/236‑1 (MQ=255)                 tCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACAAGTTAAGAGGGTTGATCGC‑GAGCTCCAAGGCTTCCACCGGt                                                                                                                                                                                                                              >  1:307917/1‑239 (MQ=255)                    cGTTCCAGCGTAAACATCCCCTTCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCCTATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTCCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTCTTGCCGACGATAAAAGAACCACACCCTGCAACGCCAGAAATCACCAGTTCAGCGGGTAGTAAGA‑GGTCACCAAGGATGCCACCgggaca                                                                                                                                                                                                                           >  2:975937/1‑235 (MQ=255)                     gTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCATTG‑CCATTCTTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTAATAGA‑AGGCCCCAACGCTTACCCCGGgcacc                                                                                                                                                                                                                          >  2:451251/1‑237 (MQ=255)                                     ccccATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTc                                                                                                                                                                                                                                     <  1:1054300/212‑1 (MQ=255)                                     ccccATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCCCCAAGGCTTc                                                                                                                                                                                                                                     >  2:1054300/1‑212 (MQ=255)                                               cgcgATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgag                                                                                                                                                                                                >  1:1040954/1‑239 (MQ=255)                                               cgataTTTCGCAATCAACAGAGAAAGCAGCGCCAACAGACCCATATTCGTTGCTTTGGTGTCTTTGGGGCGACGAAGCAGCAACGGAA‑GCGCCTTCAGCTT‑GC‑CTTTACTCGCTTTG‑CCATTATGGCGGACGATAACCGACCCAATCCCGGCACCGGCAAAAAGCACCATTTCAGCGGGTTGATAGA‑GTGCTCAAAGCCTGCCGCCGGGCATCAAATAACCGCCGAAAACTGTACCgag                                                                                                                                                                                                <  2:675005/235‑1 (MQ=255)                                                                             acgacgaaatcctaaacgaTTGTTTTGTTGTATTGGGGACGAGAAAACAGCAGCGGCG‑CCGCCGTCAGCGT‑GC‑ATTTAAGCTCTTTG‑CCTTTTTTCCCGACGATAACCCACCCACTCCCTGCCCCGGAAAGAATCACCAGTCCAGCGGGTTGATAGA‑GTGCCCCAAGCCTTCCACCGGTCATAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCa                                                                                                                                                                   <  2:669608/220‑1 (MQ=255)                                                                                 gcagATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGCTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGACCCGAGAACAACCAGGAAAACTAATAAAAAAAGAAAGAc                                                                                                                                                               >  2:58410/1‑238 (MQ=255)                                                                                    gctccatatactgttgttgtgtgtTTTTGGGGCGACGAAACAGCAACGGCG‑TCACCTTCGTAGT‑CC‑TGTTACTGTCGTAC‑CCATAATTGCCGACAATTACCGCCCCACTCCCTGCACCGGTAATAATCCCCAGTTCAGCGGGTGAAAAGG‑GTGCCCCAAGGCTCCCCCCGTTTAAAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGAcat                                                                                                                                                             <  2:166073/217‑1 (MQ=255)                                                                                     accctatatctttgttgtGGTTTTTTGGGAGCCAAGAAACCACAGGGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTGT‑CCATTATTGCGGCCGAAAAAGCTCCCAACCCCGGCCCGGGCAATAATCGCCAGTACAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGTTAATAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGAcatca                                                                                                                                                           <  2:131294/221‑1 (MQ=255)                                                                                          tttctttttttttttgttttttGGGGCGCCAAAAC‑GCAAAGGCACGCCCCTTC‑GCGTGGC‑CTTAAATCGCTTTT‑CATTTATGCCCGACGAAAAACCCACCAACCCCTGCCCCGGCAATAATCACCAGTTCAGCGGGTTAAAAGA‑GTGCTCCAAGGTTTCCACGGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTc                                                                                                                                                     <  2:441616/220‑1 (MQ=255)                                                                                            taCATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAACACGACATCATCCTTcc                                                                                                                                                    >  1:589697/1‑238 (MQ=255)                                                                                               gtttttttGTTGTTTTTGGCGCGACGAAACAACAACGGGC‑GCGCC‑TCCGGGC‑GCTCTTAAATCCTTTTG‑CCATTACTGCCGACGAAAAACGCCCCAATCCCTGCACCGCCAAACACACCCAGTGCAGCGGGTTGAAAGA‑TGGCCCAAAGGCTCCCACGGTTCATAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTg                                                                                                                                                <  2:191755/234‑1 (MQ=255)                                                                                               aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTcc                                                                                                                                                    >  1:742821/1‑235 (MQ=255)                                                                                               aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACt                                                                                                                                                 >  1:131294/1‑238 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACCACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGtt                                                                                                                                              >  1:990203/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAgg                                                                                                                                                                                      >  1:190150/1‑199 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTACTAAGATAAGCACGAAATCATCCTTCCAATGtt                                                                                                                                              >  1:522678/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCCCTGtt                                                                                                                                              >  1:441616/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGGAACCTAAAAAGAAAAGAAAGACAACATCCTTCCACTGtt                                                                                                                                              >  2:996415/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGAAACCTAAAAAGAAAAAAAAAAAATAACCCTTCCACTGtt                                                                                                                                              >  2:365367/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACGGTACCGAGAACAACCAGGTAACCTAAAAAGATAACCACGAAATCAACCTTCCACGGtt                                                                                                                                              >  2:689892/1‑239 (MQ=255)                                                                                                 tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GGGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGCAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGCCATCATCCTTCCACTGtt                                                                                                                                              >  1:1093734/1‑239 (MQ=255)                                                                                                   cTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTg                                                                                                                                             >  1:267259/1‑238 (MQ=255)                                                                                                      tGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGAAACCTAATAAGATAAGAACGACATCATCCTTCCACTATTAACCa                                                                                                                                         >  2:979118/1‑239 (MQ=255)                                                                                                           tATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCAc                                                                                                                                                                                                                                  >  1:391829/1‑145 (MQ=255)                                                                                                           tATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCAc                                                                                                                                                                                                                                  <  2:391829/145‑1 (MQ=255)                                                                                                             tttGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCCGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATACCCGCCGAAAACTGTACCGAGCACAACCAGGTAACCTAATAAGATAAGCACGACATCACCCTTCCACTTTTggaaagaaaagg                                                                                                                                  >  2:821283/1‑228 (MQ=255)                                                                                                             tttGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGGGAACAACCAGGAAACCCAAAAAGATAAGCACAACATCATCCTTCCACTGTTGACCAtaaacaa                                                                                                                                  >  2:875400/1‑233 (MQ=255)                                                                                                             ttgggAGCGACGAAACGGCACGGGCA‑CGCCCTCCGGGGT‑CC‑CTTTAATCGTTTTG‑CAATTATTGCGGACGATAAACGACCCAATCCCGCCCCCGGAAATAATCCCCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAgg                                                                                                                                  <  1:712125/236‑1 (MQ=255)                                                                                                                   gCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTAGTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGC‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCa                                                                                                                                                                   >  1:451548/1‑200 (MQ=255)                                                                                                                      acgaAACAGCACCGGCA‑GCGCCTTCAGCGG‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGCCCCAACCCCGGCACCGGCAATAATCCCCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGTTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAg                                                                                                                          <  1:689892/238‑1 (MQ=255)                                                                                                                                        gcgcCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCAACCTTCCACTGTTGACCATGAAAGGATTTTAAGTAGGAAGGCGTTAAcgcg                                                                                                        >  2:127630/1‑239 (MQ=255)                                                                                                                                                   gcttccttttAATC‑CTTTGTCCATTTTTGCCGAGGATAAACCACCAAACCCCTGCACCGGCAAAAATCACCAGTTCAGCGGGTTGAAAGA‑GTGCTCCAAGGCTTCCACCGGTCAACAAAAAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTgggg                                                                                               <  2:267259/232‑1 (MQ=255)                                                                                                                                                       c‑cTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCACTAATCACCAGTTCAGCGGGTTGCTAGA‑GTGCTCCCAGCCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGACCAACCCGGTAACCTAATCAGATCAGCACGACCTAAACCTTCCACTCTCAACCATAACGGAGTGTTCATGCGTCAAGCTGTACCCCGCGATTagggaacaa                                                                                          >  2:870815/1‑234 (MQ=255)                                                                                                                                                              aaaggacagccaattattGCCGACGATAAACGACCCAACCCCTGCCCCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAag                                                                                     <  1:35828/227‑1 (MQ=255)                                                                                                                                                                      tccctttgttgcgagcgaAAAACGACACAACCCCTGCCCGGGCAATAACCACAACTACAGCGGGTTGATAGA‑TTGCTCCAAGGCTTCAACGGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTa                                                                             <  1:699490/224‑1 (MQ=255)                                                                                                                                                                                  cgacgaTAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTTACCATGACAGGATGTTAAGTCGTCAGGCGTTAACGCGGGATTGGGGCAAAAAAAAAAAGGGGTAcg                                                                           >  2:35828/1‑229 (MQ=255)                                                                                                                                                                                     cgaTAAACGATCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGACATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATACGCAAGACATCATCCTTCCACTGTTGACCAGGACAGAATGTTCAGCCGTCAGCGGTTAACGCCGGATTGGGGCAAAAACAAAAAGCGGAACTGCGTTACCgaggt                                                              >  2:698608/1‑235 (MQ=255)                                                                                                                                                                                                ccAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTACCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATATTCAGTAGTCATGCGTTAACGACCGATTGGGGCAAAAAAAAGAAACGGAACGGCGTTACCGCAGCTGGAATGGTgg                                                   >  2:149116/1‑238 (MQ=255)                                                                                                                                                                                                         gCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGCAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGCTACCGCTGCTGGAATGTTGCGCCTCACCg                                          >  1:786267/1‑239 (MQ=255)                                                                                                                                                                                                         gCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGCCGTTACCGCTGCTGGAATGTTGCGACCCACCg                                          >  2:267473/1‑239 (MQ=255)                                                                                                                                                                                                          cACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCCGACCga                                         >  1:207045/1‑238 (MQ=255)                                                                                                                                                                                                          cACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTAAACGCGCGATTGGGGCAAAAAAAAGCACCGGTACGCCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                         >  1:827009/1‑239 (MQ=255)                                                                                                                                                                                                           ccgggaaataatCACCAGTTAAGCGGGTTGATGGG‑TTTCCACAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                         <  1:64761/232‑1 (MQ=255)                                                                                                                                                                                                           cacggaaaaaaaccacAGTTAAGGGGGTTGAAAGA‑GTCCCCAAAGGCTCCCACCGGTAACAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                         <  1:996415/223‑1 (MQ=255)                                                                                                                                                                                                             tggCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAag                                                                                     >  1:569047/2‑192 (MQ=255)                                                                                                                                                                                                                   tcacccccAGTTCAGGGGGTTTATAGA‑GTGCCCCAAGGCTCCCCCGGGTCATAAAATACCCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                <  1:127630/233‑1 (MQ=255)                                                                                                                                                                                                                   taatCACCATTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCAACAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                <  1:875400/239‑1 (MQ=255)                                                                                                                                                                                                                   taacCACCAGTCCAGGGGGGTGATAGA‑GTGCCCAACGGCTTCCACCTGTCAAAAAAAACCCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                <  1:267473/239‑1 (MQ=255)                                                                                                                                                                                                                     atcccacgTTCAGGGGTTTGATAGA‑GGGCTCCAAGCCTTCACCCGGTCAACAAATATCCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                              <  1:574574/232‑1 (MQ=255)                                                                                                                                                                                                                     atCACCAGTCCGGGGGTTTAATAGA‑GGTCCAAAAGGTTTCCCCGGGACATCAAAAACCCGACGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                              <  1:979118/239‑1 (MQ=255)                                                                                                                                                                                                                     accACCAGTCCAG‑GGGTTAAAAGAGGTCCCCCAGGGCTCCCCCCGTTCACAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                              <  1:451251/237‑1 (MQ=255)                                                                                                                                                                                                                         ccAGTTAAGCGGTTTGATAGA‑GTGCTCCAAGGCTTCCCCCGGTCAACAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                          <  1:365367/239‑1 (MQ=255)                                                                                                                                                                                                                                  cGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCt                                               >  2:679253/1‑209 (MQ=255)                                                                                                                                                                                                                                  cGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCt                                               <  1:679253/209‑1 (MQ=255)                                                                                                                                                                                                                                       tGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTGGGCCCTCACCGTATCAGTTAAACAGCCTTTACTCCCGGtt             >  1:719914/1‑238 (MQ=255)                                                                                                                                                                                                                                                tGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCCACAAAAAGCAGCGGCACGACGTTACCCAGGCAGGAAAGTTGCCCCCCACCAAATAAGTTAAAAAGCCTGTaatacatgttacacaagcg    >  2:294901/1‑221 (MQ=255)                                                                                                                                                                                                                                                tGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAACAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCCCCGTATCAGTTAAACAGCCTGAACTCTCTGTTCATCcagca    >  1:425389/1‑239 (MQ=255)                                                                                                                                                                                                                                                tGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACGGCCTGTACTCTCTGTTCATc         >  1:51567/1‑234 (MQ=255)                                                                                                                                                                                                                                                tGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGACAGGCGTTAACGCGCGATTGGGGCAAAAAAAAACAGCGGTACGTCGTTACCCCTGCTGGAATGTGGCGCCCCACCGTATCAGTTAAACACACTGTACTCCCTGATCATCcagca    >  2:993020/1‑239 (MQ=255)                                                                                                                                                                                                                                                  cTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGt  >  1:543184/1‑239 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              ACGGGGATTTTCAATATCACGTTCCAGCGAAAACATCCCCATCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCA‑GCGCCTTCAGCGT‑GC‑CTTTAATCGCTTTG‑CCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGA‑GTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTAACCTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGT  >  NC_000913/1976825‑1977299

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