Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,066,013	G→T	6.6%	H1017Q (CAC→CAA)	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,013	0	G	T	6.6%	186.1 / 3.9	61	H1017Q (CAC→CAA)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base T (2/2);  total (29/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 = 1.65e-01

AGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAATGCG  >  NC_002947/4065866‑4066160                                                                                                                                                     |                                                                                                                                                    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                                                                                                                                              >  3:212223/1‑141 (MQ=255)               ccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCTGCGCCACCGTTTCAACCGTTGgcgcg                                                                                                                                              <  4:212223/142‑2 (MQ=255)                cGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                                       >  5:149758/1‑147 (MQ=255)                    aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                                 >  5:234502/1‑149 (MQ=255)                    aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                                 <  6:390028/149‑1 (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:120088/1‑149 (MQ=255)                                                                                                   gggTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGc                                                  >  1:115142/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                             <  4:219590/149‑1 (MQ=255)                                                                                                                        tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa                             <  3:162244/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)                                                                                                                                             ttGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACgg         <  6:91850/149‑1 (MQ=255)                                                                                                                                                aaaCGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACAGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCaa      >  5:372976/1‑149 (MQ=255)                                                                                                                                                 aaCGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAAt     <  2:323651/149‑1 (MQ=255)                                                                                                                                                   cGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAATgc   <  2:120088/149‑1 (MQ=255)                                                                                                                                                   cGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAATgc   <  2:115142/149‑1 (MQ=255)                                                                                                                                                    ggTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAATgcg  <  5:52778/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAATGCG  >  NC_002947/4065866‑4066160

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