Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,066,004	G→T	7.9%	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	7.9%	192.3 / 6.1	63	N1020K (AAC→AAA)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base T (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAG  >  NC_002947/4065856‑4066147                                                                                                                                                     |                                                                                                                                                gCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACgg                                                                                                                                                 >  4:10956/1‑149 (MQ=255) gCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACgg                                                                                                                                                 >  3:82472/1‑149 (MQ=255)            ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTg                                                                                                                                      <  7:29295/149‑1 (MQ=255)              ccACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgc                                                                                                                                    <  5:154703/149‑1 (MQ=255)              ccACGCACAACCCGACACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgc                                                                                                                                    >  5:331841/1‑149 (MQ=255)               cACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGaaa                                                                                                                                           >  5:73483/1‑141 (MQ=255)               cACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGaaa                                                                                                                                           <  6:73483/141‑1 (MQ=255)               cACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGgcg                                                                                                                                   >  8:79403/1‑149 (MQ=255)                      aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGTGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            <  7:234993/149‑1 (MQ=255)                      aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            >  7:324045/1‑149 (MQ=255)                      aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            >  8:39704/1‑149 (MQ=255)                      aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            <  4:82472/149‑1 (MQ=255)                         ccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                                          >  7:42574/1‑148 (MQ=255)                          cGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGtt                                                                                                                                               <  6:345402/126‑1 (MQ=255)                          cGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGtt                                                                                                                                               >  5:345402/1‑126 (MQ=255)                            gCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcagg                                                                                                                      <  8:115437/149‑1 (MQ=255)                                   gCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacaggt                                                                                                               <  3:10956/149‑1 (MQ=255)                                     aGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTaa                                                                                                             <  8:42574/149‑1 (MQ=255)                                            aGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGGCGGTACGCACGCCGCGACAACGGATGAAACGGTGGCGCAGGCGGTAGAGGTACAGGTAAACCACCg                                                                                                      >  6:13380/1‑149 (MQ=255)                                              ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCACCGTTTCAACCGTTGgcg                                                                                                                                   >  7:334401/1‑118 (MQ=255)                                              ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCACCGTTTCAACCGTTGgcg                                                                                                                                   <  8:334401/118‑1 (MQ=255)                                              ccacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGc                                                                                                    >  8:259443/1‑149 (MQ=255)                                                    ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccac                                                                                                        <  1:181196/139‑1 (MQ=255)                                                    ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAaccac                                                                                                        >  2:181196/1‑139 (MQ=255)                                                    ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                                              <  6:331841/149‑1 (MQ=255)                                                      acggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGtat                                                                                             >  6:9784/1‑148 (MQ=255)                                                        gggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                    <  1:312125/123‑1 (MQ=255)                                                        gggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggta                                                                                                                    >  2:312125/1‑123 (MQ=255)                                                         ggTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAg                                                                                         <  7:259443/149‑1 (MQ=255)                                                           tCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGtt                                                                                                                                               <  8:184353/93‑1 (MQ=255)                                                           tCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGtt                                                                                                                                               >  7:184353/1‑93 (MQ=255)                                                             aGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGTTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTg                                                                                     <  5:114776/149‑1 (MQ=255)                                                              gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            <  8:252461/109‑1 (MQ=255)                                                              gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCgg                                                                                                                            >  7:252461/1‑109 (MQ=255)                                                              gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                    <  7:306224/149‑1 (MQ=255)                                                                      aGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTTGTATAAAGGGTGAGGATCTgg                                                                            >  4:323658/1‑149 (MQ=255)                                                                       gggCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTTGTATAAAGGGTGAGGATCTGGc                                                                           <  3:323658/149‑1 (MQ=255)                                                                           aCGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGaa                                                                       >  2:112210/1‑149 (MQ=255)                                                                             ggTGAAGTGGGGTCTGGGCAAAAGTCATAGGGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAaca                                                                     >  1:29706/1‑149 (MQ=255)                                                                             ggTGAAGTGGGGTCTGGGCAAAAATCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcgt                                                                                               >  6:39114/1‑123 (MQ=255)                                                                             ggTGAAGTGGGGTCTGGGCAAAAATCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtcgt                                                                                               <  5:39114/123‑1 (MQ=255)                                                                                     ggggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCa                                                             >  2:270769/1‑149 (MQ=255)                                                                                      gggTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCa                                                             <  7:79403/148‑1 (MQ=255)                                                                                             ggCAAAAGTCATTGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgatgac                                                     <  2:29706/149‑2 (MQ=255)                                                                                             ggCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCgat                                                        <  7:39704/146‑1 (MQ=255)                                                                                              gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGatgatg                                                    <  8:324045/149‑1 (MQ=255)                                                                                               cAAAAGTCATAGCGGGCTTTCCAGAGCGGCGTCGCTACGCACTCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATgg                                                   <  4:58306/149‑1 (MQ=255)                                                                                                 aaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTc                                                 <  5:301027/149‑1 (MQ=255)                                                                                                     tCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCTGCGCCACCGTTTCAACCGTTgg                                                                                                                                     <  7:306854/61‑1 (MQ=21)                                                                                                     tCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCTGCGCCACCGTTTCAACCGTTgg                                                                                                                                     >  8:306854/1‑61 (MQ=21)                                                                                                       aTAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                    >  3:202954/1‑108 (MQ=255)                                                                                                       aTAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                                    <  4:202954/108‑1 (MQ=255)                                                                                                         aGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAg                                         >  8:16315/1‑149 (MQ=255)                                                                                                          gCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGc                                        <  5:9784/149‑1 (MQ=255)                                                                                                           cGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCg                                       <  7:16315/149‑1 (MQ=255)                                                                                                                      gagCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCa                              >  1:128278/1‑147 (MQ=255)                                                                                                                                     gcacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGcc             >  7:201082/1‑149 (MQ=255)                                                                                                                                     gcacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGcc             >  7:226196/1‑149 (MQ=255)                                                                                                                                      cacgcCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCt            <  4:314419/149‑1 (MQ=255)                                                                                                                                        cgccggGCAAACGGTGTAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGCGTATCTGGCTCAACACCAGGCCTCCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTgg          <  1:270769/149‑1 (MQ=255)                                                                                                                                         cctgcgcCACCGTTTCAACCGTTGGCGCAGGCAGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGc         >  1:43468/4‑149 (MQ=255)                                                                                                                                           cgcgGCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCATGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCt        <  8:226196/148‑1 (MQ=255)                                                                                                                                            gcgcCAACGGTTGAAACGGTGGCGGAGGCGGTCCAGGTAGAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTca      <  5:13380/149‑1 (MQ=255)                                                                                                                                                cAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagcag  >  3:250842/1‑149 (MQ=255)                                                                                                                                                 aaCGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTcagcag  >  7:143063/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCCGCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTTCAGCCTGGCTCAGCAG  >  NC_002947/4065856‑4066147

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