Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,890,433	C→G	12.0%	E59D (GAG→GAC)	PP_2544 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,890,433	0	C	G	12.0%	201.4 / 17.1	75	E59D (GAG→GAC)	PP_2544	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/27);  new base G (5/4);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

ACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGG  >  NC_002947/2890289‑2890578                                                                                                                                                  |                                                                                                                                                     aCCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTccc                                                                                                                                                  >  2:238159/1‑149 (MQ=255)       gCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTggg                                                                                                                                             <  1:111055/148‑1 (MQ=255)        ccTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGa                                                                                                                                           >  3:318687/1‑149 (MQ=255)             aTGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCgg                                                                                                                                      >  4:107154/1‑149 (MQ=255)                 cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                                  <  4:69450/149‑1 (MQ=255)                 cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                                  <  4:74536/149‑1 (MQ=255)                 cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                                  >  3:74536/1‑149 (MQ=255)                 cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                                  >  3:69450/1‑149 (MQ=255)                 cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGAGGCGATGGCTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCACCTCTGGAAACGGCCGGCTCCCTTGGGGCCCCGGCCGt                                                                                                                                  >  5:260132/1‑149 (MQ=255)                  tGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGtt                                                                                                                                 >  1:198904/1‑149 (MQ=255)                    aCTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGGGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTc                                                                                                                               <  3:187255/149‑1 (MQ=255)                      tACCAGATCCTGGTCATTTCCCCCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGAGCACCAGTGACCGCAATGGTAGCAAGCAGCACTAAGCGTCAATCCCCTCTGGAAACGGCCGCCAACCTTGGGGACCCGGCCGtttgg                                                                                                                              >  1:293520/1‑147 (MQ=255)                        ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTccag                                                                                                                            <  2:273522/148‑4 (MQ=39)                        ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTcgc                                                                                                                             <  3:107154/147‑1 (MQ=255)                        ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCgt                                                                                                                           >  7:189819/1‑149 (MQ=255)                        ccAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCAGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTccag                                                                                                                            >  1:273522/1‑145 (MQ=38)                        ccAGATCCTGGGCTTTTCCACGGGTATCGCTGCGCCGGTGGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACAAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTcca                                                                                                                             <  5:143814/147‑3 (MQ=37)                         cAGATCCTGGTTATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTcca                                                                                                                             >  3:243072/1‑144 (MQ=39)                         cAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGTTTccagag                                                                                                                          >  5:149877/1‑144 (MQ=38)                         cAGATCCTGGTCATCTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCgtt                                                                                                                          >  7:310278/1‑149 (MQ=255)                               cTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTggg                                                                                                                    >  5:322115/1‑149 (MQ=255)                                     aTTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGttgttg                                                                                                                      >  7:158931/1‑141 (MQ=255)                                     aTTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGttgttg                                                                                                                      <  8:158931/141‑1 (MQ=255)                                         ccACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCC‑TGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCtttt                                                                                                         >  5:16371/1‑149 (MQ=255)                                          cACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCttt                                                                                                          <  8:189819/148‑1 (MQ=255)                                                 aTCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCgg                                                                                                  >  6:284356/1‑149 (MQ=255)                                                  tCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCgg                                                                                                  <  8:77786/148‑1 (MQ=255)                                                    tcTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑Gc                                                                                              <  6:260132/147‑1 (MQ=255)                                                        cggcggTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCAtgtg                                                                                         <  6:322115/149‑1 (MQ=255)                                                            ggTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGaaa                                                                                     >  1:121547/1‑149 (MQ=255)                                                                 tggtggCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTg                                                                                <  8:310278/149‑1 (MQ=255)                                                                  ggtggCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTgg                                                                               <  1:238159/149‑1 (MQ=255)                                                                      gCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑Gc                                                                                              <  5:340931/130‑1 (MQ=255)                                                                      gCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑Gc                                                                                              >  6:340931/1‑130 (MQ=255)                                                                           gggTGTGGTGCACCAG‑ATCCGCTCTCGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCAc                                                                     >  5:103339/1‑149 (MQ=255)                                                                                 ggtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGaa                                                               >  5:323145/1‑149 (MQ=255)                                                                                  gtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGt                                                                              <  4:345952/134‑1 (MQ=255)                                                                                  gtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGt                                                                              >  3:345952/1‑134 (MQ=255)                                                                                  gtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAAc                                                              <  5:75441/149‑1 (MQ=255)                                                                                  gtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAAc                                                              >  6:115343/1‑149 (MQ=255)                                                                                  gtgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCC‑TGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAAc                                                              <  6:16371/148‑1 (MQ=255)                                                                                   tgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACg                                                                    <  7:319683/142‑1 (MQ=255)                                                                                   tgCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACg                                                                    >  8:319683/1‑142 (MQ=255)                                                                                     cACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCGGGGGTTTTTCTGCGG‑‑GCATGTGGAACCGCTGGTTC‑ATCCACGGTGAAcggc                                                           >  4:67770/1‑146 (MQ=255)                                                                                       ccAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCttt                                                                                                          >  8:218188/1‑103 (MQ=255)                                                                                       ccAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCttt                                                                                                          <  7:218188/103‑1 (MQ=255)                                                                                        cAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCATTTTCTGCGG‑‑GCATGTGGAAACGCTGCTTC‑ATCCACGGTGAACCGAg                                                          >  3:274939/1‑147 (MQ=255)                                                                                         aG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAg                                                       >  2:75695/1‑149 (MQ=255)                                                                                          g‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTgg                                                                                                               >  5:216591/1‑95 (MQ=255)                                                                                          g‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTgg                                                                                                               <  6:216591/95‑1 (MQ=255)                                                                                          g‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAgg                                                      >  6:78826/1‑149 (MQ=255)                                                                                          g‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAgg                                                      >  6:265076/1‑149 (MQ=255)                                                                                            aTCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAggg                                                     <  6:323145/149‑1 (MQ=255)                                                                                            aTCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTCGAAACGCTGGTTC‑ATACACGGTGAACCCAGCAggg                                                     >  4:146840/1‑149 (MQ=255)                                                                                             tCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGc                                                    >  7:239722/1‑149 (MQ=255)                                                                                              ccGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCt                                                   >  5:301779/1‑149 (MQ=255)                                                                                                gCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGc                                                 >  4:318914/1‑149 (MQ=255)                                                                                                    tGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCg                                               <  3:318914/147‑1 (MQ=255)                                                                                                    tGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTACCg                                               <  4:318687/147‑1 (MQ=255)                                                                                                        agcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCAtgtg                                                                                         >  3:287569/1‑102 (MQ=255)                                                                                                        agcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCAtgtg                                                                                         <  4:287569/102‑1 (MQ=255)                                                                                                        agcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGa                                         >  2:115806/1‑149 (MQ=255)                                                                                                             gcagcaCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAgg                                    <  1:75695/149‑1 (MQ=255)                                                                                                               agcaCTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGc                                  <  2:121547/149‑1 (MQ=255)                                                                                                                    tAAGCGTCACTCGCCTCTGGAAACGGCCCGCTACCTTGGGGACCCGGCCCTTTCGCGTTGTTTGGCCTGGGGGTTTTGTGCGG‑‑GCATGTGGAAAAGCAGGTTCAAAACCCGGTAAACC‑AGCAGGGCTGCCGCATCCATCAGCGCCacg                              >  4:205051/1‑146 (MQ=255)                                                                                                                      aGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATg                           >  7:9242/1‑149 (MQ=255)                                                                                                                         gTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGc                         >  7:209217/1‑148 (MQ=255)                                                                                                                            aCTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGGGAgg                     <  2:293520/149‑1 (MQ=255)                                                                                                                            aCTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAg                      <  2:322298/148‑1 (MQ=255)                                                                                                                                gCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGGGAGGATGTGGTAACCCATGGTG‑ATGCACGGTGGACCGGGCAGGGCTGCCGGCGGGATGAGGGCCAGGATgggggcggg                   >  1:111894/1‑142 (MQ=255)                                                                                                                                  ctctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTTCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTTAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGTATGGCGAGGGCGATg               <  6:103339/149‑1 (MQ=255)                                                                                                                                  ctctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATg               >  6:347550/1‑149 (MQ=255)                                                                                                                                  ctctGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATg               >  8:107289/1‑149 (MQ=255)                                                                                                                                     tGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTAGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATACACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg            >  2:203824/1‑149 (MQ=255)                                                                                                                                       gAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAg            <  5:265076/147‑1 (MQ=255)                                                                                                                                           cggccggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTc      <  8:9242/149‑1 (MQ=255)                                                                                                                                             gccggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCAc    >  1:50764/1‑149 (MQ=255)                                                                                                                                              ccggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACg   <  5:347550/149‑1 (MQ=255)                                                                                                                                               cggcTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACgg  <  5:115343/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     ACCATGGCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGG  >  NC_002947/2890289‑2890578

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: ‑