Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,235,310	C→G	6.6%	S215W (TCG→TGG)	fecA →	outer membrane ferric citrate porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,235,310	0	C	G	6.6%	246.6 / 6.6	76	S215W (TCG→TGG)	fecA	outer membrane ferric citrate porin
Reads supporting (aligned to +/- strand):  ref base C (30/41);  new base G (2/3);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

ACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGGCCGCCAT  >  NC_002947/5235163‑5235457                                                                                                                                                    |                                                                                                                                                     aCGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCg                                                                                                                                                     <  5:247602/149‑1 (MQ=255) aCGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCg                                                                                                                                                     <  7:226160/149‑1 (MQ=255)   gCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCt                                                                                                                                                   <  7:278389/149‑1 (MQ=255)       ccGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGcc                                                                                                                                               <  7:145199/149‑1 (MQ=255)       ccGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCGCGACCGTGCGGATTGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCGACAGCTGCGGCTGGcc                                                                                                                                               <  3:263852/149‑1 (MQ=255)        cGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCGACAGCTGCGGCTGGccg                                                                                                                                              <  7:201600/149‑2 (MQ=255)        cGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGccc                                                                                                                                              <  8:154266/149‑1 (MQ=255)         gACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGcccc                                                                                                                                             >  3:290974/1‑149 (MQ=255)            aacaacGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCGACAGCTGCGGCTGGccgtacg                                                                                                                                          <  7:154681/149‑6 (MQ=255)               aacGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAg                                                                                                                                       <  1:16355/149‑1 (MQ=255)                acGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGc                                                                                                                                      <  6:306176/149‑1 (MQ=255)                    tACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATg‑g                                                                                                                                 <  5:62874/149‑1 (MQ=255)                    tACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATg‑g                                                                                                                                 <  4:290974/149‑1 (MQ=255)                           aGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACAt                                                                                                                          <  2:153742/149‑1 (MQ=255)                               aTGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGAc                                                                                                                      >  2:83540/1‑149 (MQ=255)                                     tGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCa                                                                                                                   >  2:3455/1‑146 (MQ=255)                                     tGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCa                                                                                                                   <  1:3455/146‑1 (MQ=255)                                       gcgcTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGCAATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTAAGCATG‑GGACACATCGAAGCAGTGGa                                                                                                              >  6:264578/1‑149 (MQ=255)                                         gcTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACg                                                                                                            <  4:32706/149‑1 (MQ=255)                                            gAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAg                                                                                                         >  1:315668/1‑149 (MQ=255)                                             aaCTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCaa                                                                                                                             >  4:45971/1‑128 (MQ=255)                                             aaCTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCaa                                                                                                                             <  3:45971/128‑1 (MQ=255)                                             aaCTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGt                                                                                                        <  8:274338/149‑1 (MQ=255)                                             aaCTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGt                                                                                                        <  6:194222/149‑1 (MQ=255)                                               cTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTgc                                                                                                      >  4:316450/1‑149 (MQ=255)                                               cTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTgc                                                                                                      >  4:311255/1‑149 (MQ=255)                                                ttCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGACGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTgcg                                                                                                     >  4:286422/1‑149 (MQ=255)                                                   ggCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCgg                                                                                                  >  3:188403/1‑149 (MQ=255)                                                     cATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTg                                                                                                <  2:222041/149‑1 (MQ=255)                                                     cATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTg                                                                                                >  2:244275/1‑149 (MQ=255)                                                      aTCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTgg                                                                                               <  1:331917/149‑1 (MQ=255)                                                      aTCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTgg                                                                                               <  2:315668/149‑1 (MQ=255)                                                      aTCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTCAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTgg                                                                                               <  5:264578/149‑1 (MQ=255)                                                         cgcgGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTggcgg                                                                                            <  3:44015/149‑1 (MQ=255)                                                             gCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGcgcg                                                                                        <  8:48751/149‑1 (MQ=255)                                                             gCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGcgcg                                                                                        >  8:311670/1‑149 (MQ=255)                                                                  aaCCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTgcg                                                                                   >  6:48911/1‑149 (MQ=255)                                                                  aaCCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTg                                                                                     <  2:209109/147‑1 (MQ=255)                                                                  aaCCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTg                                                                                     <  7:201010/147‑1 (MQ=255)                                                                   aCCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATGCGGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTgcg                                                                                   <  1:83540/149‑1 (MQ=255)                                                                     acgcgCCTGGCATCGCGCTCGACAGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTa                                                                                <  6:1301/148‑1 (MQ=255)                                                                          ccTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTggc                                                                                              <  4:49794/130‑1 (MQ=255)                                                                          ccTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTggc                                                                                              >  3:49794/1‑130 (MQ=255)                                                                                   gcgcTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAAc                                                                   <  2:80416/148‑1 (MQ=255)                                                                                              tGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGaa                                                                    <  7:63371/136‑1 (MQ=255)                                                                                              tGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGaa                                                                    >  8:63371/1‑136 (MQ=255)                                                                                                 tGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGt                                                     >  2:297761/1‑148 (MQ=255)                                                                                                 tGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGt                                                     >  7:190335/1‑148 (MQ=255)                                                                                                   atggatggCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCaa                                                  >  7:313845/1‑149 (MQ=255)                                                                                                       atggCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCCTCGTCCACTTc                                              >  7:322272/1‑149 (MQ=255)                                                                                                       atggCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTc                                              <  3:286422/149‑1 (MQ=255)                                                                                                           cATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTgc                                                                                    <  8:49204/107‑1 (MQ=255)                                                                                                           cATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTgc                                                                                    >  7:49204/1‑107 (MQ=255)                                                                                                            aTCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCAc                                         <  1:297761/149‑1 (MQ=255)                                                                                                                 ggTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCggg                                    <  1:244275/149‑1 (MQ=255)                                                                                                                    gCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCggcggc                                                          <  4:149814/124‑1 (MQ=255)                                                                                                                    gCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGACCGCAGAACGTCggcggc                                                          >  3:149814/1‑124 (MQ=255)                                                                                                                          cgcgCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCt                            <  8:190335/148‑1 (MQ=255)                                                                                                                             gcCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGaa                         >  1:199906/1‑148 (MQ=255)                                                                                                                               cGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAg                      <  2:199906/149‑1 (MQ=255)                                                                                                                                gTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATc                                                       <  3:73827/115‑1 (MQ=255)                                                                                                                                gTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATc                                                       >  4:73827/1‑115 (MQ=255)                                                                                                                                gTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAgg                     >  6:186746/1‑149 (MQ=255)                                                                                                                                   gggCCAGCGACAGCTGCGGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCa                  >  4:173046/2‑149 (MQ=255)                                                                                                                                    ggCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCAc                 <  8:286901/149‑1 (MQ=255)                                                                                                                                      ccAGCCGCAGCTGTGGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATc                                                       >  7:324977/1‑109 (MQ=255)                                                                                                                                      ccAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATc                                                       <  8:324977/109‑1 (MQ=255)                                                                                                                                      ccAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGt               >  2:346300/1‑149 (MQ=255)                                                                                                                                      ccAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGt               >  2:33942/1‑149 (MQ=255)                                                                                                                                        aGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTc             <  4:188403/149‑1 (MQ=255)                                                                                                                                          ccGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCaa           <  5:48911/149‑1 (MQ=255)                                                                                                                                          ccGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCaa           <  7:311670/149‑1 (MQ=255)                                                                                                                                             cAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGgc        >  7:77133/1‑149 (MQ=255)                                                                                                                                                cTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGgccgc     >  5:242761/1‑149 (MQ=255)                                                                                                                                                  gTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGGCCGCCa   >  2:60674/1‑149 (MQ=255)                                                                                                                                                   tCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGGCCGCCAt  >  6:281568/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACGCCCCCGACAACAACGGTACCGGCAGCCATGACCTGGCGCTGAACTTCGGCATCCGCGGCCTCAACCCGCGCCTGGCATCGCGCTCGACCGTGCTGATGGATGGCATCCCGGTGCCCTTCGCGCCGTACGGCCAGCCGCAGCTGTCGCTGGCCCCGCTGAGCATG‑GGCAACATGGACGCAGTGGACGTAGTGCGCGGTGGCGGCGCGGTGCGCTACGGCCCGCAGAACGTCGGCGGCATCGTCAACTTCGTCACCCGGGCGATCCCTGAACAGGCCACGTTCAAGGCCGCCAT  >  NC_002947/5235163‑5235457

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