Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,912 A→G 6.7% intergenic (+49/‑278) 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,9120AG 6.7% 218.7 / 5.1 74intergenic (+49/‑278)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (2/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

GCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAG  >  NC_002947/1568766‑1569059
                                                                                                                                                  |                                                                                                                                                   
gCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgg                                                                                                                                                   >  3:328287/1‑149 (MQ=255)
gCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgg                                                                                                                                                   >  8:110545/1‑149 (MQ=255)
   tCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGagc                                                                                                                                                >  2:255636/1‑149 (MQ=255)
    cGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgggg                                                                                                                                                 <  7:110545/147‑1 (MQ=255)
        cAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGGCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGt                                                                                                                                           <  3:254882/149‑1 (MQ=255)
         aGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGtt                                                                                                                                          <  4:262268/149‑1 (MQ=255)
          ggTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtg                                                                                                                                         <  4:312613/149‑1 (MQ=255)
             gCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgc                                                                                                                                      <  1:356154/149‑1 (MQ=255)
             gCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACAGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGTCGTTGtgc                                                                                                                                      >  4:235036/1‑149 (MQ=255)
              cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGttacgct                                                                                                                                     <  1:255636/149‑6 (MQ=38)
              cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGttacgct                                                                                                                                     <  5:297773/149‑6 (MQ=38)
                     cGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAAt                                                                                                                                >  4:78871/1‑147 (MQ=255)
                     cGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCa                                                                                                                              <  4:328287/149‑1 (MQ=255)
                              tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                       <  1:45745/131‑1 (MQ=255)
                              tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                       >  2:45745/1‑131 (MQ=255)
                                ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAg                                                                                                                    <  4:242580/148‑1 (MQ=255)
                                  ggcggcGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTg                                                                                                                 <  3:78871/149‑1 (MQ=255)
                                        gcTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgc                                                                                                           <  8:19487/149‑1 (MQ=255)
                                              gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                     >  2:286219/1‑149 (MQ=255)
                                                           gtcgtcAAGGGGCTGGTCAGCCACGCCGCCGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACCGCGTAACGCCGCTCCCGGAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGGGTGGGGTGGATTGTCGt                                                                                        >  6:230243/1‑149 (MQ=255)
                                                              gtcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGc                                                                                      >  8:181201/1‑148 (MQ=255)
                                                               tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                               >  4:330657/1‑138 (MQ=255)
                                                               tcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                               <  3:330657/138‑1 (MQ=255)
                                                                cAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGccc                                                                                                              >  4:78339/1‑122 (MQ=255)
                                                                cAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGccc                                                                                                              <  3:78339/122‑1 (MQ=255)
                                                                 aaGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaac                                                                                   <  7:181201/148‑1 (MQ=255)
                                                                   ggggCTTGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACaaa                                                                                <  3:235036/149‑1 (MQ=255)
                                                                    gggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGt                                                                                        <  6:275173/140‑1 (MQ=255)
                                                                    gggCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTgagctcgg                                                                                >  5:275173/1‑141 (MQ=255)
                                                                     gtcTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGGGTTGCGGACCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATg                                                                              <  5:330685/147‑1 (MQ=255)
                                                                         ggTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                     >  2:128131/1‑122 (MQ=255)
                                                                         ggTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                     <  1:128131/122‑1 (MQ=255)
                                                                            cAGCCACGCCGACGGACACTGAGTGCGAGGGCCGCGTTGCGGCCCATTCACCGCGTTACGCCTCTACCTCAGTGGCGGCGTTGTGCTTCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                        <  4:205026/148‑1 (MQ=255)
                                                                               ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAAt                                                                    >  5:289770/1‑149 (MQ=255)
                                                                               ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAAt                                                                    >  2:305758/1‑149 (MQ=255)
                                                                                cACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                                                   <  5:360272/149‑1 (MQ=255)
                                                                                  cGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCg                                                                  >  3:167899/1‑148 (MQ=255)
                                                                                   gCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                                                >  3:13789/1‑149 (MQ=255)
                                                                                        cGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGa                                                            <  6:289770/148‑1 (MQ=255)
                                                                                           aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGt                                                        >  1:209265/1‑149 (MQ=255)
                                                                                                gagggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGAt                                                    >  7:303804/1‑148 (MQ=255)
                                                                                                   ggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCAGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATc                                                >  1:77788/1‑149 (MQ=255)
                                                                                                       agggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAt                                            <  6:348732/149‑1 (MQ=255)
                                                                                                           cctcgTTGCGGCCCATTCACAGCGTAACGCCTCTCCAGCAGGTTCTGCGTTGTGCTGCAATCATGTAGCACAGTTGCAATGCGCGTGGGGTGGATAGTCGTGAAACAAATGGTTCGt                                                                        <  5:381292/114‑1 (MQ=255)
                                                                                                           ccGCGTTGCGGCCAATTCACAGCGTAACGCCGATCCCGCAGGGGCGGAGTTGTGCTGCAATCAGGGCGCACAGTTGCACAGCTCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                        >  6:381292/1‑117 (MQ=255)
                                                                                                            cgcgTTGCGGCCGATTCACAGCGTAAGGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGt                                         >  1:203225/1‑147 (MQ=255)
                                                                                                              cgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAgg                                     <  1:305758/149‑1 (MQ=255)
                                                                                                               gTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAgg                                     >  8:207463/1‑148 (MQ=255)
                                                                                                                ttGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGta                                   <  2:77788/149‑1 (MQ=255)
                                                                                                                   cGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTa                                        <  2:259462/141‑1 (MQ=255)
                                                                                                                   cGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTa                                        >  1:259462/1‑141 (MQ=255)
                                                                                                                       ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacctac                              >  6:81790/1‑147 (MQ=255)
                                                                                                                       ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCa                            <  2:209265/149‑1 (MQ=255)
                                                                                                                        cATTCACAGCGTAACGCCGCTCCCGCAGGGGCAGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacctac                              >  6:318150/1‑146 (MQ=255)
                                                                                                                        cATTAACAGCGTAACGCCGCTCCCGCAGGGGAGGAGTTGTGCTTCAATCAGGGCGCACAGTTGCCCTTCGCGTGGGGTGGATAGTAGTGCAACAAATGGTTCGTTAATCGCACAGATTGTCTATTATCGTATTGTAAGGTACCTACCAg                           >  6:283512/1‑149 (MQ=255)
                                                                                                                            cacaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtct                        <  4:13789/148‑1 (MQ=255)
                                                                                                                             acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                        <  8:306076/99‑1 (MQ=255)
                                                                                                                             acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                        >  7:306076/1‑99 (MQ=255)
                                                                                                                             acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTcc                      >  4:296535/1‑149 (MQ=255)
                                                                                                                                   tAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATc                <  4:367138/149‑1 (MQ=255)
                                                                                                                                   caaCGCCGCCCCTGCGGGAGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATc                >  3:277402/2‑149 (MQ=255)
                                                                                                                                    aaCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCt               >  7:12648/1‑149 (MQ=255)
                                                                                                                                     aCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCt               <  3:296535/148‑1 (MQ=255)
                                                                                                                                       gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAg                                      <  7:308690/123‑1 (MQ=255)
                                                                                                                                       gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAg                                      >  8:308690/1‑123 (MQ=255)
                                                                                                                                       gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTggg            >  7:226650/1‑149 (MQ=255)
                                                                                                                                       gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTggg            <  8:303804/149‑1 (MQ=255)
                                                                                                                                         cgcTCCCGCAGGGGCAGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGtt          <  5:318150/149‑1 (MQ=255)
                                                                                                                                            tCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGt                                    >  3:158231/1‑120 (MQ=255)
                                                                                                                                            tCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGt                                    <  4:158231/120‑1 (MQ=255)
                                                                                                                                            tCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTAGATTATCGGATTGTAACGTACCTACCAGTCTCCTTAATCAGGGATTAc       >  4:39965/1‑149 (MQ=255)
                                                                                                                                              ccGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTAc       <  1:286219/147‑1 (MQ=255)
                                                                                                                                                gCAGGGGCGGCGTTGTGCTGCAATCAGGGAGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCg     >  5:145444/1‑147 (MQ=255)
                                                                                                                                                 cAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAg  <  6:145444/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAG  >  NC_002947/1568766‑1569059

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