Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,066,001	A→C	6.2%	R1021R (CGT→CGG)	mdtC ←	multidrug transporter membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,066,001	0	A	C	6.2%	195.6 / 2.6	64	R1021R (CGT→CGG)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (2/2);  total (31/33)
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.49e-01

AGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCA  >  NC_002947/4065854‑4066146                                                                                                                                                     |                                                                                                                                                  aGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAAc                                                                                                                                                  >  1:322438/1‑149 (MQ=255)  ggCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACg                                                                                                                                                 <  3:192508/149‑1 (MQ=255)    cACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGACAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGt                                                                                                                                               >  3:366014/1‑149 (MQ=255)     aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                              <  1:103950/149‑1 (MQ=255)     aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                              <  4:163250/149‑1 (MQ=255)     aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                              <  7:193184/149‑1 (MQ=255)     aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                              <  1:14823/149‑1 (MQ=255)         aGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGaaa                                                                                                                                          >  2:54913/1‑149 (MQ=255)             aGGCCACGCTCAACCCGGCACCGCGCAGGCGCGTGACACGGCCACGGGTCTGTGCTTGCAGGGCACGGGGATGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAATCGGt                                                                                                                                      <  6:26326/149‑1 (MQ=255)             aGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGTCACGGTGAAGTGGGGTGTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGt                                                                                                                                      <  4:264128/149‑1 (MQ=255)              ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTg                                                                                                                                     <  8:147218/149‑1 (MQ=255)                  aCGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTAGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGgcg                                                                                                                                  >  5:390757/1‑148 (MQ=255)                    gCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAg                                                                                                                               >  5:169024/1‑149 (MQ=255)                     cacaACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCT‑GCGCCACCGTTTCAACCGTTGgcgcggc                                                                                                                              <  3:171773/149‑4 (MQ=255)                           ccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCTGCGCCACCGTTTCAACCGTTGgcgcg                                                                                                                                <  4:212223/142‑2 (MQ=255)                           ccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCTGCGCCACCGTTTCAACCGTTGgcgcg                                                                                                                                >  3:212223/1‑141 (MQ=255)                            cGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                         >  5:149758/1‑147 (MQ=255)                                aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                   <  6:390028/149‑1 (MQ=255)                                aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                   >  5:234502/1‑149 (MQ=255)                                aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                   <  6:234502/149‑1 (MQ=255)                                   gAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacag                                                                                                                <  8:123734/149‑1 (MQ=255)                                    aGCAGGCGCGAGGCACGGCCACGGGTCAGTGCTTGCAGGGCACTGTGAAGTGGGGTCTGGGCAAATGTCATAGCGGGTTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacagg                                                                                                               <  4:81623/149‑1 (MQ=255)                                            cgAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccac                                                                                                       >  4:82152/1‑149 (MQ=255)                                                    ggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcg                                                                                               <  2:96728/149‑1 (MQ=255)                                                     gccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAGAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcg                                                                                               >  7:230157/1‑148 (MQ=255)                                                      ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                                             >  2:66614/1‑149 (MQ=255)                                                      ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGACTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTGCAGGTACAGGTAAACCACCGGCGTCGta                                                                                             >  4:162244/1‑149 (MQ=255)                                                       cacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGtat                                                                                            <  3:234603/149‑1 (MQ=255)                                                          gggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGACAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATaaa                                                                                         <  4:366014/149‑1 (MQ=255)                                                                 tGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAg                                                                                  <  6:169024/149‑1 (MQ=255)                                                                    ttGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGtat                                                                                            <  6:3743/136‑1 (MQ=255)                                                                    ttGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGATTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGtat                                                                                            >  5:3743/1‑136 (MQ=255)                                                                        aGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTgg                                                                           <  2:322438/149‑1 (MQ=255)                                                                           gCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTg                                                                        <  2:109963/149‑1 (MQ=255)                                                                           gCACGGTGAAGT‑GGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTAGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCATGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGATGATCTGGCTg                                                                        >  3:219590/1‑148 (MQ=255)                                                                               ggTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAaca                                                                    >  1:365763/1‑149 (MQ=255)                                                                                        gggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGAC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcag                                                                                                                      <  2:23494/90‑1 (MQ=255)                                                                                        gggTCTGGGCAAAAGTCATAGCGGGGTTACCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcag                                                                                                                      >  1:23494/1‑90 (MQ=255)                                                                                             tGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCGACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgat                                                       >  3:162129/1‑148 (MQ=255)                                                                                             tGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgat                                                       >  8:277417/1‑148 (MQ=255)                                                                                               ggCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgatgat                                                    >  1:323651/1‑149 (MQ=255)                                                                                               ggCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgatga                                                     <  6:390757/148‑1 (MQ=255)                                                                                                 cAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTCGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATgg                                                  >  5:267134/1‑149 (MQ=255)                                                                                                  aaaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGt                                                 >  3:259461/1‑149 (MQ=255)                                                                                                   aaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTc                                                <  1:66614/149‑1 (MQ=255)                                                                                                      gTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGt                                                                                     <  6:74120/109‑1 (MQ=255)                                                                                                               gggTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGc                                    >  1:115142/1‑149 (MQ=255)                                                                                                               gggTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGc                                    >  1:120088/1‑149 (MQ=255)                                                                                                                ggTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCt                                   >  6:52778/1‑149 (MQ=255)                                                                                                                  tttCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAAGACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGa                                 <  8:230157/149‑1 (MQ=255)                                                                                                                  tttCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGa                                 <  3:82152/149‑1 (MQ=255)                                                                                                                      cAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCa                             >  7:321359/1‑149 (MQ=255)                                                                                                                      cAGAGCGGCGTCGGTACGAACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACgaa                             >  5:91850/1‑147 (MQ=255)                                                                                                                        gagCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCAtt                           <  6:149758/149‑1 (MQ=255)                                                                                                                          gcggcgTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTc                         >  3:130069/1‑149 (MQ=255)                                                                                                                                   gTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCtt                 >  6:235253/1‑148 (MQ=255)                                                                                                                                   gTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGACACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCtt                 >  6:134256/1‑148 (MQ=255)                                                                                                                                    tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa               <  3:162244/149‑1 (MQ=255)                                                                                                                                    tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa               <  4:219590/149‑1 (MQ=255)                                                                                                                                    tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa               <  2:365763/149‑1 (MQ=255)                                                                                                                                       gcacgcC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGcc            >  8:229135/1‑149 (MQ=255)                                                                                                                                            cctgcgcCACCGTTTCAACCGTTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGc        >  8:85061/4‑149 (MQ=255)                                                                                                                                                 gcCTACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagc   <  4:162129/149‑1 (MQ=255)                                                                                                                                                  ccAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagca  <  5:235253/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCA  >  NC_002947/4065854‑4066146

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