Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,819	T→A	6.8%	Q176L (CAG→CTG)	PP_4911 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,582,819	0	T	A	6.8%	251.3 / 7.4	88	Q176L (CAG→CTG)	PP_4911	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base A (4/2);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

GTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCGG  >  NC_002947/5582676‑5582966                                                                                                                                                |                                                                                                                                                    gTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTc                                                                                                                                                 >  1:83651/1‑148 (MQ=255)    cAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAgg                                                                                                                                             <  5:391813/149‑1 (MQ=255)    cAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAg                                                                                                                                              <  6:141136/148‑1 (MQ=255)       cGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTgg                                                                                                                                          <  5:399522/149‑1 (MQ=255)       cGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTg                                                                                                                                           >  7:333337/1‑148 (MQ=255)          gCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgat                                                                                                                                        >  6:71778/1‑148 (MQ=255)           cTTTTTCGGTGCCTCTTCGGCATTGGCCGGGTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgatga                                                                                                                                      <  8:333337/149‑1 (MQ=255)            tttttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGCCTCACTGCCCTCGTGCTGCTGCTCGAGGTgactac                                                                                                                                      >  5:234203/1‑143 (MQ=255)            tttttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgatga                                                                                                                                      >  3:214417/1‑148 (MQ=255)            tttttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGgatga                                                                                                                                      >  6:416682/1‑148 (MQ=255)              tttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATa                                                                                                                                    <  1:409091/148‑1 (MQ=255)                 cGGTGCCTCTTCGGCAATCGCCGGCTGATTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGt                                                                                                                                >  1:133284/1‑149 (MQ=255)                  ggTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGtt                                                                                                                               <  1:273620/149‑1 (MQ=255)                  ggTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGtt                                                                                                                               <  2:40281/149‑1 (MQ=255)                   gTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATCCGGTTg                                                                                                                              >  2:141625/1‑149 (MQ=255)                   gTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTg                                                                                                                              >  7:292641/1‑149 (MQ=255)                   gTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGAGCAGCAGCACGAgg                                                                                                                                             >  6:137707/1‑134 (MQ=255)                   gTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGAGCAGCAGCACGAgg                                                                                                                                             <  5:137707/134‑1 (MQ=255)                       ctctTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCAc                                                                                                                          <  1:244189/149‑1 (MQ=255)                        tctTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACg                                                                                                                         >  5:334517/1‑149 (MQ=255)                          ttCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGca                                                                                                                       >  5:63064/1‑149 (MQ=255)                             ggCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTc                                                                                                                            <  7:214656/141‑1 (MQ=255)                             ggCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTc                                                                                                                            >  8:214656/1‑141 (MQ=255)                              gCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATg                                                                                                                   >  8:193046/1‑149 (MQ=255)                                aaTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAgg                                                                                                                                             <  3:36593/121‑1 (MQ=255)                                aaTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAgg                                                                                                                                             >  4:36593/1‑121 (MQ=255)                                aaTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTATGCAACAGTTGGCCCTGAGACGGCACCGGAGACGTCTCACTGCCCTCGTGCGGCGTCTCGAGGTGGATGCTACGCTTGTAACGcccaggca                                                                                                                 >  4:310322/1‑145 (MQ=255)                                     ccGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATa                                                                                                                                    <  3:225264/125‑1 (MQ=255)                                     ccGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATa                                                                                                                                    >  4:225264/1‑125 (MQ=255)                                       ggCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcggcg                                                                                                          <  2:83651/149‑1 (MQ=255)                                        gCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcggcgg                                                                                                         <  5:71778/149‑1 (MQ=255)                                            gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTgg                                                                                                                                          >  8:151449/1‑112 (MQ=255)                                            gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTgg                                                                                                                                          <  7:151449/112‑1 (MQ=255)                                                gCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcggcg                                                                                                          >  4:65286/1‑140 (MQ=255)                                                gCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcggcg                                                                                                          <  3:65286/140‑1 (MQ=255)                                                 ccAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGa                                                                                                >  1:206353/1‑149 (MQ=255)                                                 ccAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGa                                                                                                >  1:299146/1‑149 (MQ=255)                                                       gCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCAc                                                                                          >  8:117713/1‑149 (MQ=255)                                                             ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCa                                                                                    >  1:114502/1‑149 (MQ=255)                                                              cgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcg                                                                                                             >  8:84056/1‑123 (MQ=255)                                                              cgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcg                                                                                                             <  7:84056/123‑1 (MQ=255)                                                                   tGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAc                                                                                                              >  4:324568/1‑117 (MQ=255)                                                                   tGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAc                                                                                                              <  3:324568/117‑1 (MQ=255)                                                                      tCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCATGAGCCTGCACATGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCGCGGCGGTCAAGCAGACCCCACTGGCCATCCAGGATg                                                                           >  6:47365/1‑149 (MQ=255)                                                                         gCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGcc                                                                        <  5:416682/149‑1 (MQ=255)                                                                            tcagggCCAACTGTTGCAGAAGCTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     >  5:276174/5‑149 (MQ=255)                                                                            gTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     >  3:201915/1‑149 (MQ=255)                                                                               cGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     <  4:71309/146‑1 (MQ=255)                                                                               cGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     >  3:71309/1‑146 (MQ=255)                                                                               cGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTgg                                                                  <  4:214417/149‑1 (MQ=255)                                                                                ggCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGt                                                                 <  8:292641/149‑1 (MQ=255)                                                                                 gCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGtt                                                                >  1:220548/1‑149 (MQ=255)                                                                                 gCCAGCTTCTGCAACAGTTGGCCATGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGtt                                                                <  8:224204/149‑1 (MQ=255)                                                                                     gCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTg                                                            <  4:201915/149‑1 (MQ=255)                                                                                         cTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCa                                                         <  2:299146/148‑1 (MQ=255)                                                                                           gCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATgg                                                                          <  6:118870/129‑1 (MQ=255)                                                                                           gCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATgg                                                                          >  5:118870/1‑129 (MQ=255)                                                                                           gCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGaa                                                      >  3:232476/1‑149 (MQ=255)                                                                                            cAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAAt                                                     >  8:231438/1‑149 (MQ=255)                                                                                                  ctGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCGCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGc                                               <  6:276174/148‑1 (MQ=255)                                                                                                     gCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgt                                             <  3:56433/148‑1 (MQ=255)                                                                                                      cccTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtgt                                           <  6:148523/149‑1 (MQ=255)                                                                                                       ccTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtgt                                           >  6:352273/1‑148 (MQ=255)                                                                                                        cTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGa                                         >  7:40845/1‑149 (MQ=255)                                                                                                             ccTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCAAGGCGTGTGATCGa                                     >  5:133873/1‑148 (MQ=255)                                                                                                                 cACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGacaaca                                <  7:117713/149‑1 (MQ=255)                                                                                                                    cTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCa                             >  6:232072/1‑149 (MQ=255)                                                                                                                       cgacgaCTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACt                           <  8:40845/148‑1 (MQ=255)                                                                                                                         acgacTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGc                        <  6:63064/149‑1 (MQ=255)                                                                                                                             cTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGt                    >  5:157403/1‑149 (MQ=255)                                                                                                                              tCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGCCGCCACTGGGCATCCAGGATGGCCGGCTCGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTc                   >  3:80493/1‑149 (MQ=255)                                                                                                                              tCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGa                                     >  6:278233/1‑131 (MQ=255)                                                                                                                              tCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGa                                     <  5:278233/131‑1 (MQ=255)                                                                                                                              tCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTc                   >  1:383914/1‑149 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGTATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACtt                          <  6:267254/140‑1 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     <  6:133838/97‑1 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGc                                                                     >  5:133838/1‑97 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCa                             >  5:37084/1‑137 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCa                             <  6:37084/137‑1 (MQ=255)                                                                                                                                aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACtt                          >  5:267254/1‑140 (MQ=255)                                                                                                                                 cTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCgg                >  4:557/1‑149 (MQ=255)                                                                                                                                  tGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCggg               >  6:406249/1‑149 (MQ=255)                                                                                                                                  ccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGc                        >  5:92744/4‑140 (MQ=255)                                                                                                                                  ccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGc                        <  6:92744/137‑1 (MQ=255)                                                                                                                                   caccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGAGATCGACAACAGCACTTGCACGTCCGGGc              >  3:125849/3‑149 (MQ=255)                                                                                                                                    cccTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACAT‑GCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCt             >  8:35843/1‑148 (MQ=255)                                                                                                                                     ccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCt             >  6:25046/1‑148 (MQ=255)                                                                                                                                      cTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTg            >  2:322206/1‑148 (MQ=255)                                                                                                                                               ctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCgg  >  7:22722/1‑149 (MQ=255)                                                                                                                                               ctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCg   <  1:141625/148‑1 (MQ=255)                                                                                                                                               ctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCg   <  7:231438/148‑1 (MQ=255)                                                                                                                                               ctgctCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTCGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCgg  >  3:173311/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGCCTGCGG  >  NC_002947/5582676‑5582966

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