Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,896 G→A 6.3% intergenic (+33/‑294) pcaK → / → pcaF‑I 4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,568,8960GA 6.3% 187.3 / 2.1 64intergenic (+33/‑294)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (31/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.48e-01

CCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTC  >  NC_002947/1568755‑1569038
                                                                                                                                               |                                                                                                                                              
cctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAcgcc                                                                                                                                         >  3:51380/1‑149 (MQ=255)
       tGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGcccctg                                                                                                                                  >  8:58441/1‑147 (MQ=255)
         gAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGc                                                                                                                                 <  7:96305/148‑1 (MQ=255)
            cTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAggg                                                                                                                             <  6:266433/149‑1 (MQ=255)
             ttCGAGCAGGTGCTGACGTCGTTGCTGGTACCGGCGTCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTTGTCAGCCACGCGGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCTCAgggg                                                                                                                            <  2:241595/149‑1 (MQ=255)
                gAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCg                                                                                                                                  <  4:141022/140‑1 (MQ=255)
                gAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCg                                                                                                                                  >  3:141022/1‑140 (MQ=255)
                gAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGagcg                                                                                                                          <  7:58441/148‑1 (MQ=255)
                gAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAAAGCCGACCCTGCGGGagcgg                                                                                                                         >  7:5000/1‑149 (MQ=255)
                      gTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgt                                                                                                                   >  2:90982/1‑149 (MQ=255)
                         cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGttacgct                                                                                                                <  4:184725/149‑6 (MQ=38)
                                        gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                  <  7:223653/132‑1 (MQ=255)
                                        gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                  >  8:223653/1‑132 (MQ=255)
                                        gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                  >  1:120862/1‑148 (MQ=255)
                                         tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                  >  5:299003/1‑147 (MQ=255)
                                           ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCCAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAg                                                                                               >  1:206115/1‑149 (MQ=255)
                                               cggcgCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGcc                                                                                          <  2:120862/149‑1 (MQ=255)
                                                 gcgcTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCt                                                                                        <  8:164198/149‑1 (MQ=255)
                                                 gcgcTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCt                                                                                        <  5:194769/149‑1 (MQ=255)
                                                     tGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgcgc                                                                                    >  7:276853/1‑149 (MQ=255)
                                                        cGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTg                                                                                 <  1:165028/149‑1 (MQ=255)
                                                          acggTGGGAGTGGTCTT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGCCGCACAGTTGCCCTGCGCGTggg                                                                               <  5:175961/146‑1 (MQ=255)
                                                              tGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                >  7:193015/1‑144 (MQ=255)
                                                              tGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                <  8:193015/144‑1 (MQ=255)
                                                               gggAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                          >  8:61092/1‑149 (MQ=255)
                                                                ggAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                             <  5:126626/145‑1 (MQ=255)
                                                                ggAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                             >  6:126626/1‑145 (MQ=255)
                                                                ggAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                          <  7:61092/148‑1 (MQ=255)
                                                                 gAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGACGGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                             >  1:53688/1‑145 (MQ=255)
                                                                 gAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                             <  2:53688/144‑1 (MQ=255)
                                                                      gtcgt‑cAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                   >  7:170351/1‑149 (MQ=255)
                                                                         gt‑cAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGc                                                                 >  6:7128/1‑148 (MQ=255)
                                                                               ggggCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAggg                                                                                                      <  1:293109/106‑1 (MQ=255)
                                                                               ggggCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAggg                                                                                                      >  2:293109/1‑106 (MQ=255)
                                                                                  gCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATgg                                                        >  4:137847/1‑149 (MQ=255)
                                                                                     ggTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGAACAGTTTCCCTGCGCGTGGGCTGGATAGTCGTGCAACAAATGGTTc                                                     >  4:313031/1‑149 (MQ=255)
                                                                                      ggcAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTCGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                                     <  1:194530/146‑1 (MQ=255)
                                                                                               gCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                           <  1:310919/149‑1 (MQ=255)
                                                                                               gCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                           >  4:91982/1‑149 (MQ=255)
                                                                                                 cGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCTTTAATCGcac                                          >  8:137929/1‑148 (MQ=255)
                                                                                                 cGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                           >  6:160275/1‑147 (MQ=255)
                                                                                                      gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATgg                                                        >  5:279860/1‑130 (MQ=255)
                                                                                                      gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATgg                                                        <  6:279860/130‑1 (MQ=255)
                                                                                                      gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGCTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATa                                     >  2:76642/1‑149 (MQ=255)
                                                                                                            tGAGGGCGAGGGCCGCGTTGCGGCCCATTAAAAGCCTAACGCCGCTCCCCCAGGGGCGGCGTTTTGCTGCAATAAGGGCGCACAGTTGACCTGCGCCTGGGGTGGATAGTCGTGCAAAAAATAGTTAGTTAATCGCACAAATAGTCGAt                               >  2:136069/1‑149 (MQ=255)
                                                                                                             gagggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                      >  6:150009/1‑125 (MQ=255)
                                                                                                             gagggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                      <  5:150009/125‑1 (MQ=255)
                                                                                                             gagggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAt                               >  7:60269/1‑148 (MQ=255)
                                                                                                               gggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTa                             >  5:210046/1‑148 (MQ=255)
                                                                                                               gggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTAt                            <  6:299003/149‑1 (MQ=255)
                                                                                                                ggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATc                           <  3:313031/149‑1 (MQ=255)
                                                                                                                ggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATc                           <  4:51380/149‑1 (MQ=255)
                                                                                                                ggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATc                           <  2:206115/149‑1 (MQ=255)
                                                                                                                   gacggcCGCGCTGCGGCCCATTCACAGCGTACCGCCGCACCCGCAGGGGCGGCGCTGTGCTGCAATCACGGCGCACAGTTGCCCTGCGCGTGGCGTGGATACTCGTGCAACAAATGGTGCGTTAATCGCACAGATAGTCGAATATCGGa                        <  8:276853/146‑1 (MQ=255)
                                                                                                                     gggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAtt                              <  4:145806/141‑1 (MQ=255)
                                                                                                                     gggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAtt                              >  3:145806/1‑141 (MQ=255)
                                                                                                                     gggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt                      >  7:148432/1‑149 (MQ=255)
                                                                                                                          gcgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAg                 <  5:160275/149‑1 (MQ=255)
                                                                                                                               gCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtac             <  5:7128/148‑1 (MQ=255)
                                                                                                                                    ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCTCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCa       <  1:136069/149‑1 (MQ=255)
                                                                                                                                         cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtctc  >  6:379/1‑149 (MQ=255)
                                                                                                                                         cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtct   <  8:60269/148‑1 (MQ=255)
                                                                                                                                         cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtct   <  3:91982/148‑1 (MQ=255)
                                                                                                                                          acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGGCGATTATCGGATTGTAAGGTACCTACCAGtctc  <  2:265751/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGT‑CAAGGGGCTGGTCAGCCACGCCGAC‑GCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTC  >  NC_002947/1568755‑1569038

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