Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	352,834	C→G	16.1%	A302G (GCG→GGG)	prpE →	propionate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	352,834	0	C	G	16.1%	172.3 / 2.5	70	A302G (GCG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (15/42);  new base G (6/6);  total (21/49)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCG  >  NC_007779/352597‑353072                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 gggTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGc                                                                                                                                                                                                                                                  <  1:604111/238‑1 (MQ=255)   gTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGGGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGa                                                                                                                                                                                                                                                <  2:1263894/238‑1 (MQ=255)     aCCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAGCCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTa                                                                                                                                                                                                                                             <  1:988512/239‑1 (MQ=255)                        tCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGcc                                                                                                                                                                                                                          <  2:1689272/239‑1 (MQ=255)                          cAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGa                                                                                                                                                                                                                        <  1:1301322/239‑1 (MQ=255)                           aaCGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGAc                                                                                                                                                                                                                       <  2:1789097/239‑1 (MQ=255)                            aCGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGAcc                                                                                                                                                                                                                      <  1:1480198/239‑1 (MQ=255)                                aaCCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGTTATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTGCGCGCCGCGGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATTGCCGAGCTgg                                                                                                                                                                                                                   <  2:1298015/238‑1 (MQ=255)                                    tcctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGa                                                                                                                                                                                                              <  1:1263057/239‑1 (MQ=255)                                      ctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACt                                                                                                                                                                                                            <  1:1681094/239‑1 (MQ=255)                                      ctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACt                                                                                                                                                                                                            <  2:837783/239‑1 (MQ=255)                                      ctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACt                                                                                                                                                                                                            <  1:1564220/239‑1 (MQ=255)                                        gCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGGGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGGGACTATCGTTTACGAAAGGTTGCCCACCCGGCCGGGCTgc                                                                                                                                                                                                          >  1:302721/1‑239 (MQ=255)                                          aTTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcg                                                                                                                                                                                                         <  1:1735883/238‑1 (MQ=255)                                          aTTCTCGCCCCTGCCGGCACGCCCGGCAAACCTAAAGGCGGGCAACGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACTCCATTTTTGGCGGCAAAGCGAGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcg                                                                                                                                                                                                         <  2:2020807/238‑1 (MQ=255)                                           ttCTCTACACTTCCGGCACGACCGGCAACCCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                                       <  1:987380/239‑1 (MQ=255)                                           ttCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTACTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                                       <  1:1623036/239‑1 (MQ=255)                                             ctctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGGCCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAATGATTGCCGACCTGGCCGGACTGCGGCgt                                                                                                                                                                                                     <  1:436317/239‑1 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                                    <  1:369785/239‑1 (MQ=255)                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTAACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                                    <  2:576250/239‑1 (MQ=255)                                                     ttCCGGCCAGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGAGATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCCGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGa                                                                                                                                                                                            <  2:2104198/239‑1 (MQ=255)                                                     ttCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGCTGGCGGCTATCGTTTTTGAAGGATTGCCGCCCCGGCCCGGCTGCGGCGCGt‑ggggg                                                                                                                                                                                             >  1:1556285/1‑238 (MQ=255)                                                       ccGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCTATCGGTTACGAAGGATTTCCGACCCTGCCGGGCTGCGGCGTgg‑gggtggaa                                                                                                                                                                                          >  1:216549/1‑233 (MQ=255)                                                                ccGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCgag                                                                                                                                                                                 <  1:1734391/239‑1 (MQ=255)                                                                        ccTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGAc                                                                                                                                                                                                                       <  1:2188002/194‑1 (MQ=255)                                                                          gaaaGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAgg                                                                                                                                                                       <  1:507572/238‑1 (MQ=255)                                                                           aaaGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGt                                                                                                                                                                      <  1:1362307/239‑1 (MQ=255)                                                                             aGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGGCCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTa                                                                                                                                                                    <  1:1165603/239‑1 (MQ=255)                                                                                   gCAGCGTGGGGGTGGCGTATATGCGGGGGGGCTGGCGACCTGTATGGACACCATGTGTGGCGGCATAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCAGTCGTGTATCGTTTACGCGCCGCTGCTGGCGAGGATGGTGACTATCGTTTGCGATGGATTGTCGACTTGTGCGTGCTGCGGCCTGT‑GGTGGTCAATCGTCGAGAAATATCATGTAAGCCgg                                                                                                                                                               <  2:871976/238‑1 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGGCTATCGGTTAAGAATGGTTGCCGCCCTGGCCGGACTGCGGGGGGT‑GGGGGACAATCGTCGGGGAAAATCAGGGTAGGCggg                                                                                                                                                              >  1:1797198/1‑238 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTCTCGTTTTCGAAGGATTGTCCACCCGGGCGGGCCGGGGGGTGG‑GGGGGGCCAACGGCGGGAAAAATCAGGGTAGcccgg                                                                                                                                                              >  1:8377/1‑236 (MQ=255)                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTTCGAAGGATTGCCGACCTTGCCGGGCCGCGGCGTGT‑GGGGGGCAATCGTCGAGAGATATCTGGGTAGACggg                                                                                                                                                              >  1:985966/1‑238 (MQ=255)                                                                                     aGCGGGGCGTCGGCGGATATGCGGGGGCGCGGGCGCCCTCGATGGAAACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCAGCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGAt                                                                                                                                                             <  2:1288676/238‑1 (MQ=255)                                                                                     aGCGGGAGGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGGGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGAt                                                                                                                                                             <  2:604186/238‑1 (MQ=255)                                                                                           aCGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCTATCGTTTTCGAAGGATTGCCGACCTGGCCGGACTGCGGGGTGT‑GGGGGGCAATCGTCGAGAAATATCAGGGTAGCCGGGTGGtctcc                                                                                                                                                      >  1:1356218/1‑238 (MQ=255)                                                                                                     atatGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCgag                                                                                                                                                                                 <  1:547280/202‑1 (MQ=255)                                                                                                                  tgggCGACCTCGGTGGCCACCATTTTGGGGGGCAAAGCGGGCGGCGTGTTCTTTTGCGCTTGGGATATCGGCTGGTTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGCCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGc                                                                                                                                <  2:1797198/236‑1 (MQ=255)                                                                                                                   ggggcaCCTCGAGGGAGACCATTGTTGGCGGAAGGGCGGGCGGGGTGTTGTTTGGCGCAGCGGGTGTCGGGGGGGGGGTGGGGGTTTCGTATTGCGGTTCGGCCCCGGGGGTGGGGGGGTGGGCGTCTATCGTTTACGAAGGATTGCTGACCTGGCCGGGCTGCGGCGTGT‑GGTGGACAATCGACGAGAAATATCAGTTTAGCCGGATGTCCTCAGCGCCGACCGCCATTCGCGTGCt                                                                                                                               <  2:801695/233‑1 (MQ=255)                                                                                                                          cTCGATGGACGCCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCACCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCCCAGCGCC‑ACCGcatatccgc                                                                                                                                   <  1:1021619/226‑9 (MQ=255)                                                                                                                                        tttttGGCGGCAAAGCGGGCAGCGTTTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATTGCCGCCCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCGGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGaa                                                                                                         <  1:2049016/239‑1 (MQ=255)                                                                                                                                         ttttGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGGGGCTTTCGGTTTCCCAGGATTTCCCCCCTGGCCGGACGGCGGCGGGG‑GGGGGGTAATCTCGGCGGAATATCTAGGTAGCCGGGGGTTTTCGGCGGCCCCCGCCCCCTGCCGGCTGAAAAAAATCTCTACCCCTCaac                                                                                                        >  2:2007835/1‑238 (MQ=255)                                                                                                                                                       ggggCGGGGTGTTTTTTGGCGATTGGGATATGGGCTGGGTGGTGGGGCATTCTTATATCTTTTGCGCGCCGCGGCTGGCGGGGGTGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACg                                                                                           <  2:2086354/236‑1 (MQ=255)                                                                                                                                                              cGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTTCTGGCGGGGGGGGGGGCTTTCGTTTACGAAGGATTGCCGACCCTGCCGGGCGGGGGGGGGTGGGGGGACAA‑CGGCGGGGAAAATCAGGGTAGCCGGGGGGTTTCCGCGCCCACCCCCATTCCCGGGGCGAAAAAAATCCCTCCCCCTGTAATTCTCCAAacacacctctcc                                                                                    >  2:1123691/1‑230 (MQ=255)                                                                                                                                                                 gtTCTTTTGCGCATCGGATATCGGCGGGGTGGTGGGGCTTTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGc                                                                               <  1:1123691/239‑1 (MQ=255)                                                                                                                                                                  ttCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGc                                                                               <  1:2173780/238‑1 (MQ=255)                                                                                                                                                                      tttGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAGGGGGACTATCGTTTACGAAGGATTGCCGACCGGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGaa                                                                          <  1:13881/239‑1 (MQ=255)                                                                                                                                                                       ttGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAg                                                                         <  1:1349284/239‑1 (MQ=255)                                                                                                                                                                          gccATCGGATGTCGGCTGGGTGGGGGGGCTTTGGTGTTTGGGTTCCGGGCGGCGGCTGGCGGGGATGGGGACTATGGTTTACGAAGGATTGCCGACCGGGCGGGACTGCGGCGTGT‑GGTGGACAACCGTCGAGAAATATCAGGTTAGCCGGATGTTCCCGGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTg                                                                       <  2:2062394/236‑1 (MQ=255)                                                                                                                                                                            gggcggtgttCGGGTGGGTGGGGGGGTATTCGTATGTCGTTTACGCGCCGCGGCTGGGGGGGTTGGCGACTATCGTTTACGAGGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATAGCAGGTTGGCCGGATTTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTTCCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGct                                                                     <  2:964691/229‑1 (MQ=255)                                                                                                                                                                                  atgtCGGGTGGGTGGTGGGGCTTTAGTATATCGTTGCCCCGCCGCTGCTGGGGGGGAGGGGGCCTATCGGTTACGAGGGATTGCCGCCCTGGCGGGACGGCGGGGTGT‑GGTGGACAACCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATct                                                               <  2:400200/235‑1 (MQ=255)                                                                                                                                                                                   tagtggggggggggTGGGGCATGTGTTTAGCGTTTCCGCGCCGCTGCTGGGGGGGAGGGGGACATTCGTTTACGAAGGATTGCCGACCTGGCGGGACTGCGGGGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctg                                                              <  2:1356218/230‑1 (MQ=255)                                                                                                                                                                                          tGGGTGGTGGGGCATTCTTATATCGTTTGCGCGCCGTTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                                       <  1:1221116/238‑1 (MQ=255)                                                                                                                                                                                                      cATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGGTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGGGGACAATCGTCCGGAAATATCAGGGTAGCCGGATGTTCTCAGCCCCGACCGCCCTTCGCGTGCTGAAAAAATTCCCTCCCGCTGCAATTC‑GCAAACACCATCTCTCGTCCCTGGACGTGCTCTCTCTCGCTGTGGGACCGcctggcgc                                          >  2:137626/1‑232 (MQ=255)                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTACCAAAGATTGCCGACCTGGCCGGGCTGCGGCGGGT‑GGTGGTCAATCTGCGGGAAATATCAAGTTATCCGGGTGTTCTCAGCCCCGACCGCCCTTCGCGTGCTGAAAAAATTCCCCACCCCTGGAATTC‑GCCAACATGCACTCTCCTCGCCGGAAGAGCTCTCTCTCGCCGGCGAACCGCTGGCCggg                                         >  2:2208364/1‑236 (MQ=255)                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGAATATCGTTTACGAAGGATTGCCGACCTGGCCCGACTGGGGCGGGG‑GGGGGACAACAGTCGGGGAAAAACAGGGTAGCCGGATGTTTTCAGCGCCGACCGCCCTTGGCGTGCGGAAAAAAATTCCCACCGCTGAAATTC‑GCCAACAAGGACTCTCGCCCCCGTGACTGCTCTCTCTGGCTGGGGAACCGCTGGagcca                                         >  2:865776/1‑234 (MQ=255)                                                                                                                                                                                                         tcttataTGGTTTGCGGGCCGCGGCTGGGGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGGGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTg                                                                       <  2:1688262/204‑1 (MQ=255)                                                                                                                                                                                                           gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAgg                                                                                                                                                                       <  1:1572500/110‑1 (MQ=255)                                                                                                                                                                                                             atatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGGcc                                   >  1:310305/1‑239 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATTGTTTACTGAGGATTGCCGACCTGGCCGGACCTCCGCGCGG‑GGTGGGCAATACTCGGGGAATATCAGGTTAGCCCGGTGTTCTCAGCGGCCACCCCCCCTTCCGTGCCGAAAAAATTTCCGACCCCCGAAATTC‑TCAAACACGATCCTCCGTCGCTGGAACTGCCCTATCCGGTCCGAGAACCGCCGCACGCGCCCAcgaccact                             >  2:1330051/1‑232 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGAt                                                                                                                                                             <  2:737418/112‑1 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAACCGGAt                                                                                                                                                             >  1:737418/1‑112 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctgg                                                             >  1:1688262/1‑207 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGCGT‑GGTGGGCAATCGTCGAGGAAAATCAGGTTAGCCGGGTGTTTTCAGGGCCGAACGCCCTTCGCGGGCTGAAAAAATTCCCTACCGCCGAAATTC‑GCAAACACCCTCTCTCGTCGCTGGAAATGCCCTCTCTGGCGGGAGAACCGCTGGACGAGCCGACCCCCCGt                             >  2:1126783/1‑239 (MQ=255)                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAAGAATGCCGACCCGGGCGGGGTGCGGGGGGG‑GGGGGACAATCCGCGAGAAATATTCGGGTAGCCGGGTGGTCTCCGCGCCGGCCGCCCTTCGCGGGCTGGAAAAATTCCCCAACGCTGGAATTC‑GCAAAAACGAACCCTCGCGCCTGGAAAAGCTCTTTCTGGCCGGGAAACCGCTGGGGGGGCCGACCGcccgg                             >  2:83533/1‑236 (MQ=255)                                                                                                                                                                                                                     tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGAGTGCCGACCTTGCCGGGCTGCGGCGGGG‑GGGGGGCAAACGGCGGGAAATATCCGGGTAGCCGGAAGGTTTCCGCGCCGGACCCCCTTCGCCGGGCGGAAAAATTCCCCACCGCCGGAAATA‑CCAAACACGATCCCCCGTCGCTGGGAAGGCCCCATTTGGCCGGGGAACGCCCGGGCGAGCCcgccgcaccgtg                           >  2:1822191/1‑233 (MQ=255)                                                                                                                                                                                                                      ccgcgCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACGGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg                          <  1:1328580/238‑1 (MQ=255)                                                                                                                                                                                                                           ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCGGGACGAGCCGACCGCCAGTTGGGTGa                      >  1:1754073/1‑238 (MQ=255)                                                                                                                                                                                                                           ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGGGG‑GGGGGACAATCGGCGAGGAAAATTAGGTTAGCCGGGTGGTCTCCGGGCCGACCGCCCTTCGCGTGCTGAAAAAATTTCCTACCCCTTAAATTC‑GCAAACACCATCTCTCGTCGCCGGAAGTGCCCCATTTTGCTCGAGAAACGCTGGGCGAGCCGGCCGCCCGATGGGTggg                     >  2:1328580/1‑237 (MQ=255)                                                                                                                                                                                                                            cgcTGCTGGGGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCGGGCCGGACTGCGGCGTGT‑GGTGGCAAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAg                     <  1:2223693/238‑1 (MQ=255)                                                                                                                                                                                                                             gctgctGGCGGGGTGGGCGACTATGGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGCCGAAAAATATCAGGTTAGCCGAATTTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCa                   <  1:169893/239‑1 (MQ=255)                                                                                                                                                                                                                               tgctgGCGGGGATGGCGACTATGGTTTACGAAGGATTGCCGACCTGGCCGGCCTGCGGCGTGT‑GGTGGACAATCGTCGAAAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCaa                  <  1:1871336/238‑1 (MQ=255)                                                                                                                                                                                                                                              cGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGCGAGCAATACGCTGGATGTGCCg  >  1:653411/1‑239 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 GGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGT‑GGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTC‑GCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCG  >  NC_007779/352597‑353072

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