Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,071,280	G→A	7.9%	L100L (CTG→TTG)	phoU ←	phosphate transport‑associated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,071,280	0	G	A	7.9%	177.1 / 4.2	63	L100L (CTG→TTG)	phoU	phosphate transport‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base A (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-02

AAAGGCATCCAGCGCATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTGGGC  >  NC_002947/6071134‑6071424                                                                                                                                                   |                                                                                                                                                  aaaGGCATCCAGCGCATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGt                                                                                                                                                 <  4:179398/149‑1 (MQ=255) aaaGGCATCCAGCGCATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGt                                                                                                                                                 >  7:5400/1‑149 (MQ=255)           aGCGCATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGa                                                                                                                                      >  3:182205/1‑149 (MQ=255)            gcgcATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAAGTCgg                                                                                                                                              <  8:15053/141‑2 (MQ=255)            gcgcATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAAGTCgg                                                                                                                                              >  7:15053/1‑140 (MQ=255)               cATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTg                                                                                                                                  >  2:88159/1‑149 (MQ=255)                   ccGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACttcctc                                                                                                                               >  6:128415/1‑144 (MQ=255)                   ccGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAAGTCGGTGA‑TCGACCTGGaac                                                                                                                              <  3:42894/149‑2 (MQ=255)                   ccGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCAACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAAGTCGGTGA‑TCGACCTGGaac                                                                                                                              <  2:110287/149‑2 (MQ=255)                         cATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTAACCGATGAGTTCCAGGTCTATAA‑CCGACTTCGAAATACtga                                                                                                                        >  6:259756/1‑149 (MQ=255)                          aTGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACtgat                                                                                                                       <  7:129398/149‑1 (MQ=255)                             ttGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACtgagga                                                                                                                     >  3:8387/1‑148 (MQ=255)                                cgcACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCGGGTCGATCA‑CCGACTTGGAAATAc                                                                                                                           >  6:187166/1‑139 (MQ=255)                                cgcACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCGGGTCGATCA‑CCGACTTGGAAATAc                                                                                                                           <  5:187166/139‑1 (MQ=255)                                 gcACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGATTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATa                                                                                                                            <  7:135918/137‑1 (MQ=255)                                 gcACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGAGTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATa                                                                                                                            >  8:135918/1‑137 (MQ=255)                                      tGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAg                                                                                                           >  4:303325/1‑149 (MQ=255)                                           gCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTcc                                                                                                       >  5:107295/1‑148 (MQ=255)                                            ccAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGCCTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTcc                                                                                                       >  3:171274/1‑147 (MQ=255)                                            ccAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAg                                                                                                                 >  4:276641/1‑137 (MQ=255)                                            ccAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAg                                                                                                                 <  3:276641/137‑1 (MQ=255)                                             cAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCC‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAg                                                                                                    >  8:152196/1‑149 (MQ=255)                                              aaTGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAgg                                                                                                   >  4:269517/1‑149 (MQ=255)                                                     cGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCgctgg                                                                                            >  4:146453/1‑149 (MQ=255)                                                      gAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCgctggc                                                                                           <  6:107295/149‑1 (MQ=255)                                                              aCATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGgcggcg                                                                                   <  4:182205/149‑1 (MQ=255)                                                                aTAACCACGCGGCGACTCGCCCTCTTCGCACCGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGAGCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCggcggg                                                                                 <  5:259756/149‑1 (MQ=255)                                                                   aCCACGCGGCGACTCGCCCTCTACGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAg                                                                                                    <  5:67811/127‑1 (MQ=255)                                                                   aCCACGCGGCGACTCGCCCTCTACGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAg                                                                                                    >  6:67811/1‑127 (MQ=255)                                                                    ccACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAg                                                                             >  7:236839/1‑149 (MQ=255)                                                                    ccACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAg                                                                             <  3:269517/149‑1 (MQ=255)                                                                    ccACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAg                                                                             >  6:222768/1‑149 (MQ=255)                                                                     cACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCACCCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAg                                                                             <  1:88159/149‑1 (MQ=255)                                                                     cACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAgg                                                                            <  8:5400/149‑1 (MQ=255)                                                                     cACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAg                                                                             >  6:272832/1‑148 (MQ=255)                                                                           gcgACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATgcgc                                                                      <  7:225301/149‑1 (MQ=255)                                                                             gagtCGCCCTCTTCGCACAGGTGGATGGCGCGGCGGGCGTTCTTGGTCGACTCGTCACCGATGAGTTCTAGGTCGATCT‑CCGACTTGGAAATACTGATGATCAGGCGCTGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATgcttgc                                                                    <  4:8387/146‑5 (MQ=255)                                                                                 cGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTCGCGGCGGGCGAGCATGagcaagca                                                                 >  6:10479/1‑148 (MQ=255)                                                                                  gCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGAACAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCAct                                                               >  6:166637/1‑149 (MQ=255)                                                                                   cccTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGAAct                                                               >  7:248368/1‑148 (MQ=255)                                                                                        ttCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGAGGGCGAGGATGCTCAAGCACTCAtcg                                                          >  6:212889/1‑148 (MQ=255)                                                                                          cGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGAGTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCtcgt                                                         <  5:263527/147‑1 (MQ=255)                                                                                          cGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCtcgtc                                                        >  7:247904/1‑148 (MQ=255)                                                                                           gCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATAc                                                                                                                           >  6:213210/1‑80 (MQ=255)                                                                                           gCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGa                                                      >  7:51634/1‑149 (MQ=255)                                                                                            cacaGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGGCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGAt                                                     <  8:51634/149‑1 (MQ=255)                                                                                                gCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTg                                                 <  8:247904/149‑1 (MQ=255)                                                                                                  tGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTgcg                                               >  1:125756/1‑149 (MQ=255)                                                                                                    gacggCGCGGCGGGCGATGTAGGTCGACTCGTCACCGATGCGTTCCAGGTCCATCA‑CCGCCTTGGAAATACTGATGATCATGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTgcgc                                              <  8:85873/145‑1 (MQ=255)                                                                                                      tGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTcc                                           <  2:256694/149‑1 (MQ=255)                                                                                                          gcggcgGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCt                                       <  5:222768/149‑1 (MQ=255)                                                                                                           cggcggGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTg                                      <  3:146453/149‑1 (MQ=255)                                                                                                            ggcgggctATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTgg                                     <  5:212889/149‑1 (MQ=255)                                                                                                                      cTTGGTCGTCTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGc                            <  4:171274/148‑1 (MQ=255)                                                                                                                         ggTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTg                             <  2:125756/144‑1 (MQ=255)                                                                                                                             gACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAAc                    <  4:121990/149‑1 (MQ=255)                                                                                                                                  gTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTc                 <  5:166637/147‑1 (MQ=255)                                                                                                                                       cGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCAcc          >  4:203335/1‑149 (MQ=255)                                                                                                                                          gtatttCCAAGTCGGTGA‑TCGACCTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACctgc       >  2:74097/5‑149 (MQ=255)                                                                                                                                             gTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGGGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTgg    >  1:92203/1‑149 (MQ=255)                                                                                                                                             gTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTgg    <  3:203335/149‑1 (MQ=255)                                                                                                                                                ccAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTGGGc  <  8:236839/148‑1 (MQ=255)                                                                                                                                                 cAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTGGGc  >  3:90911/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AAAGGCATCCAGCGCATCCCGGACCATGTTGCGCACCTGGTCGCCAATGTGGCGAACCTCGACATAACCACGCGGCGACTCGCCCTCTTCGCACAGCTGGATGGCGCGGCGGGCGATCTTGGTCGACTCGTCACCGATGCGTTCCAGGTCGATCA‑CCGACTTGGAAATACTGATGATCAGGCGCAGGTCCGAGGCCGCTGGCTGGCGGCGGGCGAGGATGCGCAAGCACTCCTCGTCGATGTTGCGCTCCATCTGGTTGATCTGCTCGTCAACCTCGCGCACCTGCTGGGC  >  NC_002947/6071134‑6071424

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