Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,098,863	C→G	6.2%	V177L (GTG→CTG) ‡	acsA‑I ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,098,863	0	C	G	6.2%	204.6 / 3.8	64	V177L (GTG→CTG) ‡	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base G (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

CCACCTTGGATTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCAAGGT  >  NC_002947/5098716‑5099009                                                                                                                                                    |                                                                                                                                                   ccaccTTGGATTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACg                                                                                                                                                   >  6:280456/1‑149 (MQ=255)  caccTTGGATTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGt                                                                                                                                                  >  8:3491/1‑149 (MQ=255)     cTTGGATTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCAc                                                                                                                                               >  6:275301/1‑149 (MQ=255)          aTTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCgg                                                                                                                                          >  6:198837/1‑149 (MQ=255)              gCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGca                                                                                                                                      >  6:189703/1‑149 (MQ=255)                    gatgatACGGCCAGCCAGTGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATcc                                                                                                                                 <  2:257028/148‑1 (MQ=255)                    gatgatACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCt                                                                                                                                <  4:273971/149‑1 (MQ=255)                     atgatACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTg                                                                                                                               <  7:169357/149‑1 (MQ=255)                     atgatACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTg                                                                                                                               <  4:36343/149‑1 (MQ=255)                        atACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCCCCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGc                                                                                                                             <  5:280456/148‑1 (MQ=255)                          aCGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCAc                                                                                                                          >  3:65189/1‑149 (MQ=255)                               cagccagcGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAgg                                                                                                                     <  8:130639/149‑1 (MQ=255)                                 gccagcGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACGTGGTCAccat                                                                                                                                            >  2:78661/1‑122 (MQ=255)                                 gccagcGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACGTGGTCAccat                                                                                                                                            <  1:78661/124‑3 (MQ=255)                                    agcGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGtt                                                                                                                <  6:249181/149‑1 (MQ=255)                                     gcgcCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTg                                                                                                               <  5:254992/149‑1 (MQ=255)                                     gcgcCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTg                                                                                                               >  3:269510/1‑149 (MQ=255)                                                  gaAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACGTGGTCAccat                                                                                                                                            >  3:26778/1‑105 (MQ=255)                                                  gaAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACGTGGTCAccat                                                                                                                                            <  4:26778/107‑3 (MQ=255)                                                         ccGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCg                                                                                                            >  5:49768/1‑132 (MQ=255)                                                         ccGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCg                                                                                                            <  6:49768/132‑1 (MQ=255)                                                            aaCACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTg                                                                                        <  5:198837/149‑1 (MQ=255)                                                              caccacCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCCCGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGc                                                                                       >  7:271378/1‑148 (MQ=255)                                                                 caccGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACt                                                                                   >  6:280480/1‑149 (MQ=255)                                                                 caccGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACt                                                                                   >  8:271030/1‑149 (MQ=255)                                                                  accGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTc                                                                                  >  7:3745/1‑149 (MQ=255)                                                                   ccGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTcg                                                                                 <  7:3491/149‑1 (MQ=255)                                                                    cGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTcgc                                                                                <  5:189703/149‑1 (MQ=255)                                                                          aaTCGCACCAATGCGTGCACAGGCCAGCATGGTGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGt                                                                                         <  2:21191/134‑1 (MQ=255)                                                                          aaTCGCACCAATGCGTGCACAGGCCAGCATGGTGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGt                                                                                         >  1:21191/1‑134 (MQ=255)                                                                           aTCGCACCAATGCGTGCACAGGCCAGCTTGGCGACGACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGt                                                                         <  8:4631/149‑1 (MQ=255)                                                                            tCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACcgtc                                                                                                                                           >  6:169924/1‑80 (MQ=255)                                                                              gCACCAATGCTTGCACAGTCAAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACTTCAAAGCGGTGCACATCCTGGCAACGCATGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGa                                                                       <  1:285565/148‑1 (MQ=255)                                                                              gCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGa                                                                       <  4:65189/148‑1 (MQ=255)                                                                               cACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATg                                                                     >  2:192166/1‑149 (MQ=255)                                                                                    aTGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCg                                                                <  5:275301/149‑1 (MQ=255)                                                                                     tGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCgg                                                               <  7:271030/149‑1 (MQ=255)                                                                                        gTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTTGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGc                                                            >  4:62770/1‑149 (MQ=255)                                                                                                  cAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAgg                                                    >  7:131705/1‑147 (MQ=255)                                                                                                  cAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAgg                                                    <  8:131705/147‑1 (MQ=255)                                                                                                  cAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGgt                                                  >  4:53362/1‑149 (MQ=255)                                                                                                   aGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGgtc                                                 <  8:3745/149‑1 (MQ=255)                                                                                                      aTGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGgtcgtc                                              <  5:280480/149‑1 (MQ=255)                                                                                                        ggCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCCAGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCa                                             >  2:212427/1‑148 (MQ=255)                                                                                                              caccGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTccc                                      <  8:186261/149‑1 (MQ=255)                                                                                                               accGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGaa                                                      <  1:265060/132‑1 (MQ=255)                                                                                                               accGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGaa                                                      >  2:265060/1‑132 (MQ=255)                                                                                                                ccGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTcc                                       >  7:72115/1‑146 (MQ=255)                                                                                                                  ccTTCGGCGATCATCGGCATATAGCTGGTGACCACGTCTCCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTAGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAgat                                  <  4:269510/148‑1 (MQ=255)                                                                                                                   cTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAgat                                  >  7:38297/1‑148 (MQ=255)                                                                                                                    ttCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAgatga                                >  4:138010/1‑149 (MQ=255)                                                                                                                       ggggATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGatgatg                              <  1:192166/148‑1 (MQ=255)                                                                                                                       ggggATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATgg                             >  7:86138/1‑149 (MQ=255)                                                                                                                                  gCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCa                  >  8:310108/1‑149 (MQ=255)                                                                                                                                       taGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACcgcg             >  8:56290/1‑149 (MQ=255)                                                                                                                                         gATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGct           <  3:53362/149‑1 (MQ=255)                                                                                                                                          aTGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGctc          >  3:236250/1‑149 (MQ=255)                                                                                                                                            ggTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCtt        >  1:209809/1‑149 (MQ=255)                                                                                                                                             gTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTc       <  7:56290/149‑1 (MQ=255)                                                                                                                                              tGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTTGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCa      <  3:62770/149‑1 (MQ=255)                                                                                                                                              tGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCa      <  7:310108/149‑1 (MQ=255)                                                                                                                                              tGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCATc       >  2:173579/1‑148 (MQ=255)                                                                                                                                                 ccACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCAAgg   >  4:259170/1‑149 (MQ=255)                                                                                                                                                  cACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCAAGGt  >  2:328985/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCACCTTGGATTTGCAGTCGATGATACGGCCAGCCAGCGCCTCGGGCGAGAAGCCACCGAACACCACCGAGTGAATCGCACCAATGCGTGCACAGGCCAGCATGGCGACCACCGCTTCGGGGATCATCGGCATATAGATGGTGACCACGTCACCGCGGTGCACATCCTGGCCACGCAGGGCGTTGGCGAACTTGCAGACCTGCTCGTGCAACTCGCGGTAAGTGATGTTGCGGTGCTCGGAAGGGTCGTCACCTTCCCAGATGATGGCCAACTGGTCACCGCGCTCTTCAAGGT  >  NC_002947/5098716‑5099009

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