Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,213,031 T→A 7.1% L19M (TTG→ATG)  mexC → multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,213,0310TA 7.1% 283.1 / 12.2 98L19M (TTG→ATG) mexCmultidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (52/39);  new base A (2/5);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCGGGC  >  NC_002947/3212884‑3213179
                                                                                                                                                   |                                                                                                                                                    
aTTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTg                                                                                                                                                    >  2:1971807/1‑150 (MQ=255)
aTTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTg                                                                                                                                                    >  1:803430/1‑150 (MQ=255)
  tAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  <  2:2334980/150‑1 (MQ=255)
  tAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  <  1:1691093/150‑1 (MQ=255)
  tAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTACGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  <  1:2235293/150‑1 (MQ=255)
   aaGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGcc                                                                                                                                                 >  2:1848417/1‑150 (MQ=255)
    agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCg                                                                                                                                                >  2:1027977/1‑150 (MQ=255)
     gagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCgg                                                                                                                                               >  2:1810665/1‑150 (MQ=255)
      agTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGCTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCggg                                                                                                                                              >  1:75252/1‑150 (MQ=255)
       gTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                             <  2:806557/150‑1 (MQ=255)
          ttttATGATTTTAATGAGACAATATTGACGCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTGTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGATGGGCGCATTGGCCATGGCGATCGCCTTGACCgggcttg                                                                                                                                          >  2:2270462/1‑146 (MQ=255)
           tttATGATTTTAATGTGTGAATTTTTGCTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTCGATGGTAATGTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGa                                                                                                                                         <  1:734933/150‑1 (MQ=255)
              aTGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                                      >  2:1007465/1‑150 (MQ=255)
                  ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgca                                                                                                                                  >  2:1801570/1‑150 (MQ=255)
                  ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgca                                                                                                                                  <  1:1007465/150‑1 (MQ=255)
                         gcgaCACTATTCAGTCATGTGCGGTTTAAAGCTATCATCCGCCTTGGTTTCAATCAAGGTTCGGGTGTGTGATGGGCAATACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCGGCaccagagc                                                                                                                           >  2:385978/3‑150 (MQ=255)
                          agaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                                          <  1:1810665/150‑1 (MQ=255)
                           gaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCTCagcggagcaa                                                                                                                         >  2:1329876/1‑149 (MQ=255)
                           gaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGCTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagcag                                                                                                                         <  1:1058091/150‑1 (MQ=255)
                            aCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCGGCAGAGCAgg                                                                                                                        >  1:1631132/1‑150 (MQ=255)
                                     aCTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAg                                                                                                               >  2:1175935/1‑150 (MQ=255)
                                      cTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                              <  2:75252/150‑1 (MQ=255)
                                      cTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGc                                                                                                              <  1:1848417/150‑1 (MQ=255)
                                       tCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGcc                                                                                                             >  2:1697384/1‑150 (MQ=255)
                                         aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                                      >  2:2167726/1‑123 (MQ=255)
                                         aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                                      <  1:2167726/123‑1 (MQ=255)
                                          ggtgGGTTTGCAAAGCAGCATCCGCCCAGTTTTCAATCAAGGTTCGGTTGTGTTATGGGCAAGACAGGTGCGATAAGGGTGGGCGCATTGGCCATGGCTATCGCCTTGGCCGGGTGTGAGCCCGGAGCAGAGCAGGACGGTGCAGCCAAt                                                                                                          <  2:1631132/149‑1 (MQ=255)
                                              ggTTTGAAAAGTAGCATCCGCCCCGTCTTCAATCAAGATTCGGGCGTGTGATGGGCAAGACAGGTTCAATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                                      <  1:1329876/150‑1 (MQ=255)
                                               gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAgg                                                                                                                        <  1:1107972/131‑1 (MQ=255)
                                               gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAgg                                                                                                                        >  2:1107972/1‑131 (MQ=255)
                                               gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAg                                                                                                     <  1:352700/150‑1 (MQ=255)
                                                tttGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGc                                                                                                                                                  >  1:1941234/1‑104 (MQ=255)
                                                tttGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGc                                                                                                                                                  <  2:1941234/104‑1 (MQ=255)
                                                      aaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  >  2:257620/1‑98 (MQ=255)
                                                      aaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  <  1:257620/98‑1 (MQ=255)
                                                       aGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCgg                                                                                             <  2:107999/150‑1 (MQ=255)
                                                         tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                           <  1:1436096/150‑1 (MQ=255)
                                                               ccGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGGCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                                     >  1:850614/1‑150 (MQ=255)
                                                               ccGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  >  1:690492/1‑89 (MQ=255)
                                                               ccGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                                  <  2:690492/89‑1 (MQ=255)
                                                                 gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                                   <  1:1697384/150‑1 (MQ=255)
                                                                   cccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGa                                                                                 <  1:976777/150‑1 (MQ=255)
                                                                     cGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACa                                                                               >  2:1208579/1‑150 (MQ=255)
                                                                      gTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACAt                                                                              >  1:1621922/1‑150 (MQ=255)
                                                                        tttCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTg                                                                            <  1:1077704/150‑1 (MQ=255)
                                                                             atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGcc                                                                       >  2:493980/1‑150 (MQ=255)
                                                                                 aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                             <  2:320475/76‑1 (MQ=255)
                                                                                 aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGgt                                                                                                                                             >  1:320475/1‑76 (MQ=255)
                                                                                   gTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCa                                                                 >  1:576961/1‑150 (MQ=255)
                                                                                    ttCGGGTGTGTGATGGGCCAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGCACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAg                                                                <  1:1287906/150‑1 (MQ=255)
                                                                                        ggTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCgg                                                                                                                                               <  1:1743402/67‑1 (MQ=255)
                                                                                        ggTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCgg                                                                                                                                               >  2:1743402/1‑67 (MQ=255)
                                                                                          tgtgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGcc                                                          >  1:113575/1‑150 (MQ=255)
                                                                                            tgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGTTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTg                                                                                  <  2:1892021/124‑1 (MQ=255)
                                                                                            tgtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCAATGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGTTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTg                                                                                  >  1:1892021/1‑124 (MQ=255)
                                                                                                aTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGATGTAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCGCGAGCAGGTCGCCTTggcg                                                    >  2:142124/1‑150 (MQ=255)
                                                                                                  gggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggcggc                                                  >  1:642671/1‑150 (MQ=255)
                                                                                                   ggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                                 >  2:1979606/1‑150 (MQ=255)
                                                                                                   ggCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                                 >  2:2316518/1‑150 (MQ=255)
                                                                                                    gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                     <  1:120481/145‑1 (MQ=255)
                                                                                                    gCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTggc                                                     >  2:120481/1‑145 (MQ=255)
                                                                                                     cAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTTCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGa                                               >  2:1066895/1‑150 (MQ=255)
                                                                                                     cAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGa                                               >  1:334293/1‑150 (MQ=255)
                                                                                                      aaGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGCGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGACATCGTGACCTTCGCGCCCGAGCAGGGCTCCATGGCGGCTGaa                                              >  2:1301995/1‑150 (MQ=255)
                                                                                                      aaGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGGGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTTaa                                              <  2:803430/150‑1 (MQ=255)
                                                                                                          cAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgc                                          >  2:1707637/1‑150 (MQ=255)
                                                                                                          cAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgc                                          <  1:493980/150‑1 (MQ=255)
                                                                                                           aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgcc                                         >  2:736315/1‑150 (MQ=255)
                                                                                                           aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgcc                                         <  2:216640/150‑1 (MQ=255)
                                                                                                             gTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgg                                       >  1:1547132/1‑150 (MQ=255)
                                                                                                              ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccggg                                      >  2:1111059/1‑150 (MQ=255)
                                                                                                              ttCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGGGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccggg                                      >  1:1277647/1‑150 (MQ=255)
                                                                                                               tCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggc                                     <  1:1208579/150‑1 (MQ=255)
                                                                                                                  aTAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggccgg                                  >  2:1581034/1‑150 (MQ=255)
                                                                                                                   tAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                                     <  2:2243778/98‑1 (MQ=255)
                                                                                                                   tAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATc                                                                                     >  1:2243778/1‑98 (MQ=255)
                                                                                                                   tAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGccgggccggg                                 >  1:800154/1‑150 (MQ=255)
                                                                                                                      gggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACt                                            <  2:179512/136‑1 (MQ=255)
                                                                                                                      gggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACt                                            >  1:179512/1‑136 (MQ=255)
                                                                                                                       ggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGa                             <  1:1581034/150‑1 (MQ=255)
                                                                                                                        gTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAg                            >  2:1456775/1‑150 (MQ=255)
                                                                                                                        gTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAg                            >  2:1050231/1‑150 (MQ=255)
                                                                                                                           ggCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACgg                                                                                                                    <  1:662923/59‑1 (MQ=255)
                                                                                                                           ggCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACgg                                                                                                                    >  2:662923/1‑59 (MQ=255)
                                                                                                                           ggCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGTGCCCGAGCAGGTc                                                             >  2:2179458/1‑114 (MQ=255)
                                                                                                                           ggCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGTGCCCGAGCAGGTc                                                             <  1:2179458/114‑1 (MQ=255)
                                                                                                                               cATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc                     <  1:736315/150‑1 (MQ=255)
                                                                                                                               cATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc                     <  2:2209495/150‑1 (MQ=255)
                                                                                                                                 ttGGACCTGGCGATCGCGTTGGCCGGGTTTGCTCCGGTAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgcgc                   <  2:1038229/150‑1 (MQ=255)
                                                                                                                                    gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                           <  1:1747386/75‑1 (MQ=255)
                                                                                                                                    gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTg                                                                                           >  2:1747386/1‑75 (MQ=255)
                                                                                                                                    gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTc                >  2:1762230/1‑150 (MQ=255)
                                                                                                                                     ccATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCg                                                                                                       <  1:1137217/62‑1 (MQ=255)
                                                                                                                                     ccATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCg                                                                                                       >  2:1137217/1‑62 (MQ=255)
                                                                                                                                         ggCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGa           <  1:1111059/150‑1 (MQ=255)
                                                                                                                                                gCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTTCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCgg    <  1:1066895/150‑1 (MQ=255)
                                                                                                                                                 ccTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGCCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCACCGAGGTACggg   >  1:395784/1‑150 (MQ=255)
                                                                                                                                                  cTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCGGgc  >  1:1954478/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGCCGAGGTGCGGGC  >  NC_002947/3212884‑3213179

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