Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,948,128 Δ1 bp 5.2% coding (1570/1746 nt) PP_1749 ← acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,948,1280C. 5.2% 306.9 / 10.2 96coding (1570/1746 nt)PP_1749acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (47/44);  new base . (3/2);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

AGCAGGTCGATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGC  >  NC_002947/1947980‑1948274
                                                                                                                                                      |                                                                                                                                                   
agcagGTCGATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGC‑CAGAACCGCGTCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGcc                                                                                                                                                    >  2:1698543/1‑150 (MQ=255)
        gATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCAc                                                                                                                                           >  2:1874067/1‑150 (MQ=255)
        gATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCAc                                                                                                                                           >  2:1874750/1‑150 (MQ=255)
        gATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCAc                                                                                                                                           >  2:1874774/1‑150 (MQ=255)
            aaaCGCTTGGCGGTGGCATGCGGTTTGTTGGTGGCCAGGCCAGCACGTCAGTCGGCTTC‑‑TCTTCGGGTGCTTTCGGTCTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgca                                                                                                                                       <  1:975396/150‑1 (MQ=255)
             aaCGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgcag                                                                                                                                      >  1:1403413/1‑150 (MQ=255)
             aaCGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCACGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTG‑CGCGCGGCACgggca                                                                                                                                      <  1:535613/150‑4 (MQ=255)
                gCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgcagcac                                                                                                                                   >  1:72481/1‑149 (MQ=255)
                gCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCAt                                                                                                                                   <  1:1874774/150‑1 (MQ=255)
                gCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCAt                                                                                                                                   <  1:1874750/150‑1 (MQ=255)
                gCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCAt                                                                                                                                   <  1:1874067/150‑1 (MQ=255)
                  tCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGGATCg                                                                                                                                 >  1:75069/1‑150 (MQ=255)
                    ggcggTGGGCTTCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGc                                                                                                                               <  1:940728/150‑1 (MQ=255)
                    ggcggTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGc                                                                                                                               >  2:338252/1‑150 (MQ=255)
                         tGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCATCa                                                                                                                          <  2:72481/150‑1 (MQ=255)
                             cTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCgg                                                                                                                      <  1:1764722/150‑1 (MQ=255)
                              tGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgca                                                                                                                                       >  1:350627/1‑132 (MQ=255)
                              tGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgca                                                                                                                                       <  2:350627/132‑1 (MQ=255)
                                 ggTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGc                                                                                                                  <  2:487934/150‑1 (MQ=255)
                                   ttCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCTTGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                <  1:151269/150‑1 (MQ=255)
                                   ttCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                <  2:121113/150‑1 (MQ=255)
                                   ttCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAg                                                                                                                >  2:491442/1‑150 (MQ=255)
                                    tCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAgg                                                                                                               <  2:1042341/150‑1 (MQ=255)
                                       tGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGt                                                                                                            >  2:1878645/1‑150 (MQ=255)
                                        ggTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGtt                                                                                                           <  2:1181150/150‑1 (MQ=255)
                                            ggccaggccaGCACGTTCGTTGGCCTCGATGATCACGAAA‑CAAAA‑TAAAACTGTATTAGTCACCATTAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGt                                                                                                       <  1:1869529/150‑1 (MQ=255)
                                            ggccaggcatGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGt                                                                                                       <  1:1021809/150‑1 (MQ=255)
                                               caggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGt                                                                                                                    <  1:837705/134‑1 (MQ=255)
                                               caggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGt                                                                                                                    <  1:837754/134‑1 (MQ=255)
                                               caggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGt                                                                                                                    >  2:837754/1‑134 (MQ=255)
                                               caggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGt                                                                                                                    >  2:837705/1‑134 (MQ=255)
                                                aggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACg                                                                                                   >  2:454686/1‑150 (MQ=255)
                                                        aCGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAg                                                                                           >  1:1120994/1‑150 (MQ=255)
                                                            tcgttGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTg                                                                                       >  1:803886/1‑150 (MQ=255)
                                                             cgttGGCCTCGATGATCACGTTC‑TCCGG‑CTGCCCGCCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCTTGCCGCGC‑GCACGGCAGCCTCGGCCAGCTCCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTgc                                                                                      <  1:1714549/150‑1 (MQ=255)
                                                              gttGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGtcct                                                                                               <  2:257161/140‑1 (MQ=255)
                                                              gttGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGtcct                                                                                               >  1:257161/1‑140 (MQ=255)
                                                                tGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTgccgc                                                                                   <  1:491442/150‑1 (MQ=255)
                                                                  gCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGccgccg                                                                                 >  1:745414/1‑150 (MQ=255)
                                                                   ccTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCgg                                                                                >  1:702936/1‑150 (MQ=255)
                                                                    cTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCGCGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGgt                                                                               <  2:1880747/150‑1 (MQ=255)
                                                                    cTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGgt                                                                               >  1:1218319/1‑150 (MQ=255)
                                                                    cTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGgt                                                                               >  1:844809/1‑150 (MQ=255)
                                                                           atCACGTAC‑TCCGG‑CTGCCCGGCATCGTGAACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAgg                                                                        >  1:528552/1‑150 (MQ=255)
                                                                            tCACGTAC‑TCCGG‑CTACCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGa                                                                                                     >  1:932530/1‑120 (MQ=255)
                                                                            tCACGTAC‑TCCGG‑CTACCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGa                                                                                                     <  2:932530/120‑1 (MQ=255)
                                                                              aCGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTg                                                                     >  1:265799/1‑150 (MQ=255)
                                                                              aCGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTg                                                                     >  2:112131/1‑150 (MQ=255)
                                                                              aCGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTg                                                                     >  2:1106824/1‑150 (MQ=255)
                                                                              aCGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTGCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTg                                                                     <  2:840117/150‑1 (MQ=255)
                                                                               cGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTgg                                                                    <  2:745414/150‑1 (MQ=255)
                                                                                gTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGc                                                                   <  1:1854444/150‑1 (MQ=255)
                                                                                  aC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCg                                                                 >  2:788474/1‑150 (MQ=255)
                                                                                     tCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGGCCTCCAGGGGGCCGCCGGTGTGCAGGTTGGCCCTg                                                               >  2:435176/1‑150 (MQ=255)
                                                                                      ccGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgc                                                              >  1:1630877/1‑150 (MQ=255)
                                                                                        gg‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgcgc                                                            >  2:127319/1‑150 (MQ=255)
                                                                                        gg‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTgcgc                                                            >  2:127325/1‑150 (MQ=255)
                                                                                         g‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgcc                                                           >  1:582952/1‑150 (MQ=255)
                                                                                            tGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgccgc                                                         >  2:1826582/1‑150 (MQ=255)
                                                                                            tGCCCGGCATCACGCACCATGAAGTATAACCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGcgccgc                                                         >  1:1297421/1‑150 (MQ=255)
                                                                                               ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGc                                                                                                              <  2:558647/94‑1 (MQ=255)
                                                                                               ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGc                                                                                                              >  1:558647/1‑94 (MQ=255)
                                                                                               ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                      >  1:1340817/1‑150 (MQ=255)
                                                                                                cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                     >  2:583788/1‑150 (MQ=255)
                                                                                                cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                     >  2:1919532/1‑150 (MQ=255)
                                                                                                cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                     >  2:1919753/1‑150 (MQ=255)
                                                                                                  gCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGcc                                                   <  1:1106824/150‑1 (MQ=255)
                                                                                                      cACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGgc                                               <  2:383684/150‑1 (MQ=255)
                                                                                                         gCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCt                                            >  2:516283/1‑150 (MQ=255)
                                                                                                         gCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGATGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCt                                            >  1:570728/1‑150 (MQ=255)
                                                                                                           aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                       <  1:1023427/137‑1 (MQ=255)
                                                                                                           aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCAc                                                       >  2:1023427/1‑137 (MQ=255)
                                                                                                           aCCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                          <  2:265799/150‑1 (MQ=255)
                                                                                                                   gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGccggccgg                                  >  2:288830/1‑150 (MQ=255)
                                                                                                                   gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGccggccgg                                  >  1:741722/1‑150 (MQ=255)
                                                                                                                       aGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAac                              <  2:803886/150‑1 (MQ=255)
                                                                                                                        gCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAaca                             <  1:112131/150‑1 (MQ=255)
                                                                                                                         ccccACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAacac                            >  1:609254/1‑150 (MQ=255)
                                                                                                                         ccccACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAacac                            <  2:1072458/150‑1 (MQ=255)
                                                                                                                          cccACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTTGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcg                           >  2:1450949/1‑150 (MQ=255)
                                                                                                                            caccacCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTCCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtc                         <  2:1488794/150‑1 (MQ=255)
                                                                                                                               caccGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCAACCCAGGGGCCGGCCGGCAACACgtcgtc                      <  1:1698543/150‑1 (MQ=255)
                                                                                                                                accGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCa                     <  1:788474/150‑1 (MQ=255)
                                                                                                                                  cGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATa                   >  1:1036014/1‑150 (MQ=255)
                                                                                                                                   ggAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAg                  <  1:454686/150‑1 (MQ=255)
                                                                                                                                   ggAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAg                  >  2:1085250/1‑150 (MQ=255)
                                                                                                                                    gAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCTGTTTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGc                 <  1:127319/150‑1 (MQ=255)
                                                                                                                                    gAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGc                 <  1:127325/150‑1 (MQ=255)
                                                                                                                                          gatgctgccGTG‑CGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGc           >  1:1620349/8‑150 (MQ=255)
                                                                                                                                           cAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGcc          <  1:1878645/150‑1 (MQ=255)
                                                                                                                                            aGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCTGGTATTCAATGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGccc         <  2:528552/150‑1 (MQ=255)
                                                                                                                                               tgccGTG‑CGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGcc                                                                  >  1:349767/3‑90 (MQ=255)
                                                                                                                                               tgccGTG‑CGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGcc                                                                  <  2:349767/88‑1 (MQ=255)
                                                                                                                                               tgccGTG‑CGCGCGGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCt      >  1:1159384/3‑150 (MQ=255)
                                                                                                                                                 ccGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgc    <  1:1826582/150‑1 (MQ=255)
                                                                                                                                                   gTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCGACACGTCGTCATAGCCAACGCCCGCTgcgc  <  2:1120994/150‑1 (MQ=255)
                                                                                                                                                   gTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgcgc  <  1:583788/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGCAGGTCGATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGC  >  NC_002947/1947980‑1948274

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