Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,915 G→A 8.0% intergenic (+52/‑275) 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,9150GA 8.0% 177.3 / 4.1 62intergenic (+52/‑275)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/3);  total (29/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 = 8.34e-01

TCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACC  >  NC_002947/1568769‑1569063
                                                                                                                                                  |                                                                                                                                                    
tCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggc                                                                                                                                                    >  1:170696/1‑149 (MQ=255)
 cGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggc                                                                                                                                                    <  1:327968/148‑1 (MQ=255)
     cAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGt                                                                                                                                               >  1:23483/1‑149 (MQ=255)
        gtggtgTCGGGGTTGTTGGTACGGGCGGCGCTGGCGACCGTGGGAGTGGTGGTCAAGGGGGTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtg                                                                                                                                             <  8:287322/148‑1 (MQ=255)
         tGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                           <  4:275049/149‑1 (MQ=255)
          gCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGttacgc                                                                                                                                          <  4:286628/149‑5 (MQ=38)
            tGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGtt                                                                                                                                              <  8:305254/143‑1 (MQ=255)
            tGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGtt                                                                                                                                              >  7:305254/1‑143 (MQ=255)
            tGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctg                                                                                                                                        >  3:2150/1‑149 (MQ=255)
             gACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctg                                                                                                                                        >  3:150334/1‑148 (MQ=255)
              aCGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGt                                                                                                                                               <  6:291404/140‑1 (MQ=255)
              aCGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGt                                                                                                                                               >  5:291404/1‑140 (MQ=255)
                ggCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAAt                                                                                                                                    >  5:6256/1‑149 (MQ=255)
                        tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACTCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgcg                                                                                                                             >  5:246701/1‑148 (MQ=255)
                           tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcaca                                                                                                                         <  4:2150/149‑1 (MQ=255)
                             ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgc                                                                                                                              <  1:154192/142‑1 (MQ=255)
                             ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgc                                                                                                                              >  2:154192/1‑142 (MQ=255)
                                           gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                         <  6:6256/149‑1 (MQ=255)
                                           gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGtta                                                                                                                                             <  7:103822/113‑2 (MQ=255)
                                           gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGtta                                                                                                                                             >  8:103822/1‑112 (MQ=255)
                                            aCCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTggg                                                                                                        <  4:150334/149‑1 (MQ=255)
                                              cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                                                      >  1:344659/1‑149 (MQ=255)
                                              cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGt                                                                                                      <  2:23483/149‑1 (MQ=255)
                                               gTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTg                                                                                                     <  8:309503/149‑1 (MQ=255)
                                                  ggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAg                                                                                                                        <  1:27028/127‑1 (MQ=255)
                                                  ggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAg                                                                                                                        >  2:27028/1‑127 (MQ=255)
                                                            tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                                            <  8:298368/145‑1 (MQ=255)
                                                            tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                                            >  7:298368/1‑145 (MQ=255)
                                                            tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTg                                                                                           >  1:102869/1‑146 (MQ=255)
                                                            tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTg                                                                                           <  2:102869/146‑1 (MQ=255)
                                                                ggggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACaaa                                                                                    <  2:170696/149‑1 (MQ=255)
                                                                 gggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAAt                                                                                   >  1:24837/1‑149 (MQ=255)
                                                                   gCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTg                                                                                           <  7:208818/139‑1 (MQ=255)
                                                                   gCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTg                                                                                           >  8:208818/1‑139 (MQ=255)
                                                                     tGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                                            >  3:96533/1‑136 (MQ=255)
                                                                     tGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                                            <  4:96533/136‑1 (MQ=255)
                                                                       gTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                             <  3:158219/149‑1 (MQ=255)
                                                                       gTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                             >  2:225301/1‑149 (MQ=255)
                                                                       gTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                             >  6:156748/1‑149 (MQ=255)
                                                                         cAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                                   >  5:258019/1‑125 (MQ=255)
                                                                         cAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                                   <  6:258019/125‑1 (MQ=255)
                                                                            ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                               <  5:103335/142‑1 (MQ=255)
                                                                            ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                               >  6:103335/1‑142 (MQ=255)
                                                                              aCGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTggg                                                                                                        <  2:272271/115‑1 (MQ=255)
                                                                              aCGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTggg                                                                                                        >  1:272271/1‑115 (MQ=255)
                                                                               cGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGc                                                                     <  3:76049/149‑1 (MQ=255)
                                                                                gCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                                                    <  2:344659/149‑1 (MQ=255)
                                                                                            tGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAt                                                        >  6:125312/1‑149 (MQ=255)
                                                                                                   gagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGa                                                 <  2:24837/149‑1 (MQ=255)
                                                                                                        ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTa                                            <  5:156748/149‑1 (MQ=255)
                                                                                                            gTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGt                                        <  5:125312/149‑1 (MQ=255)
                                                                                                                cGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacc                                     >  3:9411/1‑148 (MQ=255)
                                                                                                                    ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCa                                <  1:297701/149‑1 (MQ=255)
                                                                                                                     cATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAg                               >  3:20078/1‑149 (MQ=255)
                                                                                                                      aTTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGt                                                            >  7:174442/1‑119 (MQ=255)
                                                                                                                      aTTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGt                                                            <  8:174442/119‑1 (MQ=255)
                                                                                                                         cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaac                                                                                       <  4:179570/89‑1 (MQ=255)
                                                                                                                          acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTcc                          <  4:9411/149‑1 (MQ=255)
                                                                                                                                   cgccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                 <  1:225301/149‑1 (MQ=255)
                                                                                                                                    gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtctc                           >  8:99026/1‑138 (MQ=255)
                                                                                                                                    gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtctc                           <  7:99026/138‑1 (MQ=255)
                                                                                                                                             gCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCg         >  5:257036/1‑147 (MQ=255)
                                                                                                                                             gCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCg         >  1:70259/1‑147 (MQ=255)
                                                                                                                                             gCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCg         >  8:315461/1‑147 (MQ=255)
                                                                                                                                                 gggCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAc   >  5:18408/1‑149 (MQ=255)
                                                                                                                                                  ggcggcGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAcc  >  3:122532/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACC  >  NC_002947/1568769‑1569063

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