Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	361,668	C→T	15.1%	E65K (GAA→AAA)	lacA ←	thiogalactoside acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	361,668	0	C	T	15.1%	117.6 / 3.3	65	E65K (GAA→AAA)	lacA	thiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/34);  new base T (8/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGT  >  NC_000913/361439‑361884                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       ccGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGAGGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑G‑TCT‑Acccctc                                                                                                                                                                                                                               >  1:532777/1‑235 (MQ=255)      aacaggcttttttccAATCGGTTTCGGAAAAGAGTACATTCCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑t‑t                                                                                                                                                                                                                      <  2:225093/224‑1 (MQ=255)              tATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACACTGGTTAAATGGAAATTTGCCTAAAACTTGCGGCTTATATGGATCTTGGAACCGTAAGAGAAATAGCCAGGC‑GG‑A‑TGT‑GCCC‑A‑C‑TCG‑G‑TTTTCACTAc                                                                                                                                                                                                             >  1:478862/1‑238 (MQ=255)                      aTCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TTT‑ACCC‑A‑G‑GCG‑G‑TTTTCCCTACCgctgtcct                                                                                                                                                                                                    >  1:225093/1‑232 (MQ=255)                      aTCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCTAACCC‑A‑G‑GCG‑T‑TTTCCCC‑ACCGC‑GGGaaa                                                                                                                                                                                                   >  1:392193/1‑239 (MQ=255)                                  aaaGAGTTCATTTCGCCGTTTTTTCTAAATTCAGGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAAt                                                                                                                                                                                      <  2:691327/239‑1 (MQ=255)                                        tACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TTT‑ACCC‑C‑G‑TCG‑TGTTTCCCCCACC‑C‑GGCGAACAATTC‑TTTAATTAAGC‑t                                                                                                                                                                               >  1:61876/1‑239 (MQ=255)                                         tcttctcGTCGTTTTTTCTCAATTAATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTAACGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑TG‑T‑TCT‑ACCC‑A‑T‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑t                                                                                                                                                                               <  2:307442/235‑1 (MQ=255)                                               cGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑TCCCAA‑G‑GGG‑T‑TTTCCCC‑ACCCT‑GGCCAACACTTC‑TTTTATCAGGC‑TTTTTCtt                                                                                                                                                                        >  1:195072/1‑239 (MQ=255)                                                gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GGTT‑TTT‑CCCC‑G‑G‑GCT‑T‑TTTCCCATCCCGG‑GGCAAA‑ACTTT‑TTTAATCAAGC‑Ttttttttt                                                                                                                                                                       >  1:538148/1‑239 (MQ=255)                                                gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑TGTCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCATCCCGG‑GGAAAAAAATTT‑TTAAAACAGGT‑TTTTTCtt                                                                                                                                                                        >  2:125307/1‑239 (MQ=255)                                                    agagacagcAATTTATGGTTTACAGGGTTTCCCGTAACGGAAAGAGGAACGTTGGGTTCAATCAGTTCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑c                                                                                                                                                                  <  2:819387/231‑1 (MQ=255)                                                                       gtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GGTT‑TCT‑ACCA‑C‑G‑GGG‑T‑TTTCCCC‑ACCGC‑GGCAAAAAAATT‑TTTATATAGGCTTTTCTTTTTTTT‑CAAtt‑ttttggggggcgag                                                                                                                                               >  2:659163/1‑226 (MQ=255)                                                                       gtACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TTT‑TCCC‑C‑GCGCG‑T‑TTTCCCC‑ACCCT‑GGGCAACAATTC‑TTTAATACGGC‑GTTCTCTTTTTT‑CAACT‑TTTgtt                                                                                                                                                       >  1:819387/1‑228 (MQ=255)                                                                              gtttccctttaaggAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGAt                                                                                                                                          <  2:182373/224‑1 (MQ=255)                                                                              gtgtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTTTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGTCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTa                                                                                                                                        <  1:258886/238‑1 (MQ=255)                                                                              gtgtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGAt                                                                                                                                          <  1:540063/236‑1 (MQ=255)                                                                              gtgtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTa                                                                                                                                        <  1:659163/238‑1 (MQ=255)                                                                                gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GGTT‑TTT‑ACCC‑AGG‑GTG‑T‑TTTCCCTCCCTGG‑GGCAAA‑ATTTT‑TTTATTC‑CGC‑TTTTTCTTTTTT‑TTACt‑ttttgttgtgtgtggttgtttaaa                                                                                                                                     >  2:382240/1‑216 (MQ=255)                                                                                gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑G‑TCT‑ATCC‑C‑C‑GGG‑T‑TTTTCCCTACCCG‑GGCAAACATTTT‑TTTTATTCGGC‑TTTCCTTTTTTT‑TACCT‑TCTgttggatgacggcgttttaa                                                                                                                                      >  1:86125/1‑226 (MQ=255)                                                                                 tCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑C‑G‑GCG‑T‑TTTTCCCTCACCT‑GGCCAACAATTC‑TTTAATCAGGC‑TTTTTCTTTTTT‑TAACT‑TCTtatgggtggcgcggatttaaa                                                                                                                                     >  1:468848/1‑225 (MQ=255)                                                                                  cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TTT‑ACCC‑A‑G‑GCG‑T‑TTTTCCCTCCCGT‑GGCAAACATTTC‑TTTTATCAAGG‑TTTCTCTTTTTT‑TTACT‑TTTTCTGGATGCGCGTGGTTtaaacc                                                                                                                                   >  1:182373/1‑236 (MQ=255)                                                                                  cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTTCCCTACCGT‑GGCACAAAATTC‑TTTAATCAGGC‑TTTTTCTTTTTT‑TAAAT‑TCTGAGTGATGCGAGGGGTTAAAcac                                                                                                                                   >  1:686891/1‑237 (MQ=255)                                                                                  cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑G‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACAGT‑GGGAAAAAATTT‑TTTTATCTGGA‑TTTTTCTTTTTT‑TTACT‑TCTgtttgggggggagggttaaaaaa                                                                                                                                   >  2:540063/1‑217 (MQ=255)                                                                                   ccGTAACGGAAAGGGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCa                                                                                                                                  <  1:415999/239‑1 (MQ=255)                                                                                                        tggtgtaaaTCAGTAAGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCAAAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGa                                                                                                            <  1:382240/232‑1 (MQ=255)                                                                                                        tGGGTGCAATCAGTACGTTAACACCGATTGTTACCGTGTAGTCATCGACAATTGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGa                                                                                                            <  1:139261/239‑1 (MQ=255)                                                                                                          ggTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑TCCA‑C‑A‑GCG‑C‑TTTTCCCCACCGTGGGCAAAC‑TATC‑TTTTGTCAGGC‑CTTCTCTTTTTT‑TAACA‑TCTGTTTGATGCACGGTATTTAAATCATACACTTATCTTTTTccaccaaa                                                                                                           >  2:139261/1‑233 (MQ=255)                                                                                                                aTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑G‑TCT‑ACCC‑C‑A‑GCG‑G‑TTTCCCCTACCCT‑GGCCAACATTTT‑TTTAATACAGC‑GTTCTCTTTTTT‑TAACT‑TTTTATGGATGCGAGGGATTAAAATCCTACA‑TTAACATTTTTCCACGAAGGCtttt                                                                                                    >  1:219379/1‑239 (MQ=255)                                                                                                                      aCGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAAACAGGCGGG‑T‑TTT‑ACCC‑G‑G‑GGT‑T‑TTTCCCCCACCGG‑GCCAAAAATTT‑‑TTTTTTAACGC‑CTTTTCTTTTTT‑TTCAT‑TCTTTTGGATGCGGGGTTTTAAAATTTTAAC‑TTAACGTTTTtcccccaaacattttttttcg                                                                                              >  2:55919/1‑221 (MQ=255)                                                                                                                           aTCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TTT‑ATCC‑A‑G‑GCG‑T‑TTTTCCCCACCCT‑GGCAACAATTTT‑TTTATTTCGGA‑TTTTTCTTTTTTACCAC‑‑TCTCATAGATGCTGGGTATTAAACTCATAAC‑TTTTAGTTTTTCCCCCAAGCCTTTTTTCCGGTa‑ac                                                                                         >  2:744873/1‑237 (MQ=255)                                                                                                                           aTCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GAG‑T‑TTTTCCCTACCCT‑GGCAAAACATTA‑TTTAAACAGGC‑GTTGTCTTTTTT‑TAACA‑TTTTAAGGATGCGAGTGAATAAAATCATACA‑TTTAAGTTTTTCCACCAAGTCttttttttgggaag                                                                                          >  1:561549/1‑230 (MQ=255)                                                                                                                                  aTTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑A‑GCG‑G‑TTTCCCCTACCCG‑GGCAAACATTTT‑TTTTATTAAGC‑TTTTTCTTTTTT‑CAACT‑TTTTATGGATGCCGTGGATTTAACTCATTAA‑TTAAACTTTTCCCCCAAAGACTTTTTTTCGGTA‑AACCtttcc                                                                                  >  1:289105/1‑235 (MQ=255)                                                                                                                                  aTTGGTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAACAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑G‑TCT‑ACCC‑A‑A‑GGC‑T‑TTTCTCCTACCCT‑GGCAAACATTTC‑TTTAATACGGC‑TTTTCTTTTTTT‑CAACTCTTTCATGGA‑GCTAGAGCTTAAACTACAACATTTAACGTTTTCCCCCCAACCCTTTTTTCCGTTa‑actctgca                                                                                   >  1:897583/1‑231 (MQ=255)                                                                                                                                    tGTTACCGTGTAGTCCTCGAAAATGGTTAAATTGAAATTTGCATAAAACTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGGA‑AGCCTTCGcaca                                                                               <  1:462360/239‑1 (MQ=255)                                                                                                                                      ttcccGTTTAGTTATCGACAATGGTAAAATTGAAATTTGAATAAAAATTGCTGCCTATACGGATGTTGGAACCGTAAGAGAAATAGCCAGGG‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACata                                                                             <  2:840632/236‑1 (MQ=255)                                                                                                                                        aCCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGTC‑GG‑T‑TCT‑TCCC‑A‑C‑GGG‑G‑TTTTCCCTCCCGT‑GCCAAGAATTTT‑TTTTATTAGAC‑GTTCTCTTTTTT‑TTAAC‑TCTTTTGTTTGTGCGTGATTAAAACTATTCA‑TTAATGTTTTTCCCCCAAATATTTTTTTCGTTAGAACCTTCCCACatag                                                                            >  2:915115/1‑238 (MQ=255)                                                                                                                                             gtAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTa                                                                       <  1:145660/238‑1 (MQ=255)                                                                                                                                                   aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑a                                                                                                                                                                                                                                <  1:398987/90‑1 (MQ=255)                                                                                                                                                   aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑AACC‑a                                                                                                                                                                                                                                >  2:398987/1‑90 (MQ=255)                                                                                                                                                       aaaaTTGTTAAATTTAAATTTGCATAAAAATTGCGGCCTATATGGGTGTTAGACCCGTAAGAGAAAAAGGCAGGC‑AG‑G‑TCT‑CCCC‑A‑G‑GCT‑T‑TTTCCCCTACCGT‑TGCAAACCGTTC‑TTTAATCCGGC‑TTTCTCTTTTTT‑CAACT‑TCGGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGcc                                                           <  2:478862/237‑1 (MQ=255)                                                                                                                                                       aCAATGGTTAAATTGTAATTTGCATAAAAATGTCTTGTTATATGTATGTTGGAACCGTAATAGAAATAGACAGGC‑GG‑T‑TAT‑ACCC‑A‑T‑GCG‑T‑TTTCCCCTACCCT‑GGCAAACATTTC‑GTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATTGATGCGAGTGATTAAACTCATACA‑TTGACGTTTTCCCAATAAGTCTTTTTTCCGGTA‑AGCTTTCGCACATATCGGTAAATAGC‑TTGcc                                                           <  2:523408/239‑1 (MQ=255)                                                                                                                                                          aTGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑TCCC‑A‑G‑GCC‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGc                                                        <  1:125307/239‑1 (MQ=255)                                                                                                                                                                aaaTTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑A‑TCT‑ATCC‑A‑C‑ACGCT‑TTCCCCCCACCCG‑GGC‑AACATTTT‑TTTTATTAAGC‑TTTCTCTTTTTT‑TCACT‑TCCGATGGGTGCCAGTGATTTAAACCATACA‑TTAACGTTTTCTCCCCCAGACTTTTTTTCGGTA‑AACCTTTGCCACTATTGGTaaaaacacttcctgccctttt                                                   >  2:138898/1‑221 (MQ=255)                                                                                                                                                                       aaTTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTc                                            <  1:587684/238‑1 (MQ=255)                                                                                                                                                                       aaTTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTACCCTACCCT‑GGGAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TTTCATGGATGCGAGGGATTAAAATCATACA‑TTAAATTTTTCCCACTCAGGCTTTTTTTCGGTA‑AGCCTTCTCACAAATCGGTCAATATC‑TTGCCTCCCCTTCtttttttg                                           >  1:430713/1‑237 (MQ=255)                                                                                                                                                                           tgaaTAAAAAGTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTTATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTCTCCCACGAAGTCTTTCTTCCGGTA‑AGCCTTCGCACATATCGTTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCa                                       <  1:304911/236‑1 (MQ=255)                                                                                                                                                                               ttaaaaTTGCTGCCTATTGGGATTTTGGAACCGTAAGAGAAAAAGAACGGC‑GG‑T‑TTT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATCCA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTgg                                   <  2:883266/237‑1 (MQ=255)                                                                                                                                                                                taaatTTGGGGCCTTTTTGGATTTTGGAACCGTAAGAGAAATAGACCGGC‑GG‑T‑TCT‑TCCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGc                                  <  2:392193/238‑1 (MQ=255)                                                                                                                                                                                          gCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAt                         <  1:321502/238‑1 (MQ=255)                                                                                                                                                                                           ccTATATGGATGTTGGAACCGTAAGAGAAATAGACAGTC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGcaca                                                                               <  1:916212/184‑1 (MQ=255)                                                                                                                                                                                               tatGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGAt                                                                                                                                          <  1:555325/123‑1 (MQ=255)                                                                                                                                                                                               tatGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGAt                                                                                                                                          <  1:555314/123‑1 (MQ=255)                                                                                                                                                                                                atGGGTGGTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                  <  1:818064/239‑1 (MQ=255)                                                                                                                                                                                                atGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCa                  <  1:184058/239‑1 (MQ=255)                                                                                                                                                                                                  ggATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                <  1:138898/239‑1 (MQ=255)                                                                                                                                                                                                  ggATGTTGGAACCGGAAGAGAAATAGACCGGC‑GG‑T‑TTT‑TCCC‑A‑G‑GGG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                <  2:289105/239‑1 (MQ=255)                                                                                                                                                                                                      gTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑T‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg             <  1:756615/238‑1 (MQ=255)                                                                                                                                                                                                      gTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTGCCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg             <  2:219379/238‑1 (MQ=255)                                                                                                                                                                                                      gTTGGAACCAGAAGGGAAATAGACAGGC‑GG‑T‑TCT‑GCCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg             <  1:191546/238‑1 (MQ=255)                                                                                                                                                                                                        tGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑T‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAAGCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACGCCGTTa          <  1:866476/239‑1 (MQ=255)                                                                                                                                                                                                         ggAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑AACC‑A‑A‑GCG‑T‑TTTTCCCCAACCT‑GGCAAACATTTCTTTTAATC‑GGC‑TTTTTCTTTTTT‑TAAAT‑TCTGGTTGATGCGAGTGGTTTAAATCATACA‑TTTACGATTTTCCCCGACGTCTTTTTTCCGGTT‑AGCCTTCCTCCAATACGGTAAAAAACATTGCCT‑CTCCTTTTATTTCTGGCATTTGCCGTTTCAATGCGATCACTCCGTTAt         >  1:432513/1‑239 (MQ=255)                                                                                                                                                                                                                 aaGAGAAATAGACCGGC‑TG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCGCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGt  <  1:744873/238‑1 (MQ=255)                                                                                                                                                                                                                 aaGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGc                                  <  1:326853/206‑1 (MQ=255)                                                                                                                                                                                                                      aaaTAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCtttttt                                                                                                                                                                    <  1:872007/76‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       CCGATCCAGACGTTATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGC‑GG‑T‑TCT‑ACCC‑A‑G‑GCG‑T‑TTTCCCCTACCGT‑GGCAAACATTTC‑TTTAATCAGGC‑TTTCTCTTTTTT‑CAACT‑TCTGATGGATGCGAGTGATTAAACTCATACA‑TTAACGTTTTCCCACGAAGTCTTTTTTCCGGTA‑AGCCTTCGCACATATCGGTAAATAGC‑TTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGT  >  NC_000913/361439‑361884

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