Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,598,837 C→T 8.6% intergenic (‑76/+55) dctB ← / ← bglX C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,598,8370CT 8.6% 169.7 / 3.4 58intergenic (‑76/+55)dctB/bglXC4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (31/22);  new base T (3/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAACTGG  >  NC_002947/1598704‑1598985
                                                                                                                                      |                                                                                                                                                    
ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                       >  6:226114/1‑149 (MQ=255)
  aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGttt                                                                                                                                     >  7:4888/1‑149 (MQ=255)
     aCGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAc                                                                                                                                  <  2:148697/149‑1 (MQ=255)
       gATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAccc                                                                                                                                >  4:52413/1‑149 (MQ=255)
         tCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcgc                                                                                                                              <  8:180834/149‑1 (MQ=255)
            gcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGa                                                                                                                            >  7:17801/1‑148 (MQ=255)
                    agcAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                                    <  8:243082/146‑1 (MQ=255)
                     tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTAc                                                                                                                                  >  2:264359/1‑133 (MQ=21)
                     tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTAc                                                                                                                                  <  1:264359/133‑1 (MQ=21)
                     tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgg                                                                                                                   >  2:24861/1‑148 (MQ=255)
                     tCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgg                                                                                                                   >  3:219797/1‑148 (MQ=255)
                       aGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAAc                                                                                                                <  3:52413/149‑1 (MQ=255)
                          cccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAgg                                                                                                             >  5:122033/1‑149 (MQ=255)
                           ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                       <  7:29367/122‑1 (MQ=255)
                           ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                       >  8:29367/1‑122 (MQ=255)
                           ccAGACACAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                            >  5:202908/1‑149 (MQ=255)
                            cAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGaa                                                                                                                 >  8:120739/1‑143 (MQ=255)
                            cAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGaa                                                                                                                 <  7:120739/143‑1 (MQ=255)
                             aGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                          <  5:226114/149‑1 (MQ=255)
                             aGACACAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCTCGAAAAGGCCGGAACAggcgg                                                                                                          <  6:202908/149‑1 (MQ=255)
                                   aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGTCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                              >  1:205570/1‑107 (MQ=25)
                                   aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                              <  2:205570/107‑1 (MQ=25)
                                      gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                                 <  7:78427/149‑1 (MQ=255)
                                      gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                                 <  4:219797/149‑1 (MQ=255)
                                          gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTc                                                                                              >  7:197819/1‑148 (MQ=255)
                                            aaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                           <  8:197819/149‑1 (MQ=255)
                                             aGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAg                                                                                          <  8:17801/149‑1 (MQ=255)
                                                aaaaGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACGCTGCAGCACCTGTAGGAGCGGTTTTACCCGCGAAAAGGCCGGAACAGCCGGCATTTTTCAGTCAAAGaaa                                                                                       >  7:244153/1‑149 (MQ=38)
                                                 aaaGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGa                                                                                                                  >  1:107720/1‑121 (MQ=255)
                                                       cATCGGGGGCGGATCCGGTGGTGGGATTGTCGGGGCCTGCCTCTTCCCGGGTAAACCCGCGCATACAGGTAATCTGCAGCCCCTTTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAGCCGGCATTTTTAAAACAAACCAAtaaa                                                                                   >  1:148586/1‑144 (MQ=255)
                                                       cATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTCTGCCGCCCCGGGAGGCGCGGGGTTACCCGCGCAGAGGCAggcc                                                                                                                 >  7:176797/1‑114 (MQ=25)
                                                         tCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTGTTTCCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGc                                                                               >  6:131869/1‑148 (MQ=255)
                                                                            tggGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCTCCTGTAGGCGCGGGTTTAACCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTAAAAGCTCTGCTGCTGCCCCGCCTgg                                                           >  4:138878/1‑149 (MQ=38)
                                                                            tggGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                           <  6:122033/149‑1 (MQ=255)
                                                                              ggA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTggg                                                          >  5:85110/1‑148 (MQ=255)
                                                                                           ctgactCTTCGCGGGTAATCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCATTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAc                                             <  8:244153/149‑1 (MQ=255)
                                                                                             gcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAg                                                                                          >  1:44626/1‑102 (MQ=255)
                                                                                             gcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAg                                                                                          <  2:44626/102‑1 (MQ=255)
                                                                                              cctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTg                                          >  8:200615/1‑149 (MQ=255)
                                                                                                     gcgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCATCGCCTGGGAATCCAGCCCTACCTGGACAtt                                    <  1:24861/148‑1 (MQ=255)
                                                                                                         gTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACt                               <  7:70876/149‑1 (MQ=255)
                                                                                                         gTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTCGCAATCCAGCCCTACCTGGACATTGAAc                                >  6:125670/1‑148 (MQ=255)
                                                                                                            aaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                            >  1:60469/1‑149 (MQ=255)
                                                                                                               ccGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAg                                                   <  3:180331/123‑1 (MQ=255)
                                                                                                               ccGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAg                                                   >  4:180331/1‑123 (MQ=255)
                                                                                                               ccGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTg                         >  3:67454/1‑149 (MQ=255)
                                                                                                                 gcgcCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                         <  4:164738/115‑1 (MQ=255)
                                                                                                                 gcgcCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                         >  3:164738/1‑115 (MQ=255)
                                                                                                                   gcCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTc                     <  6:85110/149‑1 (MQ=255)
                                                                                                                    ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCg                    >  7:134776/1‑149 (MQ=37)
                                                                                                                     cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg                   <  2:60469/149‑1 (MQ=255)
                                                                                                                     cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg                   >  8:74465/1‑149 (MQ=255)
                                                                                                                           gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAg                                                                                >  2:256847/1‑82 (MQ=255)
                                                                                                                           gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAg                                                                                <  1:256847/82‑1 (MQ=255)
                                                                                                                           gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAccc             >  1:13079/1‑149 (MQ=255)
                                                                                                                           gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAccc             >  7:253542/1‑149 (MQ=255)
                                                                                                                                 tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc        >  5:62865/1‑148 (MQ=255)
                                                                                                                                 tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc        >  5:2040/1‑148 (MQ=255)
                                                                                                                                     gcaCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGGTCGGCAACCCGCTGCAACTg   >  2:79712/1‑149 (MQ=255)
                                                                                                                                      caCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAACTgg  <  4:67454/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                    
CCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGA‑TGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAACTGG  >  NC_002947/1598704‑1598985

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: ‑