Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,762,307	A→G	29.4%	G42G (GGA→GGG)	rhsA →	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,762,307	0	A	G	29.4%	55.9 / 4.0	43	G42G (GGA→GGG)	rhsA	Rhs protein with putative toxin 55 domain, putative polysaccharide synthesis/export protein, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (9/19);  new base G (12/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 2.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  NC_000913/3762076‑3762495                                                                                                                                                                                                                                            |                                                                                                                                                                                                ttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTTAACCCGGGGGCGCCTTTGTGCCCCCC‑CCGCCGTGGCCTGTTTGGT‑GT‑GCCCCGGCGGAggggag                                                                                                                                                                                           >  1:121792/1‑233 (MQ=17)           ttttcctttGATGGATTTCTCACTTCAGGCTTTCCTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATcc                                                                                                                                                                               <  2:691287/231‑1 (MQ=255)                   tGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCC‑CCGGCGGTGCCGGTTCGGG‑GG‑GCCCCCGCGGGGGGACCGCCG‑GCCACCCCGCCAAt                                                                                                                                                                        >  1:387351/1‑238 (MQ=14)                   tGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAATTTTCAGGGTTCCACCCGGGGGCGCCTTGCTGCCCCCA‑CCCGCGTGGCCTCTTCTGG‑GA‑GCCCCGGCGGGGGGGCGTccccgccccccggcaac                                                                                                                                                                         >  2:615687/1‑223 (MQ=14)                          tatCTCACTTCAGGCTTTCTTATAAATCTGGAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCt                                                                                                                                                                <  1:355663/239‑1 (MQ=255)                           atCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCt                                                                                                                                                                <  1:326087/238‑1 (MQ=255)                             ctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGGGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGG‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATacccctcccgg                                                                                                                                                             >  1:692636/1‑232 (MQ=255)                               cACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GTGGCCCCGCCGGAG‑GACGTCCA‑GCCATCCGGTCAATCCCCTGCG‑CGGt                                                                                                                                                           >  1:1043958/1‑238 (MQ=255)                               cACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCA‑CCGGCGTGGCCCGTTCGGG‑GT‑GGCCCGGCGGGGGTGCGGCCG‑GCCCCCCCGCGCATCCCCCGCT‑Cggg                                                                                                                                                           >  1:748821/1‑237 (MQ=14)                                 cTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCa                                                                                                                                                        <  1:485301/239‑1 (MQ=255)                                  ttCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGGCGGGCCTGTTCGGG‑GG‑GACCCGGCGGGGGGAGGTCCG‑GCCATCCGGGCAAACCCCTCCT‑CGGTGcca                                                                                                                                                       >  1:767946/1‑237 (MQ=14)                                  ttCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGGCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCa‑a                                                                                                                                                      <  1:73859/239‑1 (MQ=255)                                   tCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCA‑CCGGCGGGGCCTGTTCGGG‑GT‑GCCCCGGCGGGGGGACGGCCG‑GCCCTCCGGTCAATCCCCCGCT‑CCGGGCa‑aa                                                                                                                                                     >  1:264740/1‑239 (MQ=14)                                     aGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGGGGCCCCC‑CCGGCGTGGCCTGTTCGGG‑GG‑GCCCCGGGGGGGTGCGGCGCC‑GCCCACCGGCCAATCCCCCGCC‑CCGGGCA‑AAgg                                                                                                                                                   >  2:326087/1‑239 (MQ=14)                                       gCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCCTTGGGGCCCCCA‑CCGGGGGGGCCCGTTCGTG‑GG‑GGCCCGGCGGGGGGGCGTCCC‑GCCCTCCGGTCCATACCCCGccggggggcaagggc                                                                                                                                                   >  2:485654/1‑225 (MQ=11)                                                taaagcgggtgtgTTTCGCCTGTCAGCAGACAAAAAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGt                                                                                                                                        <  2:121792/227‑1 (MQ=255)                                                     cTGGAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAAc                                                                                                                                   <  1:701206/239‑1 (MQ=255)                                                       gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGGGGGACGTCCG‑GGCATCCGGCCCCTCCCCTGCC‑CGGGGCA‑AAGGGCTCTTCCGGGGAAAcca                                                                                                                                 >  1:576268/1‑238 (MQ=14)                                                       gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGAAGCAATGTTCAGGGGTCAGCCCGGGGGCCCCCTGGGGCGCCCC‑CCGGCGCTGGCTGCTCCGTGGG‑GCCCCGGCGG‑GGGACCGCCC‑GCCCCCCGGCGCATCCCCCCCC‑CGGCGCA‑CAGGGCCTTccccgggacaac                                                                                                                                  >  2:162423/1‑229 (MQ=9)                                                        tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAATGTCCAGGGGTCAGCCGGGGGACGCATCGGTGGCCCCC‑CCGGCGTGGCCCGTTCCGT‑GT‑GCCCCCGCGGGGGGGCGGCCC‑GGCCCCCGGGCCAACCCCCGCC‑CGGGGCACAAGGGCC‑TCCCGGGGAAAca                                                                                                                                  >  2:928773/1‑237 (MQ=14)                                                             tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGGTCAGCCGGGGGACGCACCGGGGCCCCCA‑CCGGCGGGGCCCGGTCGGG‑GT‑GGCCCCGCGGGGGGGCGCCCG‑GCCCCCCCGGCAAACCCC‑GCTCCGGGGGA‑AAGGGCCCTCCCCGGGACAcccccac                                                                                                                             >  2:170949/1‑233 (MQ=9)                                                                cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGAAGCATTGTTCACGGGTCAGCCGGGGGGGGCATTTGTGGCCCCACCCGGCGGGGCCGGTT‑GGG‑GT‑GCCCCGGCGGGGGGGCGTCCGCGCCACCCGGGCCTCCCCCGCCC‑CGGGGCA‑AGGGCCCTCCCCGG‑GAAAccacaaccacc                                                                                                                         >  2:442723/1‑229 (MQ=9)                                                                             ggagatgggtataagatgcggggATGAG‑AAATGGGCGGAAAACCGGCACCGCGGCAGGGCGCAATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCt                                                                                                           <  1:151287/218‑1 (MQ=14)                                                                                    aCCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCCGGGTTCAGCCGGGGGGCGCATTTGGTCCCCCA‑CCGGCGTGGCCCGTTCGGT‑GT‑GCCCCGGCGGGGTGGCGTCCG‑GCCATCCGGGCAAATCCCTCCC‑CGGGGCA‑AAGGGCCCTCCCGGGGAAACCCCCATCcccccccccccccccccccacgtc                                                                                                    >  2:48929/1‑219 (MQ=14)                                                                                             aaCAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc                                                                                           <  1:111571/239‑1 (MQ=255)                                                                                                               aGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGcc                                                                                                                         <  1:640428/192‑1 (MQ=37)                                                                                                                 cGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGG‑AGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATc                                                                                                                                                                                >  2:248198/1‑136 (MQ=255)                                                                                                                 cGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGG‑AGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATc                                                                                                                                                                                <  1:248198/136‑1 (MQ=255)                                                                                                                       aCCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGt                                                                                                      <  1:803976/203‑1 (MQ=35)                                                                                                                                 cgTCAGGGCGAAATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGGGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACGTCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcg                                                         <  1:624128/238‑1 (MQ=18)                                                                                                                                        gcgcctgggcgcgggttGGCGGTGGCATGTTTCGGGGTTCAGCCGGGGGGCGCGTTGGTGCCCCCG‑CCGGCGTGGCCTGTTCGGT‑GG‑TCCCCGGCGGGGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGt                                                  <  2:576268/222‑1 (MQ=11)                                                                                                                                            cATGGCCCAGTATGGCGGTAGCATTGGTCAGGGTTCCGCCGGGGTGCCCATTGGTTCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTgggg                                              <  1:485654/238‑1 (MQ=255)                                                                                                                                                                ccgtttttCGGGGTTCGCCGGGGGGGCGCATTGGGGCCCCCC‑CCGGCGTGGCCTGTGCGGT‑TG‑GCCCCGGCGGAGTGACGGCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaaa                         <  2:692636/233‑1 (MQ=11)                                                                                                                                                                    tGTTCAGGGTTCCGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGc                     <  1:738534/239‑1 (MQ=17)                                                                                                                                                                              tcggCCGGGGTGCGCATTGGTGCCCCCC‑CCGGGGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGCCCG‑GCCATCCGGTCAACCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat           <  2:264740/236‑1 (MQ=12)                                                                                                                                                                              tccggcGGGGTGCGCATTGGGCCCCCCA‑CCGGGTGGGCCTTTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat           <  2:767946/234‑1 (MQ=255)                                                                                                                                                                              cccgCCGGGGGGCGCTTGGGTGCCCCCC‑CGGGCGTGGCCTTTTCGGT‑GT‑GCCCCGGCGGGGTGCCGCCCG‑GCCACCCGGTCAATCCCCTGCC‑CGGTGCA‑AAGGTCCTTCCCGGTGAACCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat           <  2:748821/236‑1 (MQ=11)                                                                                                                                                                                aGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCCCCTACAGCAGTTACCGGCCAAAAACCCCCCCGCCCGTGGGGGGCCCCGGCCCCGGCCGGAAAAGGCCCGCGGAAATcg         >  2:458521/1‑238 (MQ=11)                                                                                                                                                                                aGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACCAAAACGCCCGCGCCGGGGGGGGGCCTCCGCCCCGGCTGGAAAAAGCCTGCGGATATcc         >  2:991821/1‑239 (MQ=12)                                                                                                                                                                                aGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTCGGGAGCCTCGGCCCCGGCCGGAAAATGCCTCCGGATATc          >  1:73198/1‑238 (MQ=17)                                                                                                                                                                                aGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCGCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGCGGGGGCGCCCGGCCCCGGCAGGGAAAAGCCAGCGGAAATcc         >  2:278835/1‑239 (MQ=11)                                                                                                                                                                                   cGGGGTGCGGATTTGTGCCCCCA‑CCGGCGTGGCCTGTTGGGT‑GT‑GCCCGGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCt      <  1:615687/239‑1 (MQ=9)                                                                                                                                                                                       gTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGGGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTCCTCTCCCGCACCTACAGCAGTTCCCGGACCACAACGCCCCCCCCGGTGGGGGGCCCCGGCCCCCGCTGGGAAACGCCTGCGGGTAGCCCCTTaaa  >  2:546083/1‑237 (MQ=11)                                                                                                                                                                                                                                            |                                                                                                                                                                                                TTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCA‑CCGGCGTGGCCTGTTCGGT‑GT‑GCCCCGGCGGAGTGACGTCCG‑GCCATCCGGTCAATCCCCTGCT‑CGGTGCA‑AAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACA  >  NC_000913/3762076‑3762495

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