Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,590,016	G→T	5.6%	L179L (CTG→CTT)	nicP‑I →	porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,590,016	0	G	T	5.6%	241.8 / 3.5	78	L179L (CTG→CTT)	nicP‑I	porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (41/32);  new base T (2/2);  total (44/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 = 1.75e-01

TGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCT  >  NC_002947/3589883‑3590162                                                                                                                                      |                                                                                                                                                   tGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTTGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGc                                                                                                                                     >  7:269156/1‑149 (MQ=255) tGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGc                                                                                                                                     >  6:328163/1‑149 (MQ=255)  gCAACCGAACCTGCCGGTGCTCACATTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGATTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCAGGGCACTTGAAAAGCACCCACCTGAGTAACGAAGGTGGc                                                                                                                                     >  5:8683/1‑148 (MQ=255)   cAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGa                                                                                                                                   >  3:247297/1‑149 (MQ=255)    aaCCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGAt                                                                                                                                  >  5:44082/1‑149 (MQ=255)      ccGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATgg                                                                                                                                >  3:256523/1‑149 (MQ=255)       cGAACCTGCCGTTGCAGAGCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACTGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGc                                                                                                                               <  2:244959/149‑1 (MQ=255)       cGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAg                                                                                                                                          <  4:101758/138‑1 (MQ=255)       cGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAg                                                                                                                                          >  3:101758/1‑138 (MQ=255)                   tGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAAC                                                                                                                      <  6:18981‑M1/146‑1 (MQ=255)                   tGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAAC                                                                                                                      >  5:18981‑M1/1‑146 (MQ=255)                       gACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGc                                                                                                               >  8:198958/1‑149 (MQ=255)                       gACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGAGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGc                                                                                                               >  8:193704/1‑149 (MQ=255)                        aCCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATg                                                                                                                <  5:323481/147‑1 (MQ=255)                           ttCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATACACGGGCTGACTATGCAGGCCGGGCAATTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCgc                                                                                                           >  5:231364/1‑148 (MQ=255)                              aGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCa                                                                                                        >  8:173000/1‑149 (MQ=255)                                   tatTCGCCTGTTGCCACCCACCTATCAGTGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTTCGTTACGCAg                                                                                                                                          <  8:180725/110‑1 (MQ=255)                                   tatTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTTCGTTACGCAg                                                                                                                                          >  7:180725/1‑110 (MQ=255)                                   tatTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGc                                                                                                   >  2:206633/1‑149 (MQ=255)                                   tatTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGc                                                                                                   >  3:123204/1‑149 (MQ=255)                                      tCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGa                                                                                                >  5:190943/1‑149 (MQ=255)                                      tCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGa                                                                                                >  2:133190/1‑149 (MQ=255)                                        gCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTTGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGAt                                                                                               <  8:269156/148‑1 (MQ=255)                                        gCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTTACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGAt                                                                                               <  1:85473/148‑1 (MQ=255)                                        gCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATg                                                                                              >  5:65711/1‑149 (MQ=255)                                            gTTGCTACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGc                                                                                                                                                   >  2:163989/1‑92 (MQ=255)                                            gTTGCTACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGc                                                                                                                                                   <  1:163989/92‑1 (MQ=255)                                                  acccaccTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGcc                                                                                     >  1:151242/1‑148 (MQ=255)                                                    ccaccTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAgg                                                                                                                                         >  8:326954/1‑94 (MQ=255)                                                    ccaccTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAgg                                                                                                                                         <  7:326954/94‑1 (MQ=255)                                                      accTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAg                                                                                <  8:124432/149‑1 (MQ=255)                                                          aTCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGAc                                                                            <  4:3861‑M2/149‑2 (MQ=255)                                                                gCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGtcgctttcggcctgtcgcatcggcacatggccgagcatcgcgttc                                                                      <  1:216438‑M1/149‑46 (MQ=255)                                                                      gTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAt                                                                <  7:193704/149‑1 (MQ=255)                                                                      gTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAt                                                                <  7:198958/149‑1 (MQ=255)                                                                          gAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGc                                                                                      >  4:164181/1‑123 (MQ=255)                                                                          gAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGc                                                                                      <  3:164181/123‑1 (MQ=255)                                                                          gAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAGgtgg                                                            >  7:2827/1‑149 (MQ=255)                                                                                gtcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGtcgctttcggcctgtcgcatcggcacatggccgagcatcgcgttcat                                                                    <  5:197743‑M1/135‑48 (MQ=255)                                                                                gtcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGtcgctttcggcctgtcgcatcggcacatggccgagcatcgcgttcat                                                                    >  6:197743‑M1/1‑88 (MQ=255)                                                                                gtcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATgtgggtgg                                                        <  5:180628/147‑1 (MQ=255)                                                                                 tcgGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGa                                                     >  4:256025/1‑149 (MQ=255)                                                                                   ggAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGAc                                                    >  8:176078/1‑148 (MQ=255)                                                                                   ggAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGAc                                                    >  6:222443/1‑148 (MQ=255)                                                                                     tattCGAGGGGCTGTCTATTCAGGGCGGGCACTTGAAAATCACCCACCTGCTTAACGAAGGTGGAGATGGCAAGATGATCGCGATGCTCGGCCTTGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTAt                                                 <  6:231364/146‑1 (MQ=255)                                                                                     aaaTCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTAt                                                 <  6:190943/149‑1 (MQ=255)                                                                                        tCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCg                                              <  7:173000/149‑1 (MQ=255)                                                                                          gACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTAt                                                 <  7:329681/144‑1 (MQ=255)                                                                                          gACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTAt                                                 >  8:329681/1‑144 (MQ=255)                                                                                          gACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGtt                                            <  1:206633/149‑1 (MQ=255)                                                                                          gACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGt                                             <  4:247297/148‑1 (MQ=255)                                                                                             gggCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCaa                                         <  5:222443/149‑1 (MQ=255)                                                                                               gCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATg                                       >  7:201878/1‑149 (MQ=255)                                                                                                cTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGt                                                                                                                                        >  5:163585/1‑51 (MQ=255)                                                                                                cTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGt                                                                                                                                        <  6:163585/51‑1 (MQ=255)                                                                                                  gACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGcc                                     >  6:82284/1‑148 (MQ=255)                                                                                                   aCTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaa                                   >  8:193355/1‑149 (MQ=255)                                                                                                   aCTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGcca                                    <  1:133190/148‑1 (MQ=255)                                                                                                    cTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaa                                   <  3:3861/148‑1 (MQ=255)                                                                                                     tATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAtgt                                                              >  7:178112/1‑120 (MQ=255)                                                                                                     tATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTAtgt                                                              <  8:178112/120‑1 (MQ=255)                                                                                                     tATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaac                                  >  5:329502/1‑148 (MQ=255)                                                                                                       tGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGccaacca                                >  6:218912/1‑148 (MQ=255)                                                                                                         cAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAag                             >  5:16997/1‑149 (MQ=255)                                                                                                         cAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAag                             >  5:11843/1‑149 (MQ=255)                                                                                                         cAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAag                             <  4:256523/149‑1 (MQ=255)                                                                                                          aGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAagc                            <  5:218912/149‑1 (MQ=255)                                                                                                                 gCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGt                       >  4:20418/1‑147 (MQ=255)                                                                                                                   acatGAAAACCACCCACCTGCGTAACGAAGGTGGCGATGGCACGATGAACTCGATGCTCGGCCATGTGCCGATGCGACACGCCGAAATCGACGCGTTCATCTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGGAGCGTCAgg                   <  8:2827/146‑2 (MQ=255)                                                                                                                   aCTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGc                   >  7:217238/1‑149 (MQ=255)                                                                                                                   aCTAGAAAAGCACCCACCTACGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATCTGGGTGGGGACTATGCGTTCAATGCCAACCAAACCAGCGTCAGc                   >  8:121227/1‑149 (MQ=255)                                                                                                                     ttGAAAAGCACCCACCTGCGTAACGCAGATGGCGATGACAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGCAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCt                 >  8:54205/1‑149 (MQ=255)                                                                                                                        aaaaGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCttt               <  7:193355/148‑1 (MQ=255)                                                                                                                            gCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTAc          >  4:101331/1‑149 (MQ=255)                                                                                                                               gccaccTTCGTTACGCAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCa                     <  2:43689/134‑1 (MQ=255)                                                                                                                               gccaccTTCGTTACGCAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCa                     >  1:43689/2‑135 (MQ=255)                                                                                                                                ccaccTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTc      <  2:151242/149‑1 (MQ=255)                                                                                                                                    cTGCGTAACGAAGGTGGCGATGGCAAGATGNACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCt  >  1:28470/1‑149 (MQ=255)                                                                                                                                      |                                                                                                                                                   TGCAACCGAACCTGCCGGTGCTGACCTTCAGTGATATTCGCCTGTTGCCACCCACCTATCAGGGCGTGAGTTTGAGCTCGTCGGAAATCGACGGGCTGACTATGCAGGCCGGGCACTTGAAAAGCACCCACCTGCGTAACGAAGGTGGCGATGGCAAGATGAACGCGATGCTCGGCCATGTGCCGATGCGACAGGCCGAAAGCGACGCGTTCAACTATGTGGGTGGGGACTATGCGTTCAATGCCAACCAAAGCAGCGTCAGCCTTTGGTACGGTCAGCT  >  NC_002947/3589883‑3590162

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