Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,816	T→A	10.1%	Q177L (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,816	0	T	A	10.1%	187.8 / 10.7	69	Q177L (CAG→CTG)	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 = 7.21e-01

CTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582668‑5582960                                                                                                                                                      |                                                                                                                                                 cTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGtgct                                                                                                                                                  <  1:204537/149‑1 (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                                                                                                                     <  4:37133/145‑1 (MQ=255)                                  ttCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCa                                                                                                                     >  3:37133/1‑145 (MQ=255)                                    cGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGc                                                                                                                  <  7:230570/146‑1 (MQ=255)                                         aTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            >  6:332756/1‑147 (MQ=255)                                         aTCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATgg                                                                                                            <  5:332756/147‑1 (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                                                          <  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                                                   <  3:115356/144‑1 (MQ=255)                                                                                                     aaCAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCATTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCa                                                   >  4:115356/1‑144 (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                                                                                             >  4:279699/1‑79 (MQ=255)                                                                                                                            cTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGc                                                                                             <  3:279699/79‑1 (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            >  3:17540/5‑147 (MQ=255)                                                                                                                                          tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            >  7:151414/5‑147 (MQ=255)                                                                                                                                          tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  4:17540/143‑1 (MQ=255)                                                                                                                                          tccaccTCGAGCAGCAGCACGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTcc            <  8:151414/143‑1 (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)                                                                                                                                                      |                                                                                                                                                 CTGAGTCGGTGCAGCGGGCTTTTTCGGTGCCTCTTCGGCAATCGCCGGCTGGTTGGCCAGCGGCTTCACCGCCGTGCTCAGCTGTTCGGCCAGCTTCTGCAACAGTTGGCCCTGAGCCTGCACCTGCG‑ACGACTCACTGCCCTCGTGCTGCTGCTCGAGGTGGATGATACGGTTGTCACGCACATGGCCACGGCGGTCAAGCAGACGCCACTGGGCATCCAGGATGGCCGGCTGGTTCTTGCCAGAATCCAGGCGTGTGATCGACAACAGCACTTGCACGTCCGGGCTGAAGC  >  NC_002947/5582668‑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: ‑