Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,488,203	G→A	5.3%	intergenic (+59/‑343)	PP_2182 → / → PP_2183	hypothetical protein/formate dehydrogenase subunit gamma

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,488,203	0	G	A	5.3%	284.2 / 2.3	93	intergenic (+59/‑343)	PP_2182/PP_2183	hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (41/47);  new base A (3/2);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

CGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCT  >  NC_002947/2488061‑2488342                                                                                                                                                |                                                                                                                                              cGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTgg                                                                                                                                          >  1:45575/1‑148 (MQ=255)  gACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCa                                                                                                                                        <  3:255275/149‑1 (MQ=255)    cGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAgg                                                                                                                                      <  5:347896/149‑1 (MQ=255)      cAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGcc                                                                                                                                    <  6:237270/149‑1 (MQ=255)        aTCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCg                                                                                                                                  <  2:289034/149‑1 (MQ=255)         tCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCgg                                                                                                                                 <  1:347857/149‑1 (MQ=255)            tGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCaa                                                                                                                              >  1:244270/1‑149 (MQ=255)              cGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATg                                                                                                                            <  4:65993/149‑1 (MQ=255)                     aaaGGCCCCGGCCCAGGCTCCGGTCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAgg                                                                                                                    <  5:90010/149‑1 (MQ=255)                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAgg                                                                                                                    >  3:263899/1‑149 (MQ=255)                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAgg                                                                                                                    >  7:192380/1‑149 (MQ=255)                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAgg                                                                                                                    >  2:268296/1‑149 (MQ=255)                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAg                                                                                                                     <  4:53893/148‑1 (MQ=255)                       aggcccaggcccAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGtt                                                                                                                  <  2:244270/149‑1 (MQ=255)                          cccaggcccaggCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCgg                                                                                                               <  3:7057/149‑1 (MQ=255)                           ccaggcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGt                                                                                                              >  1:347732/1‑149 (MQ=255)                           ccaggcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGt                                                                                                              <  1:268296/149‑1 (MQ=255)                                cccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCga                                                                                                         <  3:271035/149‑1 (MQ=255)                                 ccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTc                                                                                                           >  2:16592/1‑146 (MQ=255)                                  catgcTCCGGCCCCCGCCGCCGACACCGACACGGTGAAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTATCCATGCTGTATGGCATCCCATTGCCGGGGCTTCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCgaga                                                                                                       <  4:193266/146‑1 (MQ=255)                                   aggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGagag                                                                                                      >  4:21381/1‑149 (MQ=255)                                    ggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAgcgc                                                                                                     >  8:170849/1‑149 (MQ=255)                                       tCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTg                                                                                                  >  5:139645/1‑149 (MQ=255)                                       tCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTg                                                                                                  >  6:352545/1‑149 (MQ=255)                                             ccccGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATCGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGtgcat                                                                                            >  8:216905/1‑149 (MQ=255)                                                cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCTTTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                           >  1:98761/1‑148 (MQ=255)                                                cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGATCTGCCCCTGAAGCCATGCTCGATGGCACCCCATTGCCGGGGCTGCCGTCTGGCAGCCCCCGCAATGGA‑CAAGGTTCGGTTTCCAGAGCCTGGtgcatg                                                                                           >  1:103155/1‑148 (MQ=255)                                                    gccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGa                                                                                     <  5:352545/149‑1 (MQ=255)                                                     ccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                                    >  7:117937/1‑149 (MQ=255)                                                     ccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                                    <  6:306478/149‑1 (MQ=255)                                                     ccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATCGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                                    <  7:216905/149‑1 (MQ=255)                                                       gACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATgg                                                                                  >  2:325471/1‑149 (MQ=255)                                                           ccGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGttt                                                                               >  3:136827/1‑148 (MQ=255)                                                           ccGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTc                                                                              <  5:198677/149‑1 (MQ=255)                                                              cAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAg                                                                           <  4:263899/149‑1 (MQ=255)                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGTCTGCCGTCTGGCAGCCCCGGCAATGGA‑CTAGCATCGGTTTAGACAGCCTGGTGCATCCAACGATGGTTTCAAGa                                                                          >  2:111716/1‑149 (MQ=255)                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGa                                                                          <  2:45575/149‑1 (MQ=255)                                                                    tGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                            >  2:255833/1‑142 (MQ=255)                                                                    tGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                            <  1:255833/142‑1 (MQ=255)                                                                    tGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                                     <  3:21381/149‑1 (MQ=255)                                                                    tGCAGGCTCCCGACTGAGCTCAGTCGCGGTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                                     <  2:318621/149‑1 (MQ=255)                                                                      cattCTACCGACTGAGCTCATTCGCGTTCCGCCCC‑GTAGCCATGCTGGATGGCATCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATTGACCTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                                     <  1:270360/145‑1 (MQ=255)                                                                      cAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGCGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCt                                                                   >  2:115641/1‑149 (MQ=255)                                                                           tACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCCACGATGGTTTCAAGACCCACCTG‑ATCg                                                             >  5:249559/1‑149 (MQ=255)                                                                            aCCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                            >  3:52758/1‑134 (MQ=255)                                                                            aCCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                            <  4:52758/134‑1 (MQ=255)                                                                             ccGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                                     >  4:65421/1‑140 (MQ=255)                                                                             ccGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                                     <  3:65421/140‑1 (MQ=255)                                                                             ccGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCg                                                             <  7:88416/147‑1 (MQ=255)                                                                             ccGACTGAGCTCAGGCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCg                                                             <  5:366017/147‑1 (MQ=255)                                                                               gACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGc                                                           >  6:28337/1‑147 (MQ=255)                                                                               gACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCaa                                                         <  8:192380/149‑1 (MQ=255)                                                                                 cTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGc                                                       <  1:325471/149‑1 (MQ=255)                                                                                     gCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTg                                                   <  2:347732/149‑1 (MQ=255)                                                                                      cTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTgg                                                  <  1:16592/149‑1 (MQ=255)                                                                                      cTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTgg                                                  >  2:343213/1‑149 (MQ=255)                                                                                      cTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑AACGGCAAGCCGTgg                                                  >  5:308295/1‑149 (MQ=255)                                                                                      cTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCAGGCAGGCACGGCAATCGA‑CAAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATAGTATCAAGACACACCTG‑ATCGGCAAGCCGtcg                                                  >  6:5821/1‑147 (MQ=255)                                                                                         aGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTgg                                                  >  8:171400/1‑146 (MQ=255)                                                                                         aGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAg                                               <  5:28337/149‑1 (MQ=255)                                                                                           tCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAgat                                             <  4:136827/149‑1 (MQ=255)                                                                                           tCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGCGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACACACCTG‑AACGGCAAGCCGAGGCAgat                                             >  7:284937/1‑149 (MQ=255)                                                                                                     cccc‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCgggg                                   >  1:38602/1‑149 (MQ=255)                                                                                                     cccc‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCggg                                    <  2:260009/148‑1 (MQ=255)                                                                                                     cccc‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCAACCTG‑ATCGGCAAGCCGTGGCAGATGATGCAgggg                                   >  6:222509/1‑149 (MQ=255)                                                                                                      ccc‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                  <  1:343213/149‑1 (MQ=255)                                                                                                       cc‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCg                                 >  7:191620/1‑149 (MQ=255)                                                                                                              ccATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                                    >  6:61944/1‑93 (MQ=255)                                                                                                               cATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTgg                                                                                                 >  3:249799/1‑79 (MQ=255)                                                                                                               cATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTgg                                                                                                 <  4:249799/79‑1 (MQ=255)                                                                                                                aTGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACg                         >  6:159855/1‑149 (MQ=255)                                                                                                                    tGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCg                     >  6:301438/1‑149 (MQ=255)                                                                                                                      gATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGt                    <  2:36841/148‑1 (MQ=255)                                                                                                                      gATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTc                   <  6:308295/149‑1 (MQ=255)                                                                                                                         ggCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGGGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCggg                >  2:229975/1‑148 (MQ=255)                                                                                                                            agtccATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGG‑CTAGGGTCGGTTTCGAGAGCCAGGTCCATGCAACGATGGTTTCAAGACCCACATG‑ATCGGCAAGGCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGgtt              >  4:277720/4‑148 (MQ=25)                                                                                                                            aGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                  <  3:9122/128‑1 (MQ=255)                                                                                                                            aGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                  >  4:9122/1‑128 (MQ=255)                                                                                                                            aGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGt                    >  2:73746/1‑142 (MQ=255)                                                                                                                            aGCCCATTGCCGGGGCGGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGTTTCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGt                    <  1:73746/143‑1 (MQ=255)                                                                                                                             gtccATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                  <  6:323577/125‑1 (MQ=37)                                                                                                                             gtccATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                  >  5:323577/3‑127 (MQ=37)                                                                                                                              tccATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGG‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGttt              >  4:334049/2‑146 (MQ=38)                                                                                                                                cATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGaa         <  2:38602/149‑1 (MQ=255)                                                                                                                                     ccGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATgcgc    <  5:159855/149‑1 (MQ=255)                                                                                                                                      cGGGGCTGCCGTCTGGCAGGCTCGGCAATGTA‑CTTGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCATAATGCGcc   <  5:5821/149‑1 (MQ=255)                                                                                                                                      cGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTgg                                                  <  4:352622/102‑1 (MQ=255)                                                                                                                                      cGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTgg                                                  >  3:352622/1‑102 (MQ=255)                                                                                                                                       ggggCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGc                                                       <  6:168677/96‑1 (MQ=255)                                                                                                                                       ggggCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGc                                                       >  5:168677/1‑96 (MQ=255)                                                                                                                                       ggggCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCt  <  2:98761/149‑1 (MQ=255)                                                                                                                                       ggggCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCt  <  2:103155/149‑1 (MQ=255)                                                                                                                                       gcgccggCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGATCGGGTTCCAGAGCCTGGTGCATGCAACGACGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATCCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCt  <  8:331581/143‑1 (MQ=255)                                                                                                                                                |                                                                                                                                              CGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCC‑GTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGA‑CTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTG‑ATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCT  >  NC_002947/2488061‑2488342

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