Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,158,037	G→T	8.7%	R19L (CGG→CTG)	opuA →	glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,158,037	0	G	T	8.7%	198.4 / 4.4	64	R19L (CGG→CTG)	opuA	glycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base T (2/4);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGAAGTG  >  NC_002947/3157891‑3158182                                                                                                                                                   |                                                                                                                                                  aCCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGa                                                                                                                                                 >  7:110152/1‑149 (MQ=255)    ccATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCa                                                                                                                                              <  5:150038/149‑1 (MQ=255)         ccGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAgcg                                                                                                                                           >  6:116327/1‑147 (MQ=255)         ccGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCa                                                                                                                                         <  2:47993/149‑1 (MQ=255)          cgcgCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAg                                                                                                                                        <  7:114584/149‑1 (MQ=255)           gcgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgt                                                                                                                                       <  2:218064/149‑2 (MQ=255)            cgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg                                                                                                                                      <  8:110152/149‑3 (MQ=255)            cgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg                                                                                                                                      <  4:63930/149‑3 (MQ=255)              cTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGc                                                                                                                                    >  2:203408/1‑149 (MQ=255)                  aTTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCtt                                                                                                                                <  1:317444/149‑1 (MQ=255)                      cgcgGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGc                                                                                                                             <  4:119862/148‑1 (MQ=255)                          ggTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGt                                                                                                                        <  2:243874/149‑1 (MQ=255)                           gTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTa                                                                                                                       <  6:330548/149‑1 (MQ=255)                            ttttGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTAt                                                                                                                      <  2:172870/149‑1 (MQ=255)                            ttttGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTAGCTCTGTAt                                                                                                                      <  3:140769/149‑1 (MQ=255)                                 gCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGct                                                                                                                            >  6:86528/1‑138 (MQ=255)                                 gCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGct                                                                                                                            <  5:86528/138‑1 (MQ=255)                                               aGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCt                                                                                                   >  3:320141/1‑149 (MQ=255)                                                gCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTg                                                                                                  <  6:93468/149‑1 (MQ=255)                                                 ccTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTg                                                                                                  >  5:187493/1‑148 (MQ=255)                                                   tGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACg                                                                                                             >  3:37024/1‑135 (MQ=255)                                                   tGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACg                                                                                                             <  4:37024/135‑1 (MQ=255)                                                     cTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGaa                                                                                               <  1:203408/147‑1 (MQ=255)                                                      ttCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTgcgc                                                                                                    >  7:288531/1‑141 (MQ=255)                                                      ttCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTgcgc                                                                                                    <  8:288531/141‑1 (MQ=255)                                                          cTGCCCTGGAATCCATTGAACGCCAGGTGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTg                                                                                        <  7:45320/149‑1 (MQ=255)                                                               cTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGt                                                                                    >  6:34714/1‑148 (MQ=255)                                                                       cATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAAc                                                                            <  8:256457/148‑1 (MQ=255)                                                                       cATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACg                                                                           <  6:154900/149‑1 (MQ=255)                                                                        aTTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACgg                                                                          >  7:142543/1‑149 (MQ=255)                                                                          tGAACGCCA‑GGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTAc                                                                       >  8:80277/1‑149 (MQ=255)                                                                            aaCGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTAc                                                                       <  5:116327/148‑1 (MQ=255)                                                                            aaCGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTAc                                                                       >  6:117857/1‑148 (MQ=255)                                                                                  aGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGa                                                                <  4:320141/149‑1 (MQ=255)                                                                                       gCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCTGTCACCAAACGGTACGAGGGGAGCCgg                                                           <  8:201894/149‑1 (MQ=255)                                                                                       gCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCgg                                                           >  7:132676/1‑149 (MQ=255)                                                                                            gTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGa                                                      <  7:80277/149‑1 (MQ=255)                                                                                              gATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACt                                                    >  1:242093/1‑149 (MQ=255)                                                                                               aTCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTc                                                                                   <  7:255375/117‑1 (MQ=255)                                                                                               aTCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTc                                                                                   >  8:255375/1‑117 (MQ=255)                                                                                                   cttaTAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTCCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTccacc                                               <  3:208988/146‑1 (MQ=255)                                                                                                           tACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGAGATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGt                                       >  7:324873/1‑149 (MQ=255)                                                                                                            aCCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGtt                                      >  3:73953/1‑149 (MQ=255)                                                                                                              ccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa                                    <  6:187493/149‑1 (MQ=255)                                                                                                              ccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa                                    <  5:34714/149‑1 (MQ=255)                                                                                                                acaGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAg                                   >  6:252581/1‑148 (MQ=255)                                                                                                                     ttGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTg                             <  5:117857/149‑1 (MQ=255)                                                                                                                         cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAAcc                         >  7:215202/1‑149 (MQ=255)                                                                                                                                ccTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGt                  >  2:333835/1‑149 (MQ=255)                                                                                                                                 cTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa                                    >  7:81469/1‑130 (MQ=255)                                                                                                                                 cTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa                                    <  8:81469/130‑1 (MQ=255)                                                                                                                                  tttGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTgc                >  1:73916/1‑149 (MQ=255)                                                                                                                                     ggAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTgcggc             >  8:58032/1‑149 (MQ=255)                                                                                                                                      gAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTgcggcg            >  3:200377/1‑149 (MQ=255)                                                                                                                                      gAAACCACTGCGCGGATCAGCGCAGACGCACTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTgcggc             <  4:73953/148‑1 (MQ=255)                                                                                                                                       aaaCCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGt           >  1:147392/1‑149 (MQ=255)                                                                                                                                        atccACTGAGATGATCAGCGAATACGCCCTTTGCTCTGTATCGATTCACGATATTGCTCTTTAATTTCGTTCTGTCACCAAACTGTACTATGGGATCAGGTCAGACTCCACCAACGCGGTTGAGCAATTGAAACTCGCGGTGCGTCGTg          <  8:11824/147‑1 (MQ=255)                                                                                                                                        aaCCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTg          >  8:200872/1‑149 (MQ=255)                                                                                                                                          cgtcTGCGCTGATCCGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGCCTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGc        >  7:101212/4‑149 (MQ=255)                                                                                                                                          ccACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGtt                                      <  6:206646/119‑1 (MQ=255)                                                                                                                                          ccACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGtt                                      <  6:201607/119‑1 (MQ=255)                                                                                                                                          ccACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGtt                                      >  5:206646/1‑119 (MQ=255)                                                                                                                                          ccACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGtt                                      >  5:201607/1‑119 (MQ=255)                                                                                                                                             cTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGaa     <  2:73916/149‑1 (MQ=255)                                                                                                                                                cgcgGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGAAGt   <  3:67627/148‑1 (MQ=255)                                                                                                                                                cgcgGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGAAGt   <  4:200377/148‑1 (MQ=255)                                                                                                                                                cgcgGATCAGCGCAGAAGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCATACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGAAGtg  <  7:162394/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGAAGTG  >  NC_002947/3157891‑3158182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑