Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,066,004	G→T	6.4%	N1020K (AAC→AAA)	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,004	0	G	T	6.4%	193.3 / 3.0	62	N1020K (AAC→AAA)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCA  >  NC_002947/4065857‑4066146                                                                                                                                                     |                                                                                                                                               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)                                                                                                                                                     |                                                                                                                                               CACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCA  >  NC_002947/4065857‑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: ‑