Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,921,094	A→T	5.8%	E501D (GAA→GAT)	PP_4330 →	toxin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,921,094	0	A	T	5.8%	207.8 / 2.6	69	E501D (GAA→GAT)	PP_4330	toxin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

TTCAAAAGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGAAACCGTCAA  >  NC_002947/4920951‑4921240                                                                                                                                                  |                                                                                                                                                    ttCAAAAGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGCGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGCCGAAGGC‑Ga                                                                                                                                               >  7:46182/1‑149 (MQ=255)   cAAAAGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCACGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GAc                                                                                                                                              <  5:335019/148‑1 (MQ=255)    aaaaGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGc                                                                                                                                                  <  5:169762/144‑1 (MQ=255)    aaaaGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGc                                                                                                                                                  >  6:169762/1‑144 (MQ=255)    aaaaGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GAccc                                                                                                                                            >  1:19032/1‑149 (MQ=255)       aGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTTTCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTGGACGAAGGC‑GACCCGGc                                                                                                                                         >  3:43238/1‑149 (MQ=255)           aCGGTTTTGTCATGCCGCGCGACGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGtt                                                                                                                                     >  8:157383/1‑149 (MQ=255)             ggTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGc                                                                                                                                   <  3:258698/149‑1 (MQ=255)                  tGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGgc                                                                                                                               <  5:50617/148‑1 (MQ=255)                   gTCATGCCGCGCGAGGTCAACGGGCAGTCTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGgcg                                                                                                                              >  1:15179/1‑149 (MQ=255)                      aTGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCA‑CACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGt                                                                                                                          >  5:328488/1‑148 (MQ=255)                           gcgcgAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATAc                                                                                                                     >  1:96638/1‑149 (MQ=255)                               gAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCtt                                                                                                                 <  4:202229/149‑1 (MQ=255)                                 ggTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGc                                                                                                               <  5:174341/149‑1 (MQ=255)                                  gTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGc                                                                                                               >  6:295913/1‑148 (MQ=255)                                   tCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCca                                                                                                             <  6:28198/149‑1 (MQ=255)                                    cAACGGGCAA‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTg                                                                                                                <  4:76535/145‑1 (MQ=255)                                    cAACGGGCAA‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTg                                                                                                                >  3:76535/1‑145 (MQ=255)                                     aaCGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCcaca                                                                                                           >  8:9040/1‑149 (MQ=255)                                     aaCGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGATTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTCATACGCTTGCcaca                                                                                                           >  4:95625/1‑149 (MQ=255)                                       cGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTg                                                                                                                         >  3:307804/1‑133 (MQ=255)                                       cGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTg                                                                                                                         <  4:307804/133‑1 (MQ=255)                                          gCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTAc                                                                                                       >  5:28830/1‑148 (MQ=255)                                           cAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGa                                                                                                     <  3:305758/149‑1 (MQ=255)                                             g‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATg                                                                                                   <  2:96638/149‑1 (MQ=255)                                               cTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATgg                                                                                                  >  3:28709/1‑149 (MQ=255)                                                tCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACAACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGc                                                                                                 >  1:101879/1‑149 (MQ=255)                                                 cGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCg                                                                                                <  7:9040/149‑1 (MQ=255)                                                   gCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGaa                                                                                              <  7:157383/149‑1 (MQ=255)                                                     tataCGCTAAGCCCGACCAGTCAACCACGGCTTCCCGAGCCAGAAGGCCATGAGTACTTCCACA‑ACCACGAGCCGATTACACCGCTTGACGAAGCC‑GACCCTGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGAAGACGAACt                                                                                            >  4:29661/1‑149 (MQ=255)                                                     tATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACt                                                                                            >  4:213324/1‑149 (MQ=255)                                                     tATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGAGAACt                                                                                            >  3:322018/1‑149 (MQ=255)                                                       tCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTg                                                                                           <  5:12093/148‑1 (MQ=255)                                                           cTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGATGGCAGACCCGGCTGTTGCCGGCGCGGTGATACGCtt                                                                                                                 <  1:39226/123‑1 (MQ=255)                                                           cTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGATGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCtt                                                                                                                 >  4:249413/1‑122 (MQ=255)                                                           cTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGATGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCtt                                                                                                                 <  3:249413/122‑1 (MQ=255)                                                           cTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGATGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCtt                                                                                                                 >  2:39226/1‑122 (MQ=255)                                                               gCCCGACCTGTCCACCAGGGCTTCCCGAGCCAGAAGGGGCGGAGTACTTCCACC‑ACCACGAGCGGATTACCCCGTTTGACGAAGGC‑GCCCCGGCTTTTGCCGGCGCGGTGATAGGCTTGCCACACTACTATGGCGAACTGgaggcccgg                                                                                  <  6:82152/149‑8 (MQ=255)                                                                 ccGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGa                                                                                <  2:19032/149‑1 (MQ=255)                                                                  cGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGaa                                                                               <  7:56550/149‑1 (MQ=255)                                                                  cGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGaa                                                                               <  6:28830/149‑1 (MQ=255)                                                                    aCCAGTCCA‑CACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATgc                                                                            <  6:328488/149‑1 (MQ=255)                                                                     ccAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATg                                                                             >  6:129539/1‑148 (MQ=255)                                                                      cAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATgca                                                                           <  4:28709/149‑1 (MQ=255)                                                                           caccacGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCAt                                                                       <  3:213324/148‑1 (MQ=255)                                                                                      cccGAGCCAGAAGGCTCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGTC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCAACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACa                                                            <  3:29661/148‑1 (MQ=255)                                                                                       ccGAGCCAGAAGGCGCGGAGTACTTCCACA‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAAt                                                          >  6:211148/1‑149 (MQ=255)                                                                                        cGAGCCAGAAGGCGCGGAGTACTTCCACCGACCCCGAGCCGATTACACCGATTGACGAAGGC‑TACCCGGATGTTGCCGGCGCGGTGATAGCCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAAACTAAAAt                                                          >  1:192797/1‑149 (MQ=255)                                                                                            ccAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTgg                                                                                 >  2:322223/1‑121 (MQ=255)                                                                                                  ggCGCGAAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTgg                                                >  8:139511/1‑148 (MQ=255)                                                                                                     gcgGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTa                                            >  1:227519/1‑149 (MQ=255)                                                                                                      cgGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTAt                                           >  5:17190/1‑149 (MQ=255)                                                                                                        gAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGa                                         <  2:15179/149‑1 (MQ=255)                                                                                                              ttCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTg                                     >  8:244224/1‑147 (MQ=255)                                                                                                                ccacc‑accaCGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGa                                   <  3:95625/147‑1 (MQ=255)                                                                                                                 cacc‑accacGAGCCGATTACACCGCTTGACGAAGGC‑GCCCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACCCTACGATGCCGCACTGGACGACCTGGAATGCAGCATCGCAACCTACAATGCACGCGAGGTCGATGACCTGGAGCa                                >  7:237431/1‑149 (MQ=255)                                                                                                                   cc‑accacGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCacg                              <  6:17190/149‑1 (MQ=255)                                                                                                                    c‑accacGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTa                                                              >  8:151287/1‑116 (MQ=255)                                                                                                                    c‑accacGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTa                                                              <  7:151287/116‑1 (MQ=255)                                                                                                                              cctaTTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCGTGCCACACTACGATGGCGAACTGTACGAACTGGAATGCAGCATCGAAACCTACAATGCAAGCGGTGTCTATGACCTGGAGCACGACGCTTtcat                    <  4:43238/146‑1 (MQ=255)                                                                                                                                gATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTtcatc                   <  8:100153/148‑1 (MQ=255)                                                                                                                                gATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTtcatc                   >  7:100153/1‑148 (MQ=255)                                                                                                                                     cacCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCaa              <  5:295913/148‑1 (MQ=255)                                                                                                                                        cGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGa           >  3:106219/1‑148 (MQ=255)                                                                                                                                          cTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAg                                  <  4:248641/123‑1 (MQ=255)                                                                                                                                          cTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAg                                  >  3:248641/1‑123 (MQ=255)                                                                                                                                              aCGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGAAACCGTc    <  7:244224/149‑1 (MQ=255)                                                                                                                                               cGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGAAACCGTc    <  7:139511/148‑1 (MQ=255)                                                                                                                                                gAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGAAACCGTCaa  >  2:164579/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TTCAAAAGCGACGGTTTTGTCATGCCGCGCGAGGTCAACGGGCAG‑CTCGGCTATCCGCTAAGCCCGACCAGTCCACCACGGCTTCCCGAGCCAGAAGGCGCGGAGTACTTCCACC‑ACCACGAGCCGATTACACCGCTTGACGAAGGC‑GACCCGGCTGTTGCCGGCGCGGTGATACGCTTGCCACACTACGATGGCGAACTGGACGACCTGGAATGCAGCATCGAAACCTACAATGCACGCGTGGTCTATGACCTGGAGCACGACGCTTTCATCATCAAGGAAACCGTCAA  >  NC_002947/4920951‑4921240

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