Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,050,675 C→T 7.9% G236G (GGC→GGT oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,050,6750CT 7.9% 184.8 / 3.6 63G236G (GGC→GGToprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CAATTACCCCAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGT  >  NC_002947/4050528‑4050821
                                                                                                                                                    |                                                                                                                                                  
cAATTACCCCAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCa                                                                                                                                                   >  6:287593/1‑149 (MQ=255)
      ccccAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTgg                                                                                                                                             >  3:127534/1‑149 (MQ=255)
       cccAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGt                                                                                                                                            <  6:208484/149‑1 (MQ=255)
         cAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTAc                                                                                                                                          >  6:25492/1‑149 (MQ=255)
             gCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGTAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGt                                                                                                                                      <  1:109947/149‑1 (MQ=255)
              cccACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGtg                                                                                                                                     <  3:157419/149‑2 (MQ=255)
                 aCCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTCGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCACCTCCAGCCGCGACTGGCACCTGGGGCCCAAGCACACCACCAGCACCGTGCTACCGGTAGaa                                                                                                                                  >  4:17850/1‑149 (MQ=255)
                    tCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTg                                                                                                                               >  8:13767/1‑149 (MQ=255)
                      ggCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTgg                                                                                                                              <  5:25492/148‑1 (MQ=255)
                           gcgcGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCATCTCCAGCCGCGCCTCGGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGtggtgt                                                                                                                                 <  1:271424/141‑6 (MQ=255)
                           gcgcGTGAAGTACGACCTGTCACCTGCCTCGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCAGAAAGCAACGGCAACTCCAGCCGCGCCTCGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGtggtgt                                                                                                                                 >  2:271424/1‑135 (MQ=255)
                             gcgTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAg                                                                                                                      <  5:287593/149‑1 (MQ=255)
                               gTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCt                                                                                                                    <  7:13767/149‑1 (MQ=255)
                                  aaGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGc                                                                                                                  <  6:43902/148‑1 (MQ=255)
                                  aaGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGc                                                                                                                  <  4:127534/148‑1 (MQ=255)
                                       cGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGt                                                                                                             <  2:308517/148‑1 (MQ=255)
                                          ccTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAAc                                                                                                         >  6:370268/1‑149 (MQ=255)
                                           cTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGt                                                                                                             <  2:8435/144‑1 (MQ=255)
                                           cTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGt                                                                                                             >  1:8435/1‑144 (MQ=255)
                                            tGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTg                                                                                                       >  6:191908/1‑149 (MQ=255)
                                              tccacggcctgTCAGTTGCGCGTGGCGGCGTTCCACGTCGACCCCGACAGCAACGGCAACTCCAGCCGCGCCTGACACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACGGGTAGAACTGGTATACAAGCTGCAGGTTGaa                                                                                                          <  3:238932/138‑1 (MQ=255)
                                              tCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGCCACCTGGGGCCCCAGCACACCACCGTCACCCTGGTACCGCTAGACCTGGAATACAAGCTGCAGGGTGaa                                                                                                          >  4:238932/1‑144 (MQ=255)
                                                aCCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAg                                                                                                   >  5:178524/1‑149 (MQ=255)
                                                  ctgcctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGc                                                                                                 <  3:304990/149‑1 (MQ=255)
                                                  ctgcctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGc                                                                                                 <  6:178524/149‑1 (MQ=255)
                                                    gcctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGc                                                                                                 >  4:292834/1‑147 (MQ=255)
                                                     cctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGa                                                                                               <  3:46150/148‑1 (MQ=255)
                                                      ctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAg                                                                                              >  4:240893/1‑148 (MQ=255)
                                                            gTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGc                                                                                       <  5:370268/149‑1 (MQ=255)
                                                                  cgTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGtt                                                                                 <  5:191908/149‑1 (MQ=255)
                                                                      gcggcgTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtacta                                                                             >  6:304901/1‑149 (MQ=255)
                                                                          cgTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactactac                                                                         >  1:208965/1‑149 (MQ=255)
                                                                                 cGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGct                                                                  <  7:155209/149‑1 (MQ=255)
                                                                                      aCCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAt                                                             <  2:208965/149‑1 (MQ=255)
                                                                                         ccGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAtgt                                                           >  5:107988/1‑148 (MQ=255)
                                                                                          cGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAtgtg                                                          >  2:314401/1‑148 (MQ=255)
                                                                                                  aCGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactactac                                                                         >  5:287234/1‑125 (MQ=255)
                                                                                                  aCGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtactactac                                                                         <  6:287234/125‑1 (MQ=255)
                                                                                                     gCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGctc                                                                 <  1:145768/130‑1 (MQ=255)
                                                                                                       aaCTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAg                                            >  1:147254/1‑149 (MQ=255)
                                                                                                            cAGCCGCTCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgat                                        <  3:240893/148‑1 (MQ=255)
                                                                                                               ccGCGCGTGGCACCTGGGGCCCAAGCACACCACCAGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGatgatg                                    <  3:17850/149‑1 (MQ=255)
                                                                                                               ccGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGatgatg                                    <  6:130893/149‑1 (MQ=255)
                                                                                                                         gCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTc                             <  3:292834/147‑1 (MQ=255)
                                                                                                                                   cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatg                                       <  4:321541/127‑1 (MQ=255)
                                                                                                                                   cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatg                                       <  4:316167/127‑1 (MQ=255)
                                                                                                                                   cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatg                                       >  3:321541/1‑127 (MQ=255)
                                                                                                                                   cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatg                                       >  3:316167/1‑127 (MQ=255)
                                                                                                                                    ccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGc                >  3:247953/1‑149 (MQ=255)
                                                                                                                                     cAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAGACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGc                >  7:157212/1‑148 (MQ=255)
                                                                                                                                        gCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAt             >  8:265261/1‑148 (MQ=255)
                                                                                                                                        cagttctaccGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAt             >  1:122222/8‑148 (MQ=255)
                                                                                                                                          acacCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCga                                         <  7:226621/118‑1 (MQ=255)
                                                                                                                                          acacCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCga                                         >  8:226621/1‑118 (MQ=255)
                                                                                                                                            acctccGGCACCGTGGTACCGGTAGATCTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAg                                 <  2:202884/124‑1 (MQ=255)
                                                                                                                                            accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACACGCTTGGTTACTACTACGACAGCTCTGATGTGAAACCCATCGGCAGCGATGATGAAGTGTCCGGACGAGGTGGCCATTACct        >  5:200230/1‑149 (MQ=255)
                                                                                                                                            accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAg                                 <  2:202888/124‑1 (MQ=255)
                                                                                                                                            accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAg                                 >  1:202888/1‑124 (MQ=255)
                                                                                                                                            accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGCGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACCACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAg                                 >  1:202884/1‑124 (MQ=255)
                                                                                                                                              caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAgt                                <  6:195797/123‑1 (MQ=255)
                                                                                                                                              caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAgt                                >  5:195797/1‑123 (MQ=255)
                                                                                                                                                 cGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTgg   <  7:265261/149‑1 (MQ=255)
                                                                                                                                                  ggCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGt  <  5:304901/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAATTACCCCAACGCCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGGT  >  NC_002947/4050528‑4050821

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