Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,819	T→A	10.2%	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	10.2%	183.1 / 10.8	69	Q176L (CAG→CTG)	PP_4911	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base A (3/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

GTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582672‑5582960                                                                                                                                                     |                                                                                                                                              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                                                                                                                     <  4:37133/145‑1 (MQ=255)                              ttCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                     >  3:37133/1‑145 (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                                                                                                            <  5:94001/119‑1 (MQ=255)                                                                 ccgccgTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            >  6:94001/1‑119 (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                                                          <  3:29246/149‑1 (MQ=255)                                                                                     gCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGtt                                                          <  7:119806/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                                                                         >  7:302228/5‑86 (MQ=38)                                                                                                                                      tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAg                                                                         <  8:302228/82‑1 (MQ=38)                                                                                                                                      tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  4:17540/143‑1 (MQ=255)                                                                                                                                      tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            >  3:17540/5‑147 (MQ=255)                                                                                                                                      tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  8:151414/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)                                                                                                                                                     |                                                                                                                                              GTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582672‑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: ‑