Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,813	T→A	10.1%	H178L (CAC→CTC)	PP_4911 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,582,813	0	T	A	10.1%	183.5 / 10.7	69	H178L (CAC→CTC)	PP_4911	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base A (3/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

ACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582666‑5582960                                                                                                                                                     |                                                                                                                                                    aCCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGt                                                                                                                                                     <  8:206644/148‑1 (MQ=255)   cTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGtgct                                                                                                                                                  <  1:204537/149‑1 (MQ=255)    tGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGagcag                                                                                                                                                 >  3:65792/1‑147 (MQ=255)       gTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGtgctgctg                                                                                                                                              <  7:168474/149‑1 (MQ=255)            tGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGa                                                                                                                                         <  8:342819/149‑1 (MQ=255)               aGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGAGCAGCAGCACGAggg                                                                                                                                      <  4:65792/149‑2 (MQ=255)               aGCGGGCGTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGt                                                                                                                                      <  7:273671/149‑1 (MQ=255)                      tttttCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGACGTGgatga                                                                                                                                >  2:110052/1‑148 (MQ=255)                               gcatctTCGGCAATCGCCGGCTGGTTGGCCAGAGGCTTCACAGCCGTGCTCAGCTGTTCGGCCAGCTTATGCAAAAGTTGGACCTGAGCCTGAACGTGCG‑ACGACTCACTGCCCTCGTTCTGAGGCTCGAGGTGGATCATACGGTTGt                                                                                                                       >  4:140829/4‑148 (MQ=255)                                  tctTCGGCAATCGCCGGCTGGTTGGCCAGCCGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACg                                                                                                                   >  1:30380/1‑149 (MQ=255)                                   ctTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTACCTGCCCGCGTGCTGCTGCTGGAGGTGGATGATACGGTTGTCACGc                                                                                                                  >  8:119806/1‑149 (MQ=255)                                    ttCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                     >  3:37133/1‑145 (MQ=255)                                    ttCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                     <  4:37133/145‑1 (MQ=255)                                      cGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGc                                                                                                                  <  7:230570/146‑1 (MQ=255)                                           aTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            <  5:332756/147‑1 (MQ=255)                                           aTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            >  6:332756/1‑147 (MQ=255)                                            tCGCCGGCTGGTTGGCAAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCa                                                                                                         <  5:230414/149‑1 (MQ=255)                                                cGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCGCTCGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCAcgg                                                                                                      <  1:110052/149‑1 (MQ=255)                                                      gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCaa                                                                                               <  2:205890/149‑1 (MQ=255)                                                      gTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACCGTTGTCACGCACATGGCCACGGCGGTc                                                                                                 >  3:280591/1‑147 (MQ=255)                                                          gCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCa                                                                                            <  5:190075/148‑1 (MQ=255)                                                           ccAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGa                                                                                          >  4:254949/1‑149 (MQ=255)                                                             aGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACg                                                                                        >  2:10223/1‑149 (MQ=255)                                                                      aCCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTggg                                                                                >  5:194796/1‑148 (MQ=255)                                                                       ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            >  6:94001/1‑119 (MQ=255)                                                                       ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            <  5:94001/119‑1 (MQ=255)                                                                       ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCa                                                                              <  4:271643/149‑1 (MQ=255)                                                                          ccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTGGAGGTGGATGCTACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATc                                                                            >  3:223915/1‑148 (MQ=255)                                                                           cgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCa                                                                          <  3:254949/149‑1 (MQ=255)                                                                             tGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAg                                                                         >  5:266265/1‑148 (MQ=255)                                                                                tCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATg                                                                     >  5:291211/1‑149 (MQ=255)                                                                                       ttCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGAGACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGc                                                                   >  7:193309/1‑144 (MQ=255)                                                                                       ttCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGAGACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGc                                                                   <  8:193309/144‑1 (MQ=255)                                                                                          ggCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGt                                                           >  4:262894/1‑149 (MQ=255)                                                                                           gCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGtt                                                          <  7:119806/149‑1 (MQ=255)                                                                                           gCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGtt                                                          <  3:29246/149‑1 (MQ=255)                                                                                            ccatgttCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTc                                                         <  3:140829/144‑1 (MQ=255)                                                                                            ccAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTc                                                         >  6:114344/1‑149 (MQ=255)                                                                                              aGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCGAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCtt                                                       >  8:226931/1‑149 (MQ=255)                                                                                              aGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCtt                                                       <  6:291211/149‑1 (MQ=255)                                                                                               gCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGCCTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTg                                                      >  6:197780/1‑149 (MQ=255)                                                                                                      cAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGc                                                     >  5:339829/1‑143 (MQ=255)                                                                                                      cAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGc                                                     <  6:339829/143‑1 (MQ=255)                                                                                                       aaCAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCATTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCa                                                   >  4:115356/1‑144 (MQ=255)                                                                                                       aaCAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCATTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCa                                                   <  3:115356/144‑1 (MQ=255)                                                                                                        aCAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATcg                                             >  8:230678/1‑148 (MQ=255)                                                                                                         cAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCa                                                                          >  8:118674/1‑119 (MQ=255)                                                                                                         cAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCa                                                                          <  7:118674/119‑1 (MQ=255)                                                                                                          aGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAg                                           >  8:1959/1‑149 (MQ=255)                                                                                                            ttGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGc                                         <  1:144133/149‑1 (MQ=255)                                                                                                               gCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtg                                      <  6:82650/149‑1 (MQ=255)                                                                                                                cccTGAGCCTGCACCTGCC‑ACGCCTCACTCCCCTCGTGCTCCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGCCCACGGCCGACAAGCAGACGCCACTGGGCATCCAGGATCCCCGCCTCGTTCTTGCCAGAATCCAGGCgtgt                                     >  4:157918/1‑149 (MQ=255)                                                                                                                 ccTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGGCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCgtgt                                     >  5:90364/1‑148 (MQ=255)                                                                                                                           ctccTGCG‑ACGACTCACTGGCCTGGTGCTGCTGCTCGTGGTGGTTGATTCGGTTGTATCGCTCATGGCCACGGCGGTGAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGacaaca                          <  3:157918/147‑1 (MQ=255)                                                                                                                              cTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGc                                                                                             <  3:279699/79‑1 (MQ=255)                                                                                                                              cTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGc                                                                                             >  4:279699/1‑79 (MQ=255)                                                                                                                                    acgacTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGATGTCACGCACAAGGCCACGGCCGTCAAGCAGACGCCACTCGGCATCCAGGAAGGCCGCCTCGTTCTTGCCAGAATCCAGGCGAGTGATCGACAACAGCACTTg                   >  5:62153/1‑148 (MQ=255)                                                                                                                                      gacTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCAc                >  2:45201/1‑149 (MQ=255)                                                                                                                                         tCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCGACGGCGGTCAAGCAGACGCCACTGGGCATACAGGATGGCCGGCTTGTTGTTGGCAGAATACAGGCGTGTGATCGACAACAGCAAATGCACGTc             >  1:285550/1‑149 (MQ=255)                                                                                                                                           aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGacaaca                          <  7:210396/134‑1 (MQ=255)                                                                                                                                           aCTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGACACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGacaaca                          >  8:210396/1‑134 (MQ=255)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAg                                                                         <  8:302228/82‑1 (MQ=38)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAg                                                                         >  7:302228/5‑86 (MQ=38)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  8:151414/143‑1 (MQ=255)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            >  3:17540/5‑147 (MQ=255)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  4:17540/143‑1 (MQ=255)                                                                                                                                            tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            >  7:151414/5‑147 (MQ=255)                                                                                                                                            cTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCgg          <  6:266265/149‑1 (MQ=255)                                                                                                                                               cccTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCt       <  6:194796/149‑1 (MQ=255)                                                                                                                                               cccTCGTGCTGCTGCTCGAGGTGGATGATACCGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCt       <  4:280591/149‑1 (MQ=255)                                                                                                                                                ccTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTg      <  2:30380/149‑1 (MQ=255)                                                                                                                                                  tCGTGCTGCTGCTCGAGGTGGATGATACGGTTGACACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTg      >  2:288061/1‑147 (MQ=255)                                                                                                                                                    gTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGc  <  1:45201/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ACCTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582666‑5582960

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