Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,319,933	A→C	5.2%	M62R (ATG→AGG)	PP_4681 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,319,933	0	A	C	5.2%	285.8 / 3.9	89	M62R (ATG→AGG)	PP_4681	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base C (3/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTG  >  NC_002947/5319786‑5320081                                                                                                                                                       |                                                                                                                                                      gCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCAt                                                                                                                                                      >  4:213435/1‑149 (MQ=255) gCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCAt                                                                                                                                                      >  6:330496/1‑149 (MQ=255)          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                                                                                                            <  5:189593/149‑1 (MQ=255)                                          gcgTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGa                                                                                                            <  3:190124/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                                                                                           <  8:237977/148‑1 (MQ=255)                                                            cTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAAc                                                                                           <  2:111212/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                                                                                       <  7:326046/131‑1 (MQ=255)                                                                                 ggTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCa                                                                                       >  8:326046/1‑131 (MQ=255)                                                                                   tttCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCtcg                                                                   >  5:28319/1‑149 (MQ=255)                                                                                         ccATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCtt                                                             >  4:113143/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)                                                                                           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:371894/1‑149 (MQ=255)                                                                                                    ggATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCg                                                    >  1:307545/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                                                                       >  4:324858/1‑111 (MQ=255)                                                                                                                       cAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGc                                                                       <  3:324858/111‑1 (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                                    >  2:281226/1‑135 (MQ=255)                                                                                                                                   ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCtt                                    <  1:281226/135‑1 (MQ=255)                                                                                                                                   ccAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTcg                      >  8:42645/1‑149 (MQ=255)                                                                                                                                      gCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                   <  2:312583/149‑1 (MQ=255)                                                                                                                                      gCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCgcgc                   <  2:307545/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)                                                                                                                                                       |                                                                                                                                                      GCGGAAGGAGCGATCAATTGATCAAAAGCTGCCAACAGGCGGCGTAGCTCACGACTTTGCTCCGGCCGGTCGGTAAACACGGTTTCGGCCATG‑‑ACAAGGATGCTCGATGCATTGGGCAGCGGGTGCC‑CCAGCTCGATGATGATCTTCATGCGGGGCAGGAAGATCCATTGCAGC‑CACTGCTCAAAACTGAGGGTATCCACGGCGAACGGCACCGTACTGGCCAGTGCCTCGTCGCTTGGCGGCTCGTCATCCCACCAACCTTGCACCTGCAGCTCGCGCTCGATCAAGAGCAGGTG  >  NC_002947/5319786‑5320081

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