Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,411 G→C 15.6% G94A (GGC→GCC)  PP_1181 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,4110GC15.6% 235.0 / 37.4 96G94A (GGC→GCC) PP_1181two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (48/33);  new base C (7/8);  total (55/41)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGC  >  NC_002947/1356263‑1356559
                                                                                                                                                    |                                                                                                                                                      
agcatcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGc                                                                                                                                                       <  1:264151/148‑1 (MQ=255)
    tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTgg                                                                                                                                                   >  1:1464187/1‑150 (MQ=255)
     cGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGATCTCAGTGATGAACTGGACCAGGTCCAGGGCCTGGa                                                                                                                                                  <  1:802796/150‑1 (MQ=255)
      gACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGAc                                                                                                                                                 <  2:2396709/150‑1 (MQ=255)
         agccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTg                                                                                                                                              <  1:504627/150‑1 (MQ=255)
         agccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTg                                                                                                                                              <  2:1386171/150‑1 (MQ=255)
          gccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTgg                                                                                                                                             <  1:1062811/150‑1 (MQ=255)
            cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGCCCAGGTCCAGGCCCTGGACCTggtc                                                                                                                                           <  2:1456460/150‑3 (MQ=255)
            cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCAGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCATGTCCAGGGCCTGGACCTGGgc                                                                                                                                           <  1:680718/150‑1 (MQ=255)
             agccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgcg                                                                                                                                          >  2:1135239/1‑150 (MQ=255)
              gccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgcgc                                                                                                                                         <  1:2273953/150‑1 (MQ=255)
                 cGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGCGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcga                                                                                                                                      <  2:509038/150‑1 (MQ=255)
                           aTCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTCCGTc                                                                                                                            >  1:843879/1‑150 (MQ=255)
                              cTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGc                                                                                                                         >  2:70221/1‑150 (MQ=255)
                              cTCCACCTGATGCGCCCCGGCTAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGc                                                                                                                         >  2:704312/1‑150 (MQ=255)
                               tCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCCGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCa                                                                                                                        >  1:822277/1‑150 (MQ=255)
                                     tGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCg                                                                                                                  >  2:1961313/1‑150 (MQ=255)
                                      gATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGt                                                                                                                 >  1:192368/1‑150 (MQ=255)
                                      gATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGt                                                                                                                 >  1:1697491/1‑150 (MQ=255)
                                      gATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGt                                                                                                                 <  2:1687248/150‑1 (MQ=255)
                                        tGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgc                                                                                                               >  2:1332150/1‑150 (MQ=255)
                                         gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgcg                                                                                                              >  1:1392313/1‑150 (MQ=255)
                                                        cGGCCTGGGCATCTTTCGCCTTGTGCGTCGCCAGTACCAGGGGCCGATACTCATGCTCACCGCGCGCAGTGATGAACCGGACCAGCTCCAGGGCCTGGACCTGGGCGCCCACGACTACGTCTGCAAACGCGTGCGCCCACG‑‑CTTgctgct                                                                                             <  2:843111/150‑1 (MQ=255)
                                                          gCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTgg                                                                                           <  1:1591731/150‑1 (MQ=255)
                                                            cTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTTCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCa                                                                                         >  2:1599292/1‑150 (MQ=255)
                                                             tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCAc                                                                                        <  1:1135239/150‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATcc                                                                                  >  1:1401754/1‑150 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGAACACAGCCAGGGCCTGGACCTGGGCGCAGCCGATTACGGCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATcc                                                                                  >  1:1740270/1‑150 (MQ=255)
                                                                          ccaggTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCct                                                                           >  1:1290796/4‑150 (MQ=255)
                                                                            gggTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCctgc                                                                         >  1:19615/1‑150 (MQ=255)
                                                                                  gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcga                                                                                                                                      >  2:1703262/1‑85 (MQ=255)
                                                                                  gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcga                                                                                                                                      <  1:1703262/85‑1 (MQ=255)
                                                                                      ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                               >  1:1780507/1‑150 (MQ=255)
                                                                                      ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCAAAG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                               >  2:1397534/1‑150 (MQ=255)
                                                                                      ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCGGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAg                                                               >  2:1554619/1‑150 (MQ=255)
                                                                                            ccAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAc                                                         >  2:1848089/1‑150 (MQ=255)
                                                                                            ccAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAc                                                         >  2:453128/1‑150 (MQ=255)
                                                                                              aGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCg                                                       >  1:1005502/1‑150 (MQ=255)
                                                                                               ggggCCGATCCTCATGCTCACCGCGCGCAGTGATGCACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                                      <  2:1290796/150‑1 (MQ=255)
                                                                                                ggcccGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTg                                                     >  1:2402158/4‑150 (MQ=255)
                                                                                                  gCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGa                                                   >  1:215837/1‑150 (MQ=255)
                                                                                                     gATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGaaa                                                          <  1:310452/140‑1 (MQ=255)
                                                                                                     gATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGaaa                                                          >  2:310452/1‑140 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgc                                           >  1:1381914/1‑150 (MQ=255)
                                                                                                            aTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                                      >  2:1148493/1‑137 (MQ=255)
                                                                                                            aTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                                      <  1:1148493/137‑1 (MQ=255)
                                                                                                            aTGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcgc                                         <  2:1780507/150‑1 (MQ=255)
                                                                                                              gCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCACACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCTCCGCAGCGAAACCGTGGACAGCAAGCGCCa                                       <  1:704312/150‑1 (MQ=255)
                                                                                                              catCACTGCGCGCGGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGcc                                                                             >  2:1381695/3‑112 (MQ=255)
                                                                                                              catCACTGCGCGCGGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGcc                                                                             <  1:1381695/110‑1 (MQ=255)
                                                                                                                tCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAgg                                     >  1:561889/1‑150 (MQ=255)
                                                                                                                    cgcgcgCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                    <  1:2256725/131‑1 (MQ=255)
                                                                                                                    cgcgcgCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                    >  2:2256725/1‑131 (MQ=255)
                                                                                                                    cgcgcgCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCt                                 <  1:1961313/150‑1 (MQ=255)
                                                                                                                     gcgcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCTCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTg                                <  1:1397534/150‑1 (MQ=255)
                                                                                                                     gcgcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTg                                <  1:70221/150‑1 (MQ=255)
                                                                                                                     gcgcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTg                                >  1:1286274/1‑150 (MQ=255)
                                                                                                                      cgcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                               <  1:453128/150‑1 (MQ=255)
                                                                                                                          cAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTAGACAGCAAGCGCCAGGTCCTGGcgct                           <  2:822277/150‑4 (MQ=255)
                                                                                                                              gatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGt                       >  2:226618/1‑150 (MQ=255)
                                                                                                                               atgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCtgctgc                                                                      >  2:751505/1‑102 (MQ=255)
                                                                                                                               atgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCtgctgc                                                                      <  1:751505/102‑1 (MQ=255)
                                                                                                                               atgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTg                      <  2:1401754/150‑1 (MQ=255)
                                                                                                                               atgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTg                      >  2:311293/1‑150 (MQ=255)
                                                                                                                                  aaCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCt                                                                                                                           <  1:39761/48‑1 (MQ=255)
                                                                                                                                  aaCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCt                                                                                                                           >  2:39761/1‑48 (MQ=255)
                                                                                                                                  aaCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTgcgc                   >  2:153659/1‑150 (MQ=255)
                                                                                                                                   aCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCt                  >  1:1288186/1‑150 (MQ=255)
                                                                                                                                   aCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGCGCGCCAACG‑‑CTTACCGGTGGATCATGTCCATGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCt                  >  1:1496856/1‑150 (MQ=255)
                                                                                                                                    cTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTg                 >  1:93265/1‑150 (MQ=255)
                                                                                                                                    cTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCAACGAACCGGCTTCTTCGCCTTCGTGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTg                 >  2:1891777/1‑150 (MQ=255)
                                                                                                                                     tGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAATCCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCTCTGc                <  2:1697491/150‑1 (MQ=255)
                                                                                                                                     tGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc                >  2:1117590/1‑150 (MQ=255)
                                                                                                                                      ggACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGcc                                                                                                            <  1:1354959/59‑1 (MQ=255)
                                                                                                                                      ggACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGcc                                                                                                            >  2:1354959/1‑59 (MQ=255)
                                                                                                                                         ccAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTga                <  1:1892912/146‑2 (MQ=255)
                                                                                                                                         ccAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc                >  2:1892912/1‑146 (MQ=255)
                                                                                                                                         ccAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                    <  2:78325/110‑1 (MQ=255)
                                                                                                                                         ccAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTgg                                                    >  1:78325/1‑110 (MQ=255)
                                                                                                                                         ccAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt            >  2:2249477/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTAt            >  2:720178/1‑150 (MQ=255)
                                                                                                                                          cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcgc                                         >  1:480342/1‑120 (MQ=255)
                                                                                                                                          cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcgc                                         <  2:480342/120‑1 (MQ=255)
                                                                                                                                          cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc           <  2:2402158/150‑1 (MQ=255)
                                                                                                                                          cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCgg                        >  2:1842279/1‑137 (MQ=255)
                                                                                                                                          cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCgg                        <  1:1842279/137‑1 (MQ=255)
                                                                                                                                          cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCCGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc           >  2:1438019/1‑150 (MQ=255)
                                                                                                                                           aGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                               <  1:813475/129‑1 (MQ=255)
                                                                                                                                           aGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTgg                               >  2:813475/1‑129 (MQ=255)
                                                                                                                                           aGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                         <  2:1858178/135‑1 (MQ=255)
                                                                                                                                           aGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                         >  1:1858178/1‑135 (MQ=255)
                                                                                                                                            ggTCCAGGCCCTGGACCTGGTCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                                    <  1:1774497/91‑1 (MQ=255)
                                                                                                                                            ggTCCAGGCCCTGGACCTGGTCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                                    >  2:1774497/1‑91 (MQ=255)
                                                                                                                                                cAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAc     >  2:290506/1‑150 (MQ=255)
                                                                                                                                                 aGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAcc    <  2:1464187/150‑1 (MQ=255)
                                                                                                                                                 aGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAAcc    >  1:2251497/1‑150 (MQ=255)
                                                                                                                                                  gggCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAg                                                                                <  1:283260/73‑1 (MQ=255)
                                                                                                                                                  gggCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAg                                                                                >  2:283260/1‑73 (MQ=255)
                                                                                                                                                   ggCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGc  >  1:1304023/1‑150 (MQ=255)
                                                                                                                                                   ggCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAAGACCTGGCCTTCGATGCGCAGCGTATCGACAACCGc  >  2:1516472/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACG‑‑CTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGACAACCGC  >  NC_002947/1356263‑1356559

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