Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,066,010	A→C	7.0%	R1018R (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,010	0	A	C	7.0%	166.7 / 4.5	57	R1018R (CGT→CGG)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base C (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

AGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAAT  >  NC_002947/4065862‑4066157                                                                                                                                                      |                                                                                                                                                    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                                                                                                                              <  6:234502/149‑1 (MQ=255)                        aCCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                              >  5:234502/1‑149 (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                          <  3:162244/149‑1 (MQ=255)                                                                                                                            tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa                          <  2:365763/149‑1 (MQ=255)                                                                                                                            tACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCa                          <  4:219590/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)                                                                                                                                                      |                                                                                                                                                    AGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAGCAACGGCAAT  >  NC_002947/4065862‑4066157

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