Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,173,251 G→A 8.3% intergenic (‑292/+40) fadD‑I ← / ← fadD‑II long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,173,2510GA 8.3% 166.2 / 4.8 60intergenic (‑292/+40)fadD‑I/fadD‑IIlong‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

AGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTT  >  NC_002947/5173108‑5173395
                                                                                                                                               |                                                                                                                                                
aGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAg                                                                                                                                             <  7:50915/149‑1 (MQ=255)
  gCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGa                                                                                                                                           <  1:189649/149‑1 (MQ=255)
    cTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                         <  8:8450/149‑1 (MQ=255)
     ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCg                                                                                                                                        >  7:281147/1‑149 (MQ=255)
        tGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAATGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCggg                                                                                                                                      <  6:365329/148‑1 (MQ=38)
        tGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCggg                                                                                                                                      <  6:243959/148‑1 (MQ=255)
        tGGTCAAAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCGGGt                                                                                                                                     <  7:316418/149‑1 (MQ=38)
          gTCACAAAAATGACTCAATGAACAATACAGGTCAGTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGtt                                                                                                                                    >  1:37013/1‑148 (MQ=255)
               aaaaaTGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcg                                                                                                                              >  3:330902/1‑149 (MQ=255)
                   atcccccAACGAACAACACAGCCCATCTTTACCCAAATAATCCACAATCCCCCAGATAAGCCTGTCTAGAGCCTGCTACACCGGTTTAAGGCATCCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGCAGCGGGTTTACCCGCgaa                                                                                                                          <  8:281147/143‑1 (MQ=16)
                      aCTCAATCAACAATACCGGTCATTTTAAGAGTAATAATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAACGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                       <  4:317547/149‑1 (MQ=38)
                          aaTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGGCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGc                                                                                                                     <  1:13571/147‑1 (MQ=255)
                              aaCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATcc                                                                                                                <  7:352774/148‑1 (MQ=255)
                              aaCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTACTAAATTGGGCAGTGTAGGAGCGGGTTTTCCCGCGAAGAAGCCATacc                                                                                                               >  5:26627/1‑149 (MQ=255)
                                 aaTACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGGGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGt                                                                                                            >  6:7015/1‑149 (MQ=255)
                                   tACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTa                                                                                                           >  3:53501/1‑148 (MQ=255)
                                      aGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTAAATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCAc                                                                                                       >  3:207036/1‑149 (MQ=255)
                                          cATTTTTATAGTAATTATCAAAAATGTCGCATTAAAGCCTGTATCGAGCGTTTTACAGCGTTTTACGCCATTCGCGGTTAACCCGGCCCCTACATTGGGCTGTTTCGGCGGGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCt                                                                                                    <  6:26627/148‑1 (MQ=255)
                                                  gagTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgcg                                                                                            <  5:197874/148‑1 (MQ=255)
                                                      aaTTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCACCt                                                                                       >  2:6720/1‑149 (MQ=255)
                                                           tCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATcc                                                                                                                >  3:70450/1‑119 (MQ=255)
                                                           tCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATcc                                                                                                                <  4:70450/119‑1 (MQ=255)
                                                           tCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATAGCGCAACTCccgc                                                                                  >  3:65175/1‑149 (MQ=255)
                                                            ccAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCccgcc                                                                                 >  5:101487/1‑149 (MQ=255)
                                                              aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGTTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCcgccgc                                                                               <  4:379838/149‑1 (MQ=255)
                                                              aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCccgccg                                                                                >  2:120448/1‑148 (MQ=255)
                                                                 aTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGGGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCaa                                                                             >  8:177128/1‑148 (MQ=255)
                                                                  tGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCaaaa                                                                           <  4:207036/149‑1 (MQ=255)
                                                                          aTAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCGCCGACAATGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATAGCGGAACAGCCGCCGCAAAATAGTGccc                                                                   >  8:312252/1‑149 (MQ=38)
                                                                             aGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCAc                                                                 <  2:37013/148‑1 (MQ=255)
                                                                              ggCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCg                                                               >  1:333316/1‑149 (MQ=255)
                                                                               gCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATAGCGCAACTCCCGCCGCAAAATCTTGCCCACCgg                                                              >  5:137651/1‑149 (MQ=255)
                                                                                cTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                             <  3:314968/149‑1 (MQ=255)
                                                                                  gTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGt                                                           >  5:169690/1‑149 (MQ=255)
                                                                                       gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATgg                                                      <  1:152161/149‑1 (MQ=255)
                                                                                        agCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                                     <  6:101487/149‑1 (MQ=255)
                                                                                           cTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCACCTACAATGGGGAGTGTAGGAGCGGGGTTAACCGCGAAGAAACCCTCCCGGTATAAAGCGAAAAGGCGCAACTCCCCCCGCAAAATCTTGCCCACCGCCGGCCTGGGCa                                                   >  8:202853/1‑148 (MQ=255)
                                                                                               tAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGa                                               <  4:53501/148‑1 (MQ=255)
                                                                                                       tttAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGGTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATAGCGCAACTACCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGAATCGCgcag                                      >  3:199025/1‑149 (MQ=255)
                                                                                                           aGGCATTCGCGGGTAAACCCGCTCCGACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATAACGCGATATCGCGCAACTACCGCCGCAAAATCTTGCACACCGGCGTCAAGGGCAACGCTTAGAGCAGAAc                                   >  4:309391/1‑148 (MQ=37)
                                                                                                             gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAt                                <  4:330902/149‑1 (MQ=255)
                                                                                                                   gcgGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                          <  6:169690/149‑1 (MQ=255)
                                                                                                                     gggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTgg                        <  5:7015/149‑1 (MQ=255)
                                                                                                                      ggTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTggg                       >  3:191367/1‑149 (MQ=255)
                                                                                                                          aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCTTCCCGGTATCACGCTATATCGCGCAACTGCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGCTGTGCTGGGGTACt                   >  8:245144/1‑149 (MQ=255)
                                                                                                                          aaCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCCGCACGATGTGCTTGGGTACt                   >  8:245145/1‑149 (MQ=255)
                                                                                                                           aCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                  >  7:185581/1‑149 (MQ=255)
                                                                                                                            cccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACttt                 >  4:2786/1‑149 (MQ=255)
                                                                                                                            cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                  <  2:197330/148‑1 (MQ=255)
                                                                                                                            cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt                  >  1:197330/1‑148 (MQ=255)
                                                                                                                               gCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg                >  7:275010/1‑147 (MQ=255)
                                                                                                                                   cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGcc           <  6:137651/148‑1 (MQ=255)
                                                                                                                                    tACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCgg         >  5:293661/1‑149 (MQ=255)
                                                                                                                                       aTTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTg       >  2:80069/1‑148 (MQ=255)
                                                                                                                                          gggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAg    <  7:312252/148‑1 (MQ=255)
                                                                                                                                           ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGTTGAAGtt  <  5:97588/149‑1 (MQ=255)
                                                                                                                                           ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGtt  >  6:97588/1‑149 (MQ=255)
                                                                                                                                           ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGtt  >  6:204736/1‑149 (MQ=255)
                                                                                                                                             cAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGcc           >  8:14891/1‑138 (MQ=255)
                                                                                                                                             cAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGcc           <  7:14891/138‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTT  >  NC_002947/5173108‑5173395

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