Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,010,931	A→T	7.9%	Y44* (TAT→TAA)	PP_0871 ←	glycine betaine/carnitine/choline ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,010,931	0	A	T	7.9%	176.2 / 2.3	61	Y44* (TAT→TAA)	PP_0871	glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base T (2/3);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGC  >  NC_002947/1010784‑1011079                                                                                                                                                     |                                                                                                                                                      gaagaaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGa                                                                                                                                                       >  5:191276/1‑148 (MQ=255)  aagaaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGAt                                                                                                                                                      >  2:118393/1‑148 (MQ=255)     aaCTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCAATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAAcg                                                                                                                                                  >  8:28943/1‑149 (MQ=255)       cTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAAcgcg                                                                                                                                                <  6:34105/149‑1 (MQ=255)           aTGAAGCGGTCGGATTTTTCTACCCGGTGCGTCCTACTTAACGAAATGCCGGCGGGTATGCCAACGGCAAGCGCCGCCGCCATCTAAGTCCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATcgcgacga                                                                                                                                            <  1:99103/149‑5 (MQ=255)             gAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCg                                                                                                                                         >  1:187906/1‑149 (MQ=255)              aaGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATctt                                                                                                                                        <  2:69247/149‑3 (MQ=255)               aGCGCTCGGCTTTGTCGACCCGGCGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtg                                                                                                                                       >  5:228521/1‑149 (MQ=255)                gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      >  5:336287/1‑149 (MQ=255)                gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      <  7:30448/149‑1 (MQ=255)                gcgcTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgtgt                                                                                                                                      <  6:63230/149‑1 (MQ=255)                      ggCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGGAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGAt                                                                                                                                  <  2:39202/147‑1 (MQ=255)                        cTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCg                                                                                                                              <  7:28943/149‑1 (MQ=255)                        cTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCg                                                                                                                              <  5:343944/149‑1 (MQ=255)                           tGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATc                                                                                                                           >  3:281436/1‑149 (MQ=255)                            gTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACTGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGTTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATc                                                                                                                           <  1:203709/148‑1 (MQ=255)                              cGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATCttg                                                                                                                                       >  2:75361/1‑134 (MQ=255)                              cGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCG‑CGCGTTATCGCGACGATCttg                                                                                                                                       <  1:75361/134‑1 (MQ=255)                                aCCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTg                                                                                                                      <  1:118393/149‑1 (MQ=255)                                     gTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGAc                                                                                                                   >  1:178045/1‑147 (MQ=255)                                         ggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACcagcag                                                                                                             <  2:187906/149‑1 (MQ=255)                                         ggTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACcagcag                                                                                                             <  6:336287/149‑1 (MQ=255)                                             gACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGc                                                                                                          <  6:191276/148‑1 (MQ=255)                                               cTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAg                                                                                                       <  4:281436/149‑1 (MQ=255)                                                  aaCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAgg                                                                                                      <  1:300216/147‑1 (MQ=255)                                                   aCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                                                                                                    <  1:21547/148‑1 (MQ=255)                                                   aCGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGAt                                                                                                    <  2:103406/148‑1 (MQ=255)                                                         tGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGc                                                                                             <  6:228521/149‑1 (MQ=255)                                                                        ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgt                                                                                                                                        <  3:322621/91‑1 (MQ=255)                                                                        ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCgt                                                                                                                                        >  4:322621/1‑91 (MQ=255)                                                                        ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc                                                                              >  6:184047/1‑149 (MQ=255)                                                                        ccACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATcc                                                                               >  1:77276/1‑148 (MQ=255)                                                                         cACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATccc                                                                              >  7:64887/1‑148 (MQ=255)                                                                               cAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGcc                                                                       >  1:151757/1‑149 (MQ=255)                                                                                       ccgccATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCg                                                               <  5:215074/149‑1 (MQ=255)                                                                                        cgccATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCGCTTCCGAACAGCCCCTTTCCg                                                               >  8:70971/1‑148 (MQ=255)                                                                                                aaG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTa                                                      <  2:155692/149‑1 (MQ=255)                                                                                                  g‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGc                                                    <  5:184047/149‑1 (MQ=255)                                                                                                     cAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCca                                                  <  8:64887/149‑1 (MQ=255)                                                                                                           cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  4:308015/1‑149 (MQ=255)                                                                                                           cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            <  2:151757/149‑1 (MQ=255)                                                                                                           cAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  4:271262/1‑149 (MQ=255)                                                                                                                 cAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTcc                                      >  3:286201/1‑149 (MQ=255)                                                                                                                   ggTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCtt                                    <  2:178045/149‑1 (MQ=255)                                                                                                                     tccgcTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGGGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTtgt                                  <  7:70971/147‑1 (MQ=255)                                                                                                                     tGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTtgt                                  <  2:77276/149‑1 (MQ=255)                                                                                                                              aGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGc                         >  2:125512/1‑149 (MQ=255)                                                                                                                                  aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            >  5:240600/1‑126 (MQ=255)                                                                                                                                  aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTg                                            <  6:240600/126‑1 (MQ=255)                                                                                                                                  aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGa                                          >  5:176009/1‑128 (MQ=255)                                                                                                                                  aGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGa                                          <  6:176009/128‑1 (MQ=255)                                                                                                                                       ccAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACAct                >  8:69675/1‑149 (MQ=255)                                                                                                                                       ccAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACAct                >  4:129130/1‑149 (MQ=255)                                                                                                                                         aaGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACActtc              >  5:312859/1‑149 (MQ=255)                                                                                                                                           gATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACActtctt            <  4:286201/149‑1 (MQ=255)                                                                                                                                            aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           <  3:271262/149‑1 (MQ=255)                                                                                                                                            aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           >  6:157958/1‑149 (MQ=255)                                                                                                                                            aTCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           <  5:70405/149‑1 (MQ=255)                                                                                                                                            aTCGTCACGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTg           >  5:252869/1‑149 (MQ=255)                                                                                                                                               cgcgcgTTATCGCGACGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGt         >  3:210059/3‑149 (MQ=255)                                                                                                                                                 cgcgATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGc      >  8:251149/1‑149 (MQ=255)                                                                                                                                                   cgATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCg    <  1:125512/149‑1 (MQ=255)                                                                                                                                                     aTAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGc  <  5:157958/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GAAGAACTGCATGAAGCGCTCGGCTTTGTCGACCCGGTGCGGTCGACTTAACGCAATGCCGGCGGGGATGCCCACGGCCAGCGCCGCCGCCATCGAAG‑CCAGCACCAATACCAGGTGCGCTTGCAGGTAGAACCCAAGATCGTCGCGATAACGCG‑CGATCGTGTCGATGCCGATCCAGTGGACCAGCAGGGCCAGGATGACGAGCACGGCGGCCCATCCCAACAGCCCCTTTCCGTAGCGTTTAGCCACAGGTGGACTCCTTGTGTAGTCGGCGAGCACACTTCTTGGTTGCCGGC  >  NC_002947/1010784‑1011079

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