Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,514,648 C→T 6.7% A130T (GCA→ACA)  ftsK ← DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,514,6480CT 6.7% 177.4 / 2.8 60A130T (GCA→ACA) ftsKDNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base T (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGAAGTAACCCA  >  NC_002947/4514502‑4514794
                                                                                                                                                  |                                                                                                                                                  
gCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTgccgc                                                                                                                                                 >  2:987042/1‑150 (MQ=255)
   gaacaggaacagGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGcccc                                                                                                                                              <  1:1140520/150‑1 (MQ=255)
   gaacaggaacagGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGcccc                                                                                                                                              >  2:1259318/1‑150 (MQ=255)
    aacaggaacaggGCAATGAACATCACCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGccccc                                                                                                                                             <  1:661854/150‑1 (MQ=255)
       aggaacaggGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGAc                                                                                                                                          >  2:1627322/1‑150 (MQ=255)
        ggaacaggGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACa                                                                                                                                         <  1:1326136/150‑1 (MQ=255)
             aggGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCAcc                                                                                                                                    >  2:316767/1‑150 (MQ=255)
              gggCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCa                                                                                                                                   >  1:955760/1‑150 (MQ=255)
               ggCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAg                                                                                                                                  >  1:708057/1‑150 (MQ=255)
               ggCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAg                                                                                                                                  <  1:22897/150‑1 (MQ=255)
                     gAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAacac                                                                                                                            >  2:681169/1‑150 (MQ=255)
                          tCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGc                                                                                                                       <  1:1259318/150‑1 (MQ=255)
                             gCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCaa                                                                                                                    <  2:708057/150‑1 (MQ=255)
                               gTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATc                                                                                                                  <  2:24413/150‑1 (MQ=255)
                                   tgcCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAgcc                                                                                                              <  2:1107168/150‑1 (MQ=255)
                                   tgcCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAgcc                                                                                                              >  1:638165/1‑150 (MQ=255)
                                     cccTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAgccgc                                                                                                            >  1:1016360/1‑150 (MQ=255)
                                       cTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGTCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCa                                                                                                          >  2:330185/1‑150 (MQ=255)
                                        tGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAg                                                                                                         <  1:316767/150‑1 (MQ=255)
                                          cACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGTCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGga                                                                                                       <  1:330185/150‑1 (MQ=255)
                                          cACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAacac                                                                                                                            <  1:414393/129‑1 (MQ=255)
                                          cACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAacac                                                                                                                            >  2:414393/1‑129 (MQ=255)
                                            cGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGgaga                                                                                                     <  1:987042/150‑1 (MQ=255)
                                                 agcagGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAAcagc                                                                                                >  1:1290183/1‑150 (MQ=255)
                                                 agcagGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTAAAAGAAGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAAcagc                                                                                                >  1:691550/1‑150 (MQ=255)
                                                  gcagGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAAcagcc                                                                                               >  2:414496/1‑150 (MQ=255)
                                                     ggTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAAcagccagc                                                                                            <  2:955760/150‑1 (MQ=255)
                                                           gcgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAgg                                                                                                                        >  2:33580/1‑116 (MQ=255)
                                                           gcgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAgg                                                                                                                        <  1:33580/116‑1 (MQ=255)
                                                           gcgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGgag                                                                                                      <  1:600648/134‑1 (MQ=255)
                                                           gcgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGgag                                                                                                      >  2:600648/1‑134 (MQ=255)
                                                            cgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCa                                                                                     <  1:213102/150‑1 (MQ=255)
                                                            cgcCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCa                                                                                     <  1:414496/150‑1 (MQ=255)
                                                                     gCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACgg                                                                            >  1:270818/1‑150 (MQ=255)
                                                                          ggCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCtggt                                                                       <  2:1290183/150‑1 (MQ=255)
                                                                             tttCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGTCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGc                                                                    >  2:440105/1‑150 (MQ=255)
                                                                             tttCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGTCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGc                                                                    <  1:440105/150‑1 (MQ=255)
                                                                                  cccAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAg                                                                                                         >  1:1514678/1‑108 (MQ=255)
                                                                                  cccAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAg                                                                                                         <  2:1514678/108‑1 (MQ=255)
                                                                                            cACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAAt                                                       <  1:1536241/148‑1 (MQ=255)
                                                                                            cACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAAt                                                       >  2:1536241/1‑148 (MQ=255)
                                                                                            cACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCt                                                     <  1:681169/150‑1 (MQ=255)
                                                                                                      gaACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGtttt                                           >  2:707303/1‑150 (MQ=255)
                                                                                                      gaACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGtttt                                           <  1:1168863/150‑1 (MQ=255)
                                                                                                          ggCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGAcc                                       <  1:1623378/150‑1 (MQ=255)
                                                                                                           gCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCa                                                                                                                 <  1:1143427/75‑1 (MQ=255)
                                                                                                           gCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCa                                                                                                                 >  2:1143427/1‑75 (MQ=255)
                                                                                                                               aaGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAata                  >  1:639339/1‑150 (MQ=255)
                                                                                                                               aaGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAata                  <  1:1627322/150‑1 (MQ=255)
                                                                                                                                  tGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAg               >  2:1463118/1‑150 (MQ=255)
                                                                                                                                    gATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGga                                                                                                       <  2:1109766/60‑1 (MQ=255)
                                                                                                                                    gATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGga                                                                                                       >  1:1109766/1‑60 (MQ=255)
                                                                                                                                    gATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGc             >  1:787258/1‑150 (MQ=255)
                                                                                                                                       atGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCAAGCAGCAGCGGGAAAATATAGGCGaa          <  1:1345502/150‑1 (MQ=255)
                                                                                                                                          cgcCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGaa          >  2:101406/1‑147 (MQ=255)
                                                                                                                                          cgcCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGaa          <  1:101406/147‑1 (MQ=255)
                                                                                                                                           gcCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGaaa                     >  2:1403846/1‑135 (MQ=255)
                                                                                                                                           gcCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGaaa                     <  1:1403846/135‑1 (MQ=255)
                                                                                                                                           gcCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGAAGTaa      >  1:1308787/1‑150 (MQ=255)
                                                                                                                                               gTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGAAGTAACCCa  <  2:691550/150‑1 (MQ=255)
                                                                                                                                               gTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGAAGTAACCCa  >  1:441465/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGCGCCAGTGCCGCCCCCGACAGCACCAGGAACACCAGGCCAATCAGCCGCCAGGAGAACAGCCAGCCGCTCCAGTCCCACGGCTGGTGGCGTTCGCGGAAAATCTGCCAGGTTTTGACCGCCAGCAGCAGCGGGAAAATATAGGCGAAGTAACCCA  >  NC_002947/4514502‑4514794

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