Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,943	G→T	5.1%	L59M (CTG→ATG)	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,943	0	G	T	5.1%	290.5 / 2.5	89	L59M (CTG→ATG)	PP_4681	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/45);  new base T (2/3);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAAT  >  NC_002947/5319795‑5320090                                                                                                                                                        |                                                                                                                                                     gCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAg                                                                                                                                                      >  6:267909/1‑149 (MQ=255)    aTCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAgg                                                                                                                                                     >  3:95340/1‑147 (MQ=255)       aaTTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCCTGCCCCGCATGAAGa                                                                                                                                                 <  1:198873/148‑1 (MQ=255)       aaTTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAgg                                                                                                                                                     >  5:206848/1‑144 (MQ=255)         ttGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATc                                                                                                                                               >  5:364090/1‑148 (MQ=255)           gATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCAt                                                                                                                                            >  4:91884/1‑149 (MQ=255)            aTCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCAtt                                                                                                                                           >  5:298634/1‑149 (MQ=255)             tCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTg                                                                                                                                          <  7:314608/149‑1 (MQ=255)              cAAAAGCTTCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCCACCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTg                                                                                                                                          <  1:76025/149‑1 (MQ=255)               aaaaGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCa                                                                                                                                        >  2:144494/1‑149 (MQ=255)               aaaaGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCa                                                                                                                                        <  4:95340/149‑1 (MQ=255)                      gCCAACACCCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTg                                                                                                                                >  4:323498/1‑149 (MQ=255)                       ccAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGc                                                                                                                               >  6:189593/1‑149 (MQ=255)                       ccAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGc                                                                                                                               >  3:361681/1‑149 (MQ=255)                        cAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCt                                                                                                                              <  3:91884/149‑1 (MQ=255)                        cAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCt                                                                                                                              >  2:323982/1‑149 (MQ=255)                        cAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGGTTCGGCGATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑GCCGCTCGATGATGCTCTTCCTGCGGGGCAGGAAGCTCCCTTGCGGC‑CACTGc                                                                                                                               >  8:254179/1‑148 (MQ=255)                         aaCAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCt                                                                                                                              >  7:237977/1‑148 (MQ=255)                             ggcggcGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCaaaa                                                                                                                         <  1:145211/149‑1 (MQ=255)                             ggcggcGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCaaaa                                                                                                                         <  5:231310/149‑1 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                                     <  3:190124/149‑1 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                                     <  5:189593/149‑1 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                                     <  6:128661/149‑1 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATCATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                                     <  1:101879/149‑1 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGAGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACt                                                                                                                       >  6:20882/1‑147 (MQ=255)                                 gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTCCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                                     <  5:164912/149‑1 (MQ=255)                                  cgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAg                                                                                                                    <  7:184947/149‑1 (MQ=255)                                   gTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAgg                                                                                                                   <  6:364090/149‑1 (MQ=255)                                    tAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAggg                                                                                                                  >  3:380394/1‑149 (MQ=255)                                     aGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATGCCCCCAGGATGCTCGATGCATTCGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAg                                                                                                                    <  1:127625/148‑1 (MQ=255)                                       cTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTAt                                                                                                               <  3:323498/149‑1 (MQ=255)                                          aCGACTTTGCTGCGGCCGGTCGGTAAACACGTTTTCGTCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGGGCG‑CCAGCTCGATGATGATCTTCATGCGGGGCATGAAGATCCATTGCAGC‑CACTGCTCAAAACTTAGGGTATCCa                                                                                                            <  7:346687/149‑1 (MQ=255)                                             ccTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CAGGCTCGATGATGATCTTCATGCGGGGCAGGAACCTCCATTGCAGC‑CACTGCTCATACCTGAGGGGATGCGCgg                                                                                                         >  4:305583/2‑149 (MQ=255)                                                   cTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAAc                                                                                                    <  2:111212/148‑1 (MQ=255)                                                   cTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAAc                                                                                                    <  8:237977/148‑1 (MQ=255)                                                      cggccggTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCGCGGCGAACgg                                                                                                  <  8:26189/147‑1 (MQ=255)                                                      cggccggTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCa                                                                                                >  7:311001/1‑149 (MQ=255)                                                      cggccggTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCa                                                                                                >  3:160478/1‑149 (MQ=255)                                                                gtAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGc                                                                                      <  4:361681/149‑1 (MQ=255)                                                                gaaaaCACGTTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGc                                                                                      <  8:348719/147‑1 (MQ=255)                                                                   aaCACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAg                                                                                   >  5:272750/1‑149 (MQ=255)                                                                    acacGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGt                                                                                  >  1:157669/1‑149 (MQ=255)                                                                      acGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGc                                                                                >  5:102958/1‑149 (MQ=255)                                                                        ggTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCa                                                                                                >  8:326046/1‑131 (MQ=255)                                                                        ggTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCa                                                                                                <  7:326046/131‑1 (MQ=255)                                                                          tttCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCtcg                                                                            >  5:28319/1‑149 (MQ=255)                                                                                ccATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCtt                                                                      <  5:330496/149‑1 (MQ=255)                                                                                ccATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCtt                                                                      <  5:20882/149‑1 (MQ=255)                                                                                ccATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCtt                                                                      <  3:213435/149‑1 (MQ=255)                                                                                ccATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCtt                                                                      >  4:113143/1‑149 (MQ=255)                                                                                  aTG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTgg                                                                    <  4:380394/149‑1 (MQ=255)                                                                                   tG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGACCAGTGCCTCGTCGCTTggc                                                                   >  5:301692/1‑149 (MQ=255)                                                                                       aCAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTggcgg                                                                 >  8:208421/1‑149 (MQ=255)                                                                                        cAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTggcggc                                                                <  8:311001/149‑1 (MQ=255)                                                                                         aaGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCt                                                               >  6:295963/1‑149 (MQ=255)                                                                                           ggATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                             >  1:312583/1‑149 (MQ=255)                                                                                           ggATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                             >  1:307545/1‑149 (MQ=255)                                                                                           ggATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                             >  1:371894/1‑149 (MQ=255)                                                                                           ggATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                             >  5:58223/1‑149 (MQ=255)                                                                                                   gATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCcca                                                     <  7:254179/149‑1 (MQ=255)                                                                                                   gATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCcca                                                     <  2:371894/149‑1 (MQ=255)                                                                                                    aTGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCccac                                                    <  6:28319/149‑1 (MQ=255)                                                                                                    aTGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CAATGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCATTGCCTCGTCGCTTGGCGGCTCGTCATCccac                                                    >  4:309018/1‑149 (MQ=255)                                                                                                            ggCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGACCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                            <  6:301692/149‑1 (MQ=255)                                                                                                              cAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGc                                                                                <  3:324858/111‑1 (MQ=255)                                                                                                              cAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGc                                                                                >  4:324858/1‑111 (MQ=255)                                                                                                              cAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCa                                          <  6:298634/149‑1 (MQ=255)                                                                                                                   ggTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTg                                      >  1:149224/1‑148 (MQ=255)                                                                                                                      gCC‑CCTGCTCGCTGATGTTGTTCATGCGGGGGAGGATGATCCATTTCAGC‑CACTGCTCAAAACTGAGGGTATCAACCGCGAACCGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCAACCCACCAACCTTGCACCTGCAGc                                  <  7:60757/149‑1 (MQ=255)                                                                                                                        c‑ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTc                                <  6:272750/149‑1 (MQ=255)                                                                                                                          ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCtt                                             <  1:281226/135‑1 (MQ=255)                                                                                                                          ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCtt                                             >  2:281226/1‑135 (MQ=255)                                                                                                                          ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                               >  8:42645/1‑149 (MQ=255)                                                                                                                             gCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                            <  2:307545/149‑1 (MQ=255)                                                                                                                             gCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                            <  2:312583/149‑1 (MQ=255)                                                                                                                                 gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAgg                                                                                                                   >  6:201886/1‑58 (MQ=255)                                                                                                                                 gatgatgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAgg                                                                                                                   <  5:201886/58‑1 (MQ=255)                                                                                                                                    atggatCTTCCTGCCCCGCATGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATc                      >  6:118171/4‑148 (MQ=255)                                                                                                                                     atgatCTTCATGCGGGGCAGGAAGATCCATTGCAGCGCACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATccc                                                      >  1:64721/1‑116 (MQ=255)                                                                                                                                     atgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATccc                                                      <  2:64721/115‑1 (MQ=255)                                                                                                                                     atgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATc                      >  5:221361/1‑147 (MQ=255)                                                                                                                                     atgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCaa                    <  1:144494/149‑1 (MQ=255)                                                                                                                                      tgatCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATc                      >  1:29248/1‑146 (MQ=255)                                                                                                                                       gatCTTCCTGCCCCGCATGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCaa                    <  7:38577/147‑1 (MQ=255)                                                                                                                                       gatCTTCCTGCCCCGCATGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCaa                    >  8:38577/1‑147 (MQ=255)                                                                                                                                           ttCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAg              <  7:208421/149‑1 (MQ=255)                                                                                                                                             cATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAg              <  6:206848/147‑1 (MQ=255)                                                                                                                                              aTGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAggtg           <  4:160478/149‑1 (MQ=255)                                                                                                                                                gCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                            <  5:330543/114‑1 (MQ=255)                                                                                                                                                gCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTg                                            >  6:330543/1‑114 (MQ=255)                                                                                                                                                gCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAggtggt         >  3:333139/1‑149 (MQ=255)                                                                                                                                                       aGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAat  <  5:267909/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     GCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTGGTCGGCAAT  >  NC_002947/5319795‑5320090

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