Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,293,767	C→G	15.8%	A123G (GCG→GGG)	lpoA →	OM lipoprotein stimulator of MrcA transpeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,293,767	0	C	G	15.8%	113.3 / 2.6	65	A123G (GCG→GGG)	lpoA	OM lipoprotein stimulator of MrcA transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (10/42);  new base G (4/7);  total (15/50)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

AGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAATATTC  >  NC_000913/3293530‑3294000                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            agggttgtttggCCTTTTTTCTTCAGGGAGTGCAGCAAAGCTCGGAGGAATCCAGGATCAACTGGCAATTACTCGCCATTCGGGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGc                                                                                                                                                                                                                                             <  2:197467/227‑1 (MQ=255)        ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGGGAAAA‑CACTGCTGGCGGGAAGGATTAAAC‑TGGCGCCAGAAGAATTTTCTGGCG‑CGCAAAACATGCTTGGGGAAAtta                                                                                                                                                                                                                                     >  1:137933/1‑237 (MQ=255)             ccTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTAACGCAAGAACCGAAACATGCGCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGAGAAAACACGGCGC‑GAAAGATTTTGCTGTGG‑CGCCAAAAATGTTGCGGAAAAAA‑ACAccc                                                                                                                                                                                                                               >  2:721472/1‑238 (MQ=255)               tttttctttCAGCATATTCAGAAAAGCTTGTATGTTACAAGGATCACCTGGAAATTACTCGCCTTTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGc                                                                                                                                                                                                                             <  2:216468/232‑1 (MQ=255)                    tCTTCAGCAGATGCAGGAAAGCTCTGATGATACCAGGATCAACTGGCAATTCCCCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGAtt                                                                                                                                                                                                                        <  1:843895/239‑1 (MQ=255)                    tCTTCAGCAGATGCAGCAAAGCTTTGATGATCCCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGAtt                                                                                                                                                                                                                        <  1:8246/239‑1 (MQ=255)                    tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGc                                                                                                                                                                                                                             <  1:743126/234‑1 (MQ=255)                                     aaaGCTTTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGc                                                                                                                                                                                                        <  1:207008/237‑1 (MQ=255)                                       aGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGGGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAgc                                                                                                                                                                                                    <  1:233271/239‑1 (MQ=255)                                         ctctGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCCCCGAAAGATTTTGATGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACCAAACC‑ACCAAgccc                                                                                                                                                                                                  >  1:449753/1‑237 (MQ=255)                                                 atcccAGGACCAACTGGAATTGTCCCCCCATTCGGGCACGGGTAAAAGAGGGTAAAACCGGGCAGGGTGGTGAGGTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCCG‑CGTCGCGAGAAAA‑CCCTGCGTGCGGTAGAGATTAACC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCAAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTgg                                                                                                                                                                                          <  2:792803/236‑1 (MQ=255)                                                        gATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTAACGCAAGAACTGAACGATGCTCAG‑CGTTGCGGGAAAA‑CAATGCTGGCGGTAAAGATTTAAC‑CTGCGCCAAAAGAATTTTGTGGCG‑CGCCAAACTTTCTTGGGGAAATC‑CCACACCCCCGTTTTGAAAAAAAAC‑AACAAACGCGCTAATgggaagaaa                                                                                                                                                                                   >  2:224150/1‑232 (MQ=255)                                                            aaCTGGCAATTTCTCGCCATTCGTGCACTGGTGAAAGAAGGTAAACCCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAAt                                                                                                                                                                               <  1:712039/239‑1 (MQ=255)                                                             aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGCTGCTCAC‑CGTCGCGAGAAAA‑CAATGCTGGCGGGAGAGATTAAAA‑TAGGGCACAAAGATTTTGTTCGCGCCGCAAACATCGGTGG‑GCAAAAC‑ACACCCGCCCGTTTTGAAAAAAACCAAAAAAGCAGGTAAATGGCAGGaaaaaaa                                                                                                                                                                               >  2:805455/1‑238 (MQ=255)                                                             aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAAATGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAA‑CTGCGCCGAAAGATTTTGCTGGCG‑CGCAAAAATTTCTGGGGAAAAACAACCCCTCCCATTTTAAAAAAAAACC‑AACAAACGCGTAACGGGAGGGCaaaaaa                                                                                                                                                                               >  2:743806/1‑238 (MQ=255)                                                                  cAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATg                                                                                                                                                                          <  1:106803/238‑1 (MQ=255)                                                                        cgcGCCATTGGTGCACGGGTAAAAGAAGGTAAAGCCGGGCAGGGGGTTGAGTTGTTTAA‑CAAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CCCTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATgccagcca                                                                                                                                                                   <  1:872951/237‑1 (MQ=255)                                                                            ttaTTCGGGCACGGGGGAAAGAAGGAAAAAGCGGGGAGGCGGTTGAGTTGTTTAA‑CCACCTCGCGAAAAAACTGAAGGATGCCCGG‑CGTCCCGAAGAAA‑CCTTGCTGGCGGGAGAGTTTAAAC‑TGGCGGAGAAAGATTTTGCTGGCG‑CGAAAAACTTGCTGGGGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAgggg                                                                                                                                                               <  2:641815/237‑1 (MQ=255)                                                                            ccttgcGTGCCCTGGTAAAAGAAGGTAAACCCGGGCAGGCGGTTGGGTTGTTTAA‑CCAACTACCGCAAGAACTAAACGATGCCCGG‑GGTCGCGAAAAAA‑CCCTGCGGGGGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGAAAAACTTGCTGGGGAAAATC‑CCACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAgggg                                                                                                                                                               <  2:20834/234‑1 (MQ=255)                                                                                tCGTGCACTGGTGAAAGAAGGTAAAGCCGGGCAGGCGGTTGAGTTTTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTc                                                                                                                                                           <  1:863684/239‑1 (MQ=255)                                                                                tCGTGCAATGGTGAAAGAAGGTAAAACCGGGCAGGGGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTc                                                                                                                                                           <  1:321545/239‑1 (MQ=255)                                                                                    gCATTGGTGAAAGAAGGTAAAACCGGGAAGGCGGTTGAGTTTTTTAA‑CCAACTACCGGAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTc                                                                                                                                                       <  1:843954/239‑1 (MQ=255)                                                                                      tgtGGTGAAAAAAGGTAAAAGGGGGCAGGGGGTTTGGTTTTTTAA‑CCATCTACCGCAAGACCTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGAGTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  2:253996/237‑1 (MQ=255)                                                                                      ccgggTGAAAAAAGGTAAACCCGGGCAGGCGTTTGGGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTAAG‑CGTCGCGAGAAAA‑CACGGCTGGAGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGGGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  2:300176/236‑1 (MQ=255)                                                                                       cTGGTGAAAGGAGGTAAAACCGGGCCGGGGGGTGGGTTGTTTAA‑CCAACTACCGGAAGAATTGAAGGATGCTCAG‑CGTAGCGAGAAAA‑CAGTGCTGGGGGTAGAGATTAAAC‑TGGCGCAGAAATATTTTGCGGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  1:737096/238‑1 (MQ=255)                                                                                          gTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAATTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGTTGGCGGTGGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAAc                                                                                                                                                                                                             <  1:511995/180‑1 (MQ=255)                                                                                                     ttaaaCGGGGCAGGGGGTTGAGTTTTTTAA‑CCAACTCCCGAAAGAACTGAACGGTGCCACG‑TGTCCGGAAAAAA‑CGCTGCGGGCGGAAGAGATTAAAC‑GGGCGCAGAAAGATTTTGCTGGCG‑CGAAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  2:648640/235‑1 (MQ=255)                                                                                                     caaaaCCGGGCCGGCGGTTGGGTTTTTTAA‑CCACATACCGCAAAAACGGAACGATGC‑CAGAGGTCCCGAGAAAA‑CCCGGCTGGCGGTAGAGATTAAAC‑TGGCGCAAAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  2:777828/236‑1 (MQ=255)                                                                                                     aaaaaCCGGGGAGGGGGTTG‑GTTTTTTAACCCAACTGCCGGAAGAAATGAACGGTGCCCAG‑GGCTGGAAGAAAA‑CCCTGCTGGCGGTAGGATTTAAAC‑TGGGGCAGAAAGATTTTGCTGGCG‑AGAAAAACTTGCTGGGGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  2:748880/236‑1 (MQ=255)                                                                                                       aaaCCGGGCAGGCGGTTGAGTTGGTTAA‑CCAACTCCCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CCCTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcgcgc                                                                                                                                     <  1:774167/238‑1 (MQ=255)                                                                                                            gggCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGACCTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTAGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCGGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAAt                                                                                                                               <  1:566816/239‑1 (MQ=255)                                                                                                                                      aaCTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        >  2:716029/1‑205 (MQ=255)                                                                                                                                      aaCTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  1:716029/205‑1 (MQ=255)                                                                                                                                        cccccGCAAGAGATGAACATGGCTCGG‑GGGCGCGAGAAAA‑CCCCGCGGGCGGTAGAGTTAAAAC‑TGCCGCAGAAAGATTTTGCGGGCG‑CGCAAAACTTGCTGGAGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCgaa                                                                                                    <  2:75633/236‑1 (MQ=255)                                                                                                                                                aaaaCTGAACGATGCCCAG‑CGTCGGGAGAAAA‑CCTGGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCGGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGaaaaa                                                                                             <  1:447540/236‑1 (MQ=255)                                                                                                                                                 aaaCTGAACGTGCCTCAG‑CGT‑GCGAAAAAACCCCGGCGGGCGGTAAAGATTAACC‑TGGCGCAAAAAGTTTTTGCTGGCC‑CGCAAAACTTGCTGGCGAAAAGC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAAca                                                                                           <  2:89338/238‑1 (MQ=255)                                                                                                                                                                ccaggggcgagagaAAC‑CCTTGCGGGCGGTAGAGATTAAAC‑TGGCGCAAAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGAt                                                                              <  1:886240/228‑1 (MQ=255)                                                                                                                                                                cAG‑CGTCGCGAGAAAA‑CCCTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGAAAACCTGGCTGGGGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGc                                                                            <  1:89285/239‑1 (MQ=255)                                                                                                                                                                    ggTCGCGAGAAAA‑CCCGGCGGGCGGTAGAGATTAAAC‑TGGCGGAGAAAGATTTTGCGGGCG‑CGAAAACCTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  2:584370/162‑1 (MQ=255)                                                                                                                                                                       cgcgAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAAATGCTCAGGAACCGCTGCTTGGCGCGCAAGAAAAAAAGAAGAatt                                                                                   >  1:532991/1‑225 (MQ=255)                                                                                                                                                                            gaAAA‑CCCTGCTGGCGGTAGAGTATAAAC‑TGGCCCAGAAAGATTTGTCTGGCG‑CCCAAAACTTGCTGGCGAAAACC‑CCACCTGCCGATTTAGAGCAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAg                                                                  <  1:781085/238‑1 (MQ=255)                                                                                                                                                                            gaAAA‑CCCTGCGGGCGTGTGATATTAAAC‑TGGCGAAGAAAGATTTTGCTGGCG‑CGAAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAg                                                                  <  2:137933/238‑1 (MQ=255)                                                                                                                                                                              aaa‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGGAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACt                                                                                                                                           >  1:584370/1‑163 (MQ=255)                                                                                                                                                                               aa‑CCCTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGgcgc                                                              <  1:763244/239‑1 (MQ=255)                                                                                                                                                                                a‑cacTGCTGGCGGTAGAGATTAAAC‑TGGCGAAGAAAGATTTTGCCGGCC‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGct                                                             <  1:707558/239‑1 (MQ=255)                                                                                                                                                                                  ccctgctgGCGGTAGAAATTAAAC‑TGGCGCAGAAAGATTTGTCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc                                                            <  1:857912/237‑1 (MQ=255)                                                                                                                                                                                     tgctgGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCCAAAATCGATGCCACCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCCACGGGACCGCGCCTTGTGCCGGACAAAAAACCACAGAAAATTGTTGCCACCCGGccagggccccc                                                          >  2:309136/1‑228 (MQ=255)                                                                                                                                                                                       ctgGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACCGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGGAAGAAAAACAGCAAAATATTGATGCCACCTGGCAGGCGCTcccct                                                        >  1:266227/1‑238 (MQ=255)                                                                                                                                                                                        tgGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCGGCTTGGCGCGAAAGAAAAACAGCCAGATATTGATGCCACCTGGCAGGCGCCCtccccc                                                      >  1:359972/1‑239 (MQ=255)                                                                                                                                                                                          gggttAGAGATTAAAC‑TGGCGAAAAAAGTTTTTGCTGGCG‑CGAAAAATTTGCTGGGGAAAACA‑CCCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  2:266227/235‑1 (MQ=255)                                                                                                                                                                                           ggttAGAGATTAAAC‑TGGCGAAGAAAGATTTTGCTGGCG‑CGCAAAGCTTGCTGGCGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  1:636199/235‑1 (MQ=255)                                                                                                                                                                                           cGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  1:116256/238‑1 (MQ=255)                                                                                                                                                                                             gTAGAGTTTAAAC‑GGGCGCAGAAAGATTTGGCTGGCC‑CGCAAACCTTCATGGCGAAAATC‑ACCCCTGCCGATTTAAAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGa                                                  <  1:876308/238‑1 (MQ=255)                                                                                                                                                                                              aagagATTAAAC‑TGGCGAAGAAAGATTTTGCTGGCG‑GGAAAAACTTGCGGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGAc                                                 <  1:506773/237‑1 (MQ=255)                                                                                                                                                                                                    ttAAAC‑GGGCGCAGAAAGATTTGGCTGGCG‑CGAAAAACTTGCTGGCCAAAATC‑ACCACGTCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAgg                                            <  1:825350/237‑1 (MQ=255)                                                                                                                                                                                                         acgggagcAAAAAGA‑TTTGCTGGAG‑AGAAAACCTTGCTGGCGAAAATC‑ACCCATGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGc                                     <  1:9137/233‑1 (MQ=255)                                                                                                                                                                                                              gccggaaaacttttGCTGGAG‑CACAAACCGTGCTGGCAAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCTTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGaa                                  <  1:186299/228‑1 (MQ=255)                                                                                                                                                                                                              cgcAAAAAGATTTTGCTGGCG‑CGAAAAACTTGCGGGCAAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGaa                                  <  1:644844/238‑1 (MQ=255)                                                                                                                                                                                                                  gAAAGTTTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACg                              <  1:594326/238‑1 (MQ=255)                                                                                                                                                                                                                           tGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGAATCAGGAACAGGCGACTAAGCTGGGGTTCa                    >  1:831162/1‑239 (MQ=255)                                                                                                                                                                                                                           tGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGAAGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCCGGGGATCa                    >  1:14414/1‑239 (MQ=255)                                                                                                                                                                                                                              gggCG‑CCGAAAACTTGCTGGCAAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACg                 <  1:721472/238‑1 (MQ=255)                                                                                                                                                                                                                                 cg‑cgcAAAACTTGCTGGCGAAAATC‑CCCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGcc               <  1:743806/238‑1 (MQ=255)                                                                                                                                                                                                                                 cg‑cgaAAATCTTGCGGGGAAAAACC‑CCCCCTGCCGATTTAAAAAAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAAGCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGcc               <  2:843835/238‑1 (MQ=255)                                                                                                                                                                                                                                           cTTGCTGGCGAAAATC‑CCCCCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAata     <  1:224150/239‑1 (MQ=255)                                                                                                                                                                                                                                              gCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTTATGCCACCTGGCAGGCGCTCCCCCCCATGACTCAGGAAAAGGGGGAAAACCGGGGGATAAAAGCCGACCaaaaaaatt  >  2:644498/1‑234 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            AGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAA‑CCAACTACCGCAAGAACTGAACGATGCTCAG‑CGTCGCGAGAAAA‑CACTGCTGGCGGTAGAGATTAAAC‑TGGCGCAGAAAGATTTTGCTGGCG‑CGCAAAACTTGCTGGCGAAAATC‑ACACCTGCCGATTTAGAACAAAACC‑AGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAATATTC  >  NC_000913/3293530‑3294000

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