Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,186,332	T→G	18.1%	G163G (GGT→GGG)	pepT →	peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,186,332	0	T	G	18.1%	144.8 / 2.1	70	G163G (GGT→GGG)	pepT	peptidase T
Reads supporting (aligned to +/- strand):  ref base T (21/34);  new base G (13/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGG  >  NC_000913/1186116‑1186555                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   tcccaaattttttgAAAATTATCGGGGGGGCGATATTGGGCGGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTCCCCGGTGCTGCATCTGCTGCTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGt                                                                                                                                                                                                                         <  2:340363/217‑1 (MQ=255)  ccGCAAATTGGTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTTTCGGGGATGAAGTTTTAA‑CCCCGGTTATGTTCCCGGTTCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGt                                                                                                                                                                                                                         <  1:79897/227‑1 (MQ=255)           gttgAAACCTTTCGCGGTGGCGGTGTCGCGCT‑GGTTTCTGGTGTGTAATGTTGTCCCCCGGTTATGTTCCCGGTTATGCATCAGCTACTGGGTCAGACCCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATTACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCAGGGTGATATTCGCGTCGCCTTTAcc                                                                                                                                                                                                               <  2:661957/228‑1 (MQ=255)            ttgAAATCTATCGGGGTGGCGATATTGCGCGGGGTATCGGCGATGAAGTTTTAT‑CACCGTTTAGGTTCCCGGTGCGGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTAcc                                                                                                                                                                                                               <  2:247507/227‑1 (MQ=255)             ttaaaaCTTTCGCGGTGGTGATTTTGCGGTGGGTATCGGGGATGAAGTTTTTT‑CCCCGGTTAGGGTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTAcccc                                                                                                                                                                                                             <  2:353232/226‑1 (MQ=255)             tggatgcttgtgcgGGGGCGATGGGGCGCTGGGTATCGGGGGTGAGGGTTTAT‑CGCTGGTTATGTCGCCGGTGGTGGAGCAGGTACGGGGTCAGGCGCTGATTAGCACCGATGGTAAAACCTTGGTAGGGGCCGATGACAAAGAAGGTATTGCAGAAATGATGGCCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATGCCGCAGGGTGATATTCGCGTCGCCTTTGcccc                                                                                                                                                                                                             <  2:231760/217‑1 (MQ=255)                  aCTATCGCGGTTGGGATATTGCGCTGGGTATCGGCGATGAAGTTTTAT‑CCCCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑Tg                                                                                                                                                                                                       <  1:619453/228‑1 (MQ=255)                       ggcctgtgCGTTATTGGGCTGGGTATGGGAGATGAATTTTTAG‑TACCTGTAATCGTCCCGGTCAGACATTAGGTACTGGGTCAGACCGTTATTACCACCGAAGGTAAACCCTTGTTAGGTGCCGATGACAAAGCAAGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑Tgaagaa                                                                                                                                                                                                  <  2:613156/221‑1 (MQ=255)                           gTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGGGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTCAAACCTTGTTCGTTGCCGAGGACAAAGCAGGGATTGCAGGAATCATGACCGCACTGGCGGGGTTGCA‑ACAAAAAAAAAAGTCCG‑AAGGGGATATTCGCGTCTGCTTTACCCCGGC‑TGAAGAAA‑Tgg                                                                                                                                                                                             >  2:429244/1‑228 (MQ=255)                                        gcTGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGGTTTGCACAC‑A‑AAAAAAATTCCCAAAGGGGAAATTCCGGTGCCCTTTTCCCCGGAGTGAA‑AAG‑AGGGCAAAgggggaaa                                                                                                                                                                                >  2:9375/1‑224 (MQ=255)                                         cTGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATCACCACCGATGGTAAAACCTTGTTAGGGGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGGATTGCA‑AA‑AGGAAAAAATTCCGAATGGGGGAATTCCGGTCGCCTTTACCCCCGA‑TGAAAAAA‑GGGGCAAAaggggcaaa                                                                                                                                                                               >  2:147100/1‑223 (MQ=255)                                            ggTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGGGATATTCCGGCCGCCTTTACCCCGGA‑GGAAGAAG‑GGGGACAACGGGCGGAACa                                                                                                                                                                             >  1:211759/1‑227 (MQ=255)                                            ggTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGAATTGCA‑AC‑AGAAAAAAATTCCGCATGGGGACATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑AGGGCAAAGGGGCGAAACa                                                                                                                                                                             >  1:579722/1‑227 (MQ=255)                                            ggTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAAACTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑CGAAAAAAAATCCGCCTGGGGATCTTCTCGTCGCCTTTACCCCCGA‑AGAAAAAA‑TGGGAAAAAGGGCGAAACAt                                                                                                                                                                            >  2:19353/1‑228 (MQ=255)                                             gTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTAcccc                                                                                                                                                                                                             <  1:522480/196‑1 (MQ=255)                                                     gATGAAGTTTTAT‑CCGCGGTTATTGTCCCGGTGCCGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGtt                                                                                                                                                                   <  1:248858/228‑1 (MQ=255)                                                      aTGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTCCTGGGTCAGACGCTGATTACCACCGATGGTCAAACCTTGTTAGGTGCCGATGACCAAGCAGGTATTGCCGAAATCATGACCGCGCTGGCGGTATTGGA‑AC‑AAGAAAAAATTCCGACTTGGGAATTTCGCGTCGCCTTTTCCCCGGG‑GTAAAACGAGGGGCCAGGGGGCGAAA‑ATTttgttgttg                                                                                                                                                                  >  2:379417/1‑228 (MQ=255)                                                      aTGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑Tgaaga                                                                                                                                                                                                   >  1:463168/1‑196 (MQ=255)                                                        gAAGTTTTAT‑CCCCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGGCGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAACGCAGGTATTGCAGAAATCATGACCGCTCTGGCGGTATTGAC‑AC‑CAAAAAAAATTCCCCCTGGTGAAATTCGCGGTGCCTTTACCCCGGG‑TGAAGAAA‑TGGGGCACAGGGCGACACCtttttgtttttac                                                                                                                                                                >  2:187817/1‑219 (MQ=255)                                                                        gTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACAACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTAcc                                                                                                                                                                                                               <  1:513832/168‑1 (MQ=255)                                                                         ttATTTTCCCGGGGTCTCATCGGCTACTGGGTCAGACGCGGATTTCCACCGATGGGAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGc                                                                                                                                                 <  2:51005/227‑1 (MQ=255)                                                                                 ccGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTac                                                                                                                                        <  1:19353/228‑1 (MQ=255)                                                                                  gggggctgcATAGGTTACGGGGTAAGACGCTGATTACCACCGAAGGTAAAACCATGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCAGGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTaca                                                                                                                                       <  2:283409/224‑1 (MQ=255)                                                                                  cggttgtgcATTAGCTTCTGGGTTAGCAGCGGATTACCAACGATGGTAAAACCTTGTTAGGTGCCGAGGAAAAAGCAGGTATTGCAGAAATTATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTaca                                                                                                                                       <  2:129778/222‑1 (MQ=255)                                                                                  cGGTGCTGCATCAGCTACTGGGTCAGACGGTGATCACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTaca                                                                                                                                       <  1:147100/228‑1 (MQ=255)                                                                                      gctgcATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCGGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGGGATATTAGCGTCGCCTTTACCCCGGA‑TGACAAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTTACGCCTTCGATGCCCGCTGGGCTTACCCTGt                                                                                                                                   >  1:413544/1‑228 (MQ=255)                                                                                         tcATTAGGTACGGGGGCAGACCATGGTTACCACCGATGGTAAAACCTTGGTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCATTCGATGCCCGCTGGGCTTACACTGTTGa                                                                                                                                <  1:117288/227‑1 (MQ=255)                                                                                          cATTAGGTACTGGGTCAGACGCTGATTACCAACGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGa                                                                                                                                <  1:405983/227‑1 (MQ=255)                                                                                          cATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCg                                                                                                                                                  >  1:72850/1‑209 (MQ=255)                                                                                           aTCAGCTACTGGGTCAGACGCTTATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGc                                                                                                                                                 <  2:72850/209‑1 (MQ=255)                                                                                              aGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGAtggtg                                                                                                                           <  1:379417/228‑1 (MQ=255)                                                                                               gCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGCAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCTGAATGGGGATATTCGCGTCGCCTTTACACCGGA‑TGAAGAAC‑TGGGCAAAGGGCGGAAAAATTTTGATGTTGGCGCCTTCGAGGCCAGCATGGCTTCCCCTGtggtggggg                                                                                                                           >  2:502158/1‑219 (MQ=255)                                                                                                    tGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGGGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGACACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTAAACTGCTGATGGTGGTGGCg                                                                                                                     >  1:172484/1‑228 (MQ=255)                                                                                                                tGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATAGGGACATTCGCGCCGCCTTTCCCCCGGA‑AGAAGAAG‑TGGGCAAAGGGGCGAAACCTTTTGATGTGTACGCCTTCGATGCCCGCTGGGGTGCAACTGTGTTTGGTGGTGGGGGAGGCGAAGT‑gg                                                                                                        >  2:213484/1‑228 (MQ=255)                                                                                                                 gATTACCACCGATGGGAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAAGTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTacac                                                                                                                                      <  1:681172/198‑1 (MQ=255)                                                                                                                          cGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGaaa                                                                                               >  1:21168/1‑227 (MQ=255)                                                                                                                            aTGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGCGTTTTAAAACt                                                                                            >  1:358832/1‑228 (MQ=255)                                                                                                                            aTGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCCCTGGGGTTAAACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGGGTTTGAAAACt                                                                                            >  1:513011/1‑228 (MQ=255)                                                                                                                                   aaCCTTGGTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAAcgc                                                                                     <  1:9375/228‑1 (MQ=255)                                                                                                                                        tGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTac                                                                                                                                        >  2:551690/1‑173 (MQ=255)                                                                                                                                        tGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACACCGGA‑TGAAAAAG‑TGGGCAAAGGGGCGAAACCTTTTTATGGTGACGCCTTCGATGCCGGCTGGGCGTACACTGGTAATGGGGGTGGGGAAGGCGAACT‑GGGGTTTGAAAAATTCTAcacccggg                                                                                >  2:614335/1‑224 (MQ=255)                                                                                                                                        tGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGAATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAAAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCGTCGATGGCCGCTGGGCTTACACTGTTGATGGTGTGGGCGTAGCAGAACT‑GGCGTTTGAAAACTTCAAcgccgggg                                                                                >  2:102201/1‑225 (MQ=255)                                                                                                                                         gTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGGTTCACCAGTTTTGGGTGGTGGCGTAGGGGAACTGGGGGTTGGAAACCTTAAAcgcgcgtg                                                                                >  2:199077/1‑223 (MQ=255)                                                                                                                                         gTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTaca                                                                                                                                       <  1:551690/173‑1 (MQ=255)                                                                                                                                            gggTGCCGTTGACAAAGCAGTTATTGCAGAAATTATGCCCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATGCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGt                                                                            <  2:211759/227‑1 (MQ=255)                                                                                                                                            aGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGCCCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGt                                                                            <  1:421684/228‑1 (MQ=255)                                                                                                                                             ggTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATGCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGt                                                                            <  1:187817/227‑1 (MQ=255)                                                                                                                                              gTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGGAAAAAATTCCGCATGGTGATATTCCCGTCGCCTTTACCCCGGC‑TGAAGAAG‑TGGGCAAAGGGGCTAACCATTTTTATGTTGACGCCTTCGTGGCCCGGCGGGCGTACCCTGTTGATGGTGATGCGGTAGGCGAACTCGGCGTTT‑AAAACTTCAATGCCGAGTGGGGCa                                                                          >  2:288316/1‑228 (MQ=255)                                                                                                                                              gTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTAATGTTGACGCCTTCGATGCCCGCTGGGATTCCACTGTTGCTGGTGGTGGCGTAGGCGAACT‑GGAGTTTGCGAACTTACATGCCCCGCCGGTCa                                                                          >  1:331961/1‑228 (MQ=255)                                                                                                                                                         aaaGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAtt                                                               >  1:422029/1‑228 (MQ=255)                                                                                                                                                         aaaGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGGGATAATCGCGTGGCCTTTACCCCGGA‑TGAAGAAG‑CGGGCAAAGGGGCGAAACATTTTGATGTTGCCGCCTTCGATGTCCGCTGGGCTTACACTCTTAATGGTGTGGGCGAAGGCGACCT‑CGAGGTTGAAAACTTCAACACCGCTGCGGGAAAGATCaaaact                                                               >  2:272889/1‑226 (MQ=255)                                                                                                                                                           aGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AAAAAAAAATTCCGCATAGGGATATTCCCGTCGCCTTTACCCCGGA‑TGAAAAAG‑GGGGCCAAGGGGCGAACCATTTTGATGATGACGCCTTCGATGCCCGCGGGGCTTACACCGTGGCTGGTGGTGGCGTAGGTGGACT‑GGTGGTTGGCAACCTCACCTCCGGGTTGGGCCATATGAAAATTGt                                                             >  2:567527/1‑228 (MQ=255)                                                                                                                                                             cagtttttGCAGAAATCAGGACCGCGCGGGCGGTATTTGA‑AC‑AGAAAAAAATTCCGCATGGTGATAGTCGCGTCGCCTTTCCCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTc                                                            <  2:358832/221‑1 (MQ=255)                                                                                                                                                             cAGGTAATGCAGAAATCAGGACCGCGCGGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTc                                                            <  1:269942/227‑1 (MQ=255)                                                                                                                                                             aaggtttttgaGAAATCATGCCCGCGCTGGGGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCg                                                           <  2:579722/218‑1 (MQ=255)                                                                                                                                                                 taatggcgAAATCATGACCGGTCTGGCGGTATTGCA‑AC‑AGAAAAACATTCCGCATGGGTATATTCGCGTCGCCTTTGCCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAGCATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTaa                                                       <  1:284502/221‑1 (MQ=255)                                                                                                                                                                  agggacagaaTCAGGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTaa                                                       <  1:372674/219‑1 (MQ=255)                                                                                                                                                                   ttGCAGTAATCATGCCCACGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGGCTTTGCCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTaa                                                       <  1:502158/226‑1 (MQ=255)                                                                                                                                                                      cAGAAATCATGACCGGGCTGGCGGGATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATg                                                  <  1:102201/228‑1 (MQ=255)                                                                                                                                                                           aTCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCAt                                             >  1:188294/1‑228 (MQ=255)                                                                                                                                                                                 tccGGGCTGGGGGTATTGAA‑AC‑AGAAAAAAATTCCGGATGGTGATATGCGCGTCGTCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGc                                       <  2:188294/227‑1 (MQ=255)                                                                                                                                                                                     cgcTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCCAAGGGGCGCAACATTTTGTTGTTGACGCCTTCGATGCCCGCTGGGCTTACACCGTTGATGGTGGGGGCGTAGGCCAACT‑GGAGTTTGAAAACTTCACACCCGCGTCGGTCAATATCAAAATTGTGGGGAACAATGTTAATcctgtctctt                                   >  2:372674/1‑220 (MQ=255)                                                                                                                                                                                     cgcTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGCCGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGGGGCGTAGGTGAACT‑GGAGTTTGGAAACTTCAACGGCGCATCGGGCAATATCAAAATTGACGGTAACAATGTTCCTCCgggaccgg                                   >  1:647958/1‑223 (MQ=255)                                                                                                                                                                                      gcTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGGAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGGGGGGTAGGCGACCT‑GGAGTTTGAAAACTTCAAGCCCGCGTCGGTCAATATCCAAAATGTCGGGACAAATTTTCATCCGGGCAAggg                                  >  1:294251/1‑227 (MQ=255)                                                                                                                                                                                         gggggTATTGCA‑AC‑AGAAAAAAATTCCGGATGGGGAATTGCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGaa                               <  2:473135/225‑1 (MQ=255)                                                                                                                                                                                          gCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GCAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATACGTGCACGGCGcaa                              >  1:187311/1‑228 (MQ=255)                                                                                                                                                                                                ttGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTAAACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGCGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCACAATTGTCGGTAACAATGTTCATCCCGGCACGGCGCAAGGAgg                         >  1:225331/1‑226 (MQ=255)                                                                                                                                                                                                ttGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCCCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAACCATTTTGATGTTGCACCCTTCGATGCCCGCTGGGCTTACAATTTTTATGGTGGTGGCGTATGCGAACT‑GGAGTTTGAAAACTTTCACGCCGGGTCGGGCCATATCAAAATTTTCGGTAAAAATGTTCATCCGGGCAAggggcaaagaaga                        >  2:628425/1‑217 (MQ=255)                                                                                                                                                                                                    a‑aC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTg                        >  1:340952/1‑224 (MQ=255)                                                                                                                                                                                                                        attgTGATATTCGTGGTGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTgg  <  1:272889/225‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   TCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTAT‑CACCGGTTATGTTCCCGGTGCTGCATCAGCTACTGGGTCAGACGCTGATTACCACCGATGGTAAAACCTTGTTAGGTGCCGATGACAAAGCAGGTATTGCAGAAATCATGACCGCGCTGGCGGTATTGCA‑AC‑AGAAAAAAATTCCGCATGGTGATATTCGCGTCGCCTTTACCCCGGA‑TGAAGAAG‑TGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACT‑GGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGG  >  NC_000913/1186116‑1186555

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