Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,890,439	T→A	11.0%	P57P (CCA→CCT) ‡	PP_2544 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,890,439	0	T	A	11.0%	219.5 / 14.4	82	P57P (CCA→CCT) ‡	PP_2544	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/31);  new base A (5/4);  total (47/35)
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.63e-01

GCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTGCG  >  NC_002947/2890295‑2890585                                                                                                                                                  |                                                                                                                                                      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:74536/149‑1 (MQ=255)           cTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGGTCCCCAAGGGAGCCGGCCGt                                                                                                                                         <  4:69450/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                                                                                                <  4:287569/102‑1 (MQ=255)                                                                                                  agcaagcaGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCAtgtg                                                                                                >  3:287569/1‑102 (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)                                                                                                                                              cccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGc                                                  >  4:299607/1‑103 (MQ=255)                                                                                                                                              cccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGc                                                  <  3:299607/103‑1 (MQ=255)                                                                                                                                               ccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTgc   >  5:370571/1‑149 (MQ=255)                                                                                                                                               ccTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTgc   >  5:123389/1‑149 (MQ=255)                                                                                                                                                 ttGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTgcg  >  6:180155/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GCCTTCATGCCTGACTACCAGATCCTGGTCATTTCCACCGGTATCGCTGCGGCGGTCGTGGTGGCGATGGGTGTGGTGCACCAG‑ATCCGCTCTGGTAGCAAGCAGCACTAAGCGTCACTCGCCTCTGGAAACGGCCGGCTCCCTTGGGGACCCGGCCGTTTCGCGTTGTTGGGCCTGGGCTTTTTCTGCGG‑‑GCATGTGGAAACGCTGGTTC‑ATCCACGGTGAACCGAGCAGGGCTGCCGCAGCGATCAGGGCCAGGATGGCGAGGGCGATGTAGGTCTTCACGGGCGTGCG  >  NC_002947/2890295‑2890585

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