Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,879,983 A→G 6.2% S592G (AGC→GGC)  mexF → multidrug RND transporter MexF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,879,9830AG 6.2% 190.5 / 2.3 64S592G (AGC→GGC) mexFmultidrug RND transporter MexF
Reads supporting (aligned to +/- strand):  ref base A (22/38);  new base G (2/2);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

GGGTCATCCGCTCCAGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCAACAGCGGC  >  NC_002947/3879834‑3880111
                                                                                                                                                     |                                                                                                                                
gggTCATCCGCTCCAGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCa                                                                                                                                  <  1:1454541/150‑1 (MQ=255)
   tCATCCGCTCCAGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGcc                                                                                                                               >  2:81677/1‑150 (MQ=255)
            ccAGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTa                                                                                                                      <  2:972864/150‑1 (MQ=255)
                  gCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAg                                                                                                                <  2:206239/150‑1 (MQ=255)
                    aTCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCg                                                                                                              <  2:902989/150‑1 (MQ=255)
                        cccTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGa                                                                                                          <  1:449828/150‑1 (MQ=255)
                        cccTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGa                                                                                                          <  1:537113/150‑1 (MQ=255)
                         ccTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTa                                                                                                                      <  1:714385/137‑1 (MQ=255)
                         ccTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTa                                                                                                                      >  2:714385/1‑137 (MQ=255)
                            gTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             >  1:609745/1‑143 (MQ=255)
                            gTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             <  2:609745/143‑1 (MQ=255)
                                  gtATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCAt                                                                                                <  2:751752/150‑1 (MQ=255)
                                  gtATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCAt                                                                                                <  1:1591403/150‑1 (MQ=255)
                                   tATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTAGACCGTACCGAAGCGGTGATCAAGCGCATg                                                                                               >  2:1070128/1‑150 (MQ=255)
                                    aTGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGa                                                                                              <  1:667193/150‑1 (MQ=255)
                                        cGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCGGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGa                                                                                          >  1:162407/1‑150 (MQ=255)
                                        cGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCGGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGa                                                                                          >  1:162686/1‑150 (MQ=255)
                                        cGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             <  1:792686/131‑1 (MQ=255)
                                        cGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             >  2:792686/1‑131 (MQ=255)
                                         ggCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACACCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGaa                                                                                         >  2:944947/1‑150 (MQ=255)
                                          gCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGt                                                                                                            >  2:1000274/1‑130 (MQ=255)
                                          gCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGt                                                                                                            <  1:1000274/130‑1 (MQ=255)
                                             tGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCGGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCg                                                                                     <  2:162686/150‑1 (MQ=255)
                                             tGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCGGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCg                                                                                     <  2:162407/150‑1 (MQ=255)
                                               aTGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGcc                                                                                   <  1:685469/150‑1 (MQ=255)
                                               aTGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCATCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGcc                                                                                   <  2:1197803/150‑1 (MQ=255)
                                                    tCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGaa                                                                              >  2:1266926/1‑150 (MQ=255)
                                                         cctacctCGGCCTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAagcagc                                                                         <  2:1082715/150‑1 (MQ=255)
                                                          ctccctCGGCTCCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCTTCCGCCCAACTGCCGGACGCGGCCAGCCCCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGcc                                                                        <  1:734412/147‑1 (MQ=255)
                                                            acctCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCTGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTg                                                                      <  1:1570404/150‑1 (MQ=255)
                                                            acctCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCt                                                                       >  1:1015692/1‑149 (MQ=255)
                                                            acctCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCt                                                                       <  2:1015692/149‑1 (MQ=255)
                                                                 ggCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTg                                                                 <  2:444011/150‑1 (MQ=255)
                                                                 ggCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTg                                                                 <  1:1178324/150‑1 (MQ=255)
                                                                   cTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGc                                                               <  1:81677/150‑1 (MQ=255)
                                                                             aCGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTa                                                     <  2:290264/150‑1 (MQ=255)
                                                                               gCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGc                                                   <  1:1267814/150‑1 (MQ=255)
                                                                                  gACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCtt                                                >  2:43886/1‑150 (MQ=255)
                                                                                      ggTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCg                                            <  2:84185/150‑1 (MQ=255)
                                                                                               ccGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCg                                   >  2:940841/1‑150 (MQ=255)
                                                                                                cGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAag                                                                             <  1:570370/107‑1 (MQ=255)
                                                                                                cGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAag                                                                             >  2:570370/1‑107 (MQ=255)
                                                                                                    ccAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATc                                >  2:161512/1‑148 (MQ=255)
                                                                                                    ccAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATc                                <  1:161512/148‑1 (MQ=255)
                                                                                                     cAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAAc                             <  2:1239222/150‑1 (MQ=255)
                                                                                                     cAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAAc                             >  2:270917/1‑150 (MQ=255)
                                                                                                     cAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGCACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAAc                             <  1:14564/150‑1 (MQ=255)
                                                                                                        gACAAGCAGTACCTGGTGGCCTTCGCCTAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGt                          >  1:1215882/1‑150 (MQ=255)
                                                                                                             gCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             >  2:1525551/1‑62 (MQ=255)
                                                                                                             gCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCgg                                                                                                             <  1:1525551/62‑1 (MQ=255)
                                                                                                             gCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCAc                     >  2:1153831/1‑150 (MQ=255)
                                                                                                              cAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCAcc                    <  1:1602956/150‑1 (MQ=255)
                                                                                                                 tACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGt                                     >  2:1523467/1‑130 (MQ=255)
                                                                                                                 tACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGt                                     <  1:1523467/130‑1 (MQ=255)
                                                                                                                    cTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGATTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGc              <  1:1065463/150‑1 (MQ=255)
                                                                                                                    cTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGc              <  1:648034/150‑1 (MQ=255)
                                                                                                                      ggtggCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGccc            <  1:43886/150‑1 (MQ=255)
                                                                                                                       gtggCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTAc                                                                                                                     >  2:432398/1‑44 (MQ=255)
                                                                                                                       gtggCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTAc                                                                                                                     <  1:432398/44‑1 (MQ=255)
                                                                                                                         ggCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGTCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCaa         <  1:1266926/150‑1 (MQ=255)
                                                                                                                         ggCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCac         <  1:940841/150‑2 (MQ=255)
                                                                                                                            cTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACa                <  2:1011869/140‑1 (MQ=255)
                                                                                                                            cTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACa                >  1:1011869/1‑140 (MQ=255)
                                                                                                                            cTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCAACAg      >  1:1127409/1‑150 (MQ=255)
                                                                                                                                gCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCAACAGCGGc  >  1:183304/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                
GGGTCATCCGCTCCAGCGGCATCGCCCTGTTCGTGTATGCCGGCCTGATGGGTCTGACCTACCTCGGCTTCTCGTCCACGCCGACCGGTTTCGTGCCGGCCCAGGACAAGCAGTACCTGGTGGCCTTCGCCCAACTGCCGGACGCGGCCAGCCTCGACCGTACCGAAGCGGTGATCAAGCGCATGAGCGAAATCGCCCTGAAGCAGCCTGGCGTGGCCGACTCGGTAGCCTTCCCGGGCCTGTCGATCAACGGTTTCACCAACAGCCCCAACAGCGGC  >  NC_002947/3879834‑3880111

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: