Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,384,527	G→C	6.1%	G331R (GGC→CGC) ‡	lldP →	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,384,527	0	G	C	6.1%	209.8 / 3.3	65	G331R (GGC→CGC) ‡	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CCGTTGTGCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCACCGCGAT  >  NC_002947/5384390‑5384672                                                                                                                                          |                                                                                                                                                   ccGTTGTGCAGCCACGCCGCACCCAGCCCAGCCCGTCCAGCTTTCGCGAGATCCTCCCCGCCTGGGCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGCCCGGCGCGaaa                                                                                                                                         <  4:211075/149‑3 (MQ=255)  cGTTGTGCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGt                                                                                                                                        <  8:2825/149‑1 (MQ=255)    ttGTGCAGCCAAGCGGCAGCCAGCCTAGCGCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTa                                                                                                                                       >  7:294613/1‑148 (MQ=255)    ttGTGCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTAc                                                                                                                                      <  7:106397/149‑1 (MQ=255)    ttGTGCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTAc                                                                                                                                      <  1:282558/149‑1 (MQ=255)        gCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGCCCGGcgcg                                                                                                                                            <  3:198320/139‑1 (MQ=255)        gCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGCCCGGcgcg                                                                                                                                            >  4:198320/1‑139 (MQ=255)         cAGCCAAGCGGCAGCCAGCCTAGCCCGTACCGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGCAACCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACtt                                                                                                                                 >  4:51241/1‑149 (MQ=255)           gCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCg                                                                                                                               <  3:52102/149‑1 (MQ=255)            ccAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCgt                                                                                                                              <  2:277128/149‑1 (MQ=255)               aGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGtt                                                                                                                           <  4:133752/149‑1 (MQ=255)                      gccagccTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTc                                                                                                                    <  2:30017/149‑1 (MQ=255)                              aGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATccc                                                                                                            >  5:204429/1‑149 (MQ=255)                                  cgcaCAGCTTTGTCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTGCAACTTCGTGTTCAACTTCGCCATCCCGCAc                                                                                                        <  2:133102/146‑1 (MQ=255)                                         cTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGAc                                                                                                  >  6:277332/1‑148 (MQ=255)                                             ggCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGgcggcg                                                                                                                                              <  7:82069/100‑1 (MQ=255)                                             ggCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGgcggcg                                                                                                                                              >  8:82069/1‑100 (MQ=255)                                               cGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTgg                                                                                          >  2:315435/1‑149 (MQ=255)                                               cGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTgg                                                                                          <  6:11334/149‑1 (MQ=255)                                               cGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTgg                                                                                          >  5:11334/1‑149 (MQ=255)                                                gagaTCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGt                                                                                         >  1:333220/1‑149 (MQ=255)                                                  gaTCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAATCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGa                                                                                       <  3:51241/149‑1 (MQ=255)                                                    tCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGAt                                                                                      <  3:172444/148‑1 (MQ=255)                                                      ttCACGGCCTGGTCGCCGTTCCTGATCCTCAACGTGCTGCTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATc                                                                                     <  4:201958/147‑1 (MQ=255)                                                         aaGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGa                                                                                 <  3:172584/148‑1 (MQ=255)                                                          aGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGAcc                                                                               >  7:77064/1‑149 (MQ=255)                                                           ggCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCg                                                                              <  7:23631/149‑1 (MQ=255)                                                            gCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTc                                                                                                     <  6:219451/126‑1 (MQ=255)                                                            gCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTc                                                                                                     >  5:219451/1‑126 (MQ=255)                                                                      cGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGtt                                                                   <  6:204429/149‑1 (MQ=255)                                                                           cTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCcgc                                                              <  4:285415/149‑1 (MQ=255)                                                                            tGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCcgcg                                                             >  5:154324/1‑149 (MQ=255)                                                                                 cTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGa                                                         <  8:77064/148‑1 (MQ=255)                                                                                  tCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTCTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGAcc                                                       >  7:61363/1‑149 (MQ=255)                                                                                    aCCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGCCCAG‑CTGGTGATCAAGACCGCACCGATCg                                                                     >  5:51430/1‑133 (MQ=255)                                                                                    aCCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCg                                                                     <  6:51430/133‑1 (MQ=255)                                                                                     ccGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGc                                                     >  6:190656/1‑148 (MQ=255)                                                                                                aCCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGgt                                         >  8:100776/1‑149 (MQ=255)                                                                                                  cATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGgtgt                                       <  7:100776/149‑1 (MQ=255)                                                                                                   aTCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCcgg                                             >  5:174767/1‑142 (MQ=255)                                                                                                   aTCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCcgg                                             <  6:174767/142‑1 (MQ=255)                                                                                                        gACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTc                                     <  8:294613/145‑1 (MQ=255)                                                                                                        gACACTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGc                                 >  2:100405/1‑149 (MQ=255)                                                                                                         aCCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCt                                >  2:214389/1‑149 (MQ=255)                                                                                                           ccTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTgg                                                                                          >  7:95485/1‑89 (MQ=255)                                                                                                           ccTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTgg                                                                                          <  8:95485/89‑1 (MQ=255)                                                                                                             tGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGAc                            >  1:277749/1‑149 (MQ=255)                                                                                                             tGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGAc                            >  2:290150/1‑149 (MQ=255)                                                                                                               aaGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGAccc                          >  1:147865/1‑149 (MQ=255)                                                                                                               aaGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGAc                            <  1:214389/147‑1 (MQ=255)                                                                                                                   cGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGAtc                      <  5:186545/149‑1 (MQ=255)                                                                                                                    gTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGAtct                     >  8:20053/1‑149 (MQ=255)                                                                                                                    gTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGAtct                     >  2:224811/1‑149 (MQ=255)                                                                                                                     ttCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCAATCCCGCACCTCGACCAG‑CTGCTGATCAAGACCGCACAGATCGATGCCGCGCCGACAGCCAACCCGGCGGTGATCAAACTCGACCCGAtct                     >  4:317410/1‑148 (MQ=255)                                                                                                                       cAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTcc                   <  7:20053/148‑1 (MQ=255)                                                                                                                        aaGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGa                                                                                 <  6:76857/85‑1 (MQ=255)                                                                                                                        aaGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGa                                                                                 >  5:76857/1‑85 (MQ=255)                                                                                                                           gcggcTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCAc              <  2:147865/149‑1 (MQ=255)                                                                                                                               cTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCg            <  3:314838/147‑1 (MQ=255)                                                                                                                               cTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCg            >  4:314838/1‑147 (MQ=255)                                                                                                                                tttCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCa         >  6:68078/1‑149 (MQ=255)                                                                                                                                 atCGCGCCGCCCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCAc        >  6:240003/2‑149 (MQ=255)                                                                                                                                     cgcCGGGCGGCGCGATGTACAACTTCGTGTTCACCTTCGCCATCCCGCACCTCGACCAGAGTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCACcgc     >  1:271953/1‑149 (MQ=255)                                                                                                                                      gcCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCACCGCGa   >  3:286182/1‑149 (MQ=255)                                                                                                                                       ccGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCACCGCGAt  >  6:310370/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                   CCGTTGTGCAGCCAAGCGGCAGCCAGCCTAGCCCGTACAGCTTTGGCGAGATCTTCAAGGCCTGGTCGCCGTTCCTGATCCTCACCGTGCTGGTCACCATCTGGACCCTGAAGCCGTTCAAGGCGGCTTTCGCGCCGGGCGGCGCGATGTACAACTTCGTGTTCAACTTCGCCATCCCGCACCTCGACCAG‑CTGGTGATCAAGACCGCACCGATCGTTGCCGCGCCGACCGCCATGCCGGCGGTGTTCAAGCTCGACCCGATCTCCGCCACCGGCACCGCGAT  >  NC_002947/5384390‑5384672

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