Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,820,205	A→T	6.5%	L3662Q (CTG→CAG)	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,820,205	0	A	T	6.5%	182.3 / 3.1	61	L3662Q (CTG→CAG)	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (2/2);  total (28/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 = 5.55e-01

GCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGCTGG  >  NC_002947/4820058‑4820352                                                                                                                                                    |                                                                                                                                                    gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                                    <  5:325695/149‑1 (MQ=255) gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                                    >  1:59556/1‑149 (MQ=255) gcTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCcagcagcag                                                                                                                                                    <  4:223347/149‑1 (MQ=255)   ttCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgc                                                                                                                                                  >  4:362597/1‑149 (MQ=255)      gAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgcgc                                                                                                                                                <  7:157545/148‑1 (MQ=255)      gAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGgcgc                                                                                                                                                <  2:193524/148‑1 (MQ=255)       aaCTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                              <  3:75349/149‑2 (MQ=255)       aaCTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCtt                                                                                                                                              >  3:331567/1‑149 (MQ=255)         cTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGa                                                                                                                                            <  5:62328/149‑1 (MQ=255)            cgcgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGaa                                                                                                                                           <  3:136863/147‑1 (MQ=255)            cgcgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAAtt                                                                                                                                         >  4:206889/1‑149 (MQ=255)              cgATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCg                                                                                                                                       <  2:184404/149‑1 (MQ=255)               gATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCgg                                                                                                                                      >  5:157734/1‑149 (MQ=255)                     gCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAACCCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                              >  3:229984/1‑134 (MQ=255)                     gCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCTGCTGCTGGCGCtg                                                                                                                                              <  4:229984/135‑2 (MQ=255)                       tctGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAgg                                                                                                                              >  5:242594/1‑149 (MQ=255)                         tGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGACcagcgctac                                                                                                                            >  4:356407/1‑144 (MQ=255)                            tCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTc                                                                                                                         >  8:45843/1‑149 (MQ=255)                              cAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGa                                                                                                                       >  5:133659/1‑149 (MQ=255)                              cAGCGGATAAACCCGCTCGGTCCGGGTGCAGGCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCGCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGa                                                                                                                       <  8:75015/149‑1 (MQ=255)                                   gATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTc                                                                                                                  <  6:242594/149‑1 (MQ=255)                                      aaaCCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCggg                                                                                                               >  1:322578/1‑149 (MQ=255)                                        aCCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCggg                                                                                                               >  5:152712/1‑147 (MQ=255)                                            gCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCg                                                                                                         <  3:206889/149‑1 (MQ=255)                                               cGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAg                                                                                                      <  8:222688/149‑1 (MQ=255)                                                      gTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTGACCTCCCCTTCAATACCTTGTAACAGCGCCAGCAGCAGGCGGTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTGGGAGTAGCGCt                                                                                               >  8:140204/1‑149 (MQ=255)                                                       ttGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTTGGGTAGTCGAAGTAGCGCTg                                                                                              <  8:360711/149‑1 (MQ=255)                                                             tcctccAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGaa                                                                                        <  4:371464/149‑1 (MQ=255)                                                               ctccAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAAc                                                                                       <  8:113153/148‑1 (MQ=255)                                                                  cagcagGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACa                                                                                   <  7:45843/149‑1 (MQ=255)                                                                     cagGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                                >  7:299283/1‑149 (MQ=255)                                                                        gAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                                >  1:58544/1‑146 (MQ=255)                                                                        gAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGt                                                                                <  2:58544/146‑1 (MQ=255)                                                                             cTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGcc                                                                         >  2:370495/1‑148 (MQ=255)                                                                             cTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTCGTCGAACGACAGGTGCAGGcc                                                                         >  5:90364/1‑148 (MQ=255)                                                                              ttGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGcc                                                                         >  3:120514/1‑147 (MQ=255)                                                                                gCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAg                                                                      <  2:59556/148‑1 (MQ=255)                                                                                 cTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAgg                                                                     <  3:231590/148‑1 (MQ=255)                                                                                  tCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGt                                                                    >  6:9640/1‑148 (MQ=255)                                                                                      gggCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcatca                                                               <  8:299283/149‑1 (MQ=255)                                                                                       ggCGCTCAGCAGCGGTAGCTCGCTCACCTGCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcatca                                                               <  3:87054/148‑1 (MQ=255)                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCgg                                                          <  2:322578/149‑1 (MQ=255)                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCgg                                                          <  1:370495/149‑1 (MQ=255)                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCg                                                           <  2:88763/148‑1 (MQ=255)                                                                                           cTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCg                                                           >  1:88763/1‑148 (MQ=255)                                                                                             cagcagCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGCCGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCggg                                                         >  6:59913/1‑148 (MQ=255)                                                                                             cagcagCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCggg                                                         >  5:80309/1‑148 (MQ=255)                                                                                               gcagcGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTaa                                                      <  5:268544/149‑1 (MQ=255)                                                                                                     gTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACAcc                                                <  7:140204/149‑1 (MQ=255)                                                                                                     gTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACAcc                                                >  6:169927/1‑149 (MQ=255)                                                                                                          tggctcACCTCCCCTTCAACACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAgg                                                                           <  8:166422/115‑1 (MQ=255)                                                                                                              tcACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGc                                       >  5:307800/1‑149 (MQ=255)                                                                                                                   tCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGATTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCggg                                    >  8:360394/1‑147 (MQ=255)                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                              <  6:80309/149‑1 (MQ=255)                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                              <  4:286898/149‑1 (MQ=255)                                                                                                                       cTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTg                              >  7:297137/1‑149 (MQ=255)                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCTGTCAGCGGGAAGTTGgt                            <  5:9640/148‑1 (MQ=255)                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTgg                             <  6:152712/147‑1 (MQ=255)                                                                                                                          cAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGgt                            <  7:360394/148‑1 (MQ=255)                                                                                                                            aaaCCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGtcat                                                                 <  6:330065/109‑1 (MQ=255)                                                                                                                             aaCCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTgg                        <  6:90364/149‑1 (MQ=255)                                                                                                                                cTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTa                     >  3:196694/1‑149 (MQ=255)                                                                                                                                 ttGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTAc                    >  7:2053/1‑149 (MQ=255)                                                                                                                                       aaGCGCCTGCTGCTGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc               >  7:63026/2‑148 (MQ=255)                                                                                                                                              agcagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGa       <  4:196694/149‑1 (MQ=255)                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc               <  7:78960/139‑1 (MQ=255)                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGccc               >  8:78960/1‑139 (MQ=255)                                                                                                                                                cagcagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGc     >  8:219875/1‑149 (MQ=255)                                                                                                                                                   cagGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGCTgg  <  4:331567/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTTCGAACTGCGCGATGTAGCTCTGCTCCAGCGGATAAACCCGCTCGGTCCGGTTGCAGTCCTCCAGCAGGAAGCCTTGCTCCTGGGCGCTCAGCAGCGGTAGCTCGCTCACCTCCCCTTCAAAACCTTGTACCAGCGCCAGCAGCAGGCGCTTGAATTCGGCCAGCAGGCGCTCGACGGTCGGGTAGTCGAAGTAGCGCTGGTCGAACGACAGGTGCAGGCCCAGGTCATCACCCGGGTAACACACCGCCGTCAGCGGGAAGTTGGTGTGGGTACGGCCCGAGTCGGAGCTGG  >  NC_002947/4820058‑4820352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑