Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,948,133:1 +G 5.2% coding (1565/1746 nt) PP_1749 ← acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,948,1331.G 5.2% 399.9 / 13.3 96coding (1565/1746 nt)PP_1749acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (47/44);  new base G (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 = 9.60e-01

GATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGT  >  NC_002947/1947988‑1948281
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                              <  2:350627/132‑1 (MQ=255)
                      tGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGgca                                                                                                                                              >  1:350627/1‑132 (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                                                                                                                           >  2:837754/1‑134 (MQ=255)
                                       caggccagCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGt                                                                                                                           <  1:837754/134‑1 (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                                                                                                                     >  1:558647/1‑94 (MQ=255)
                                                                                       ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGc                                                                                                                     <  2:558647/94‑1 (MQ=255)
                                                                                       ccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                             >  1:1340817/1‑150 (MQ=255)
                                                                                        cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                            >  2:583788/1‑150 (MQ=255)
                                                                                        cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                            >  2:1919753/1‑150 (MQ=255)
                                                                                        cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACgg                                                            >  2:1919532/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                                         >  1:741722/1‑150 (MQ=255)
                                                                                                           gTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGccggccgg                                         >  2:288830/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)
                                                                                                                                                  gc‑gcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGt  <  2:254215/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAAAACGCTCGGCGGTGGGCTGCGGTTCATGGTTGGCCAGGCCAGCACGTTCGTTGGCCTCGATGATCACGTAC‑TCCGG‑CTGCCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCCGTGCCGCGC‑GCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGT  >  NC_002947/1947988‑1948281

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