Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,713,570 A→T 5.8% L315Q (CTG→CAG)  paaK ← phenylacetate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,713,5700AT 5.8% 205.5 / 2.6 69L315Q (CTG→CAG) paaKphenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

TCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCC  >  NC_002947/3713424‑3713717
                                                                                                                                                    |                                                                                                                                                   
tctcATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGgcag                                                                                                                                                    >  7:324716/1‑148 (MQ=255)
    aTAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGgcagca                                                                                                                                                  <  3:314757/146‑1 (MQ=255)
    aTAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGcc                                                                                                                                               >  8:85698/1‑149 (MQ=255)
    aTAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGcc                                                                                                                                               >  6:3109/1‑149 (MQ=255)
    aTAAAGCTCAGAAAGCAGTTTTATTTATAGCACCTGTTCCTCGATCAGAGTCGGGAACACATTGACCCCGCGAATGATCAGCAAGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGcc                                                                                                                                               >  4:325405/1‑149 (MQ=255)
       aaGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCTGCTGCCgg                                                                                                                                             <  6:36968/148‑1 (MQ=255)
       aaGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCgg                                                                                                                                             <  7:30067/148‑1 (MQ=255)
               aaaGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTcg                                                                                                                                    >  4:176643/1‑149 (MQ=255)
                aaGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTcgc                                                                                                                                   <  8:63850/149‑1 (MQ=255)
                  gCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCgg                                                                                                                                 >  6:374725/1‑149 (MQ=255)
                      ttttATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTg                                                                                                                              <  3:331002/148‑1 (MQ=255)
                        ttATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCTGCTGCCgggc                                                                                                                                           >  8:114765/1‑132 (MQ=255)
                        ttATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCTGCTGCCgggc                                                                                                                                           <  7:114765/133‑2 (MQ=255)
                        ttATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCgg                                                                                                                           >  2:14288/1‑149 (MQ=255)
                           tttttAGCCACTGTGGGTCGAGCTGAGTCGGGAACACATTGACGCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGAATCGGCCGG‑GCGGTGCCCGGCAGCAGCGGGGTGAGGTCGCGGGTGCGGTAg                                                                                                                        >  6:18071/1‑149 (MQ=255)
                           tttttAGCACCTGTTCCTCGATCTGAGTCGGGCACCCATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAg                                                                                                                        >  7:158062/1‑149 (MQ=255)
                             tttAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAgcg                                                                                                                      >  6:54534/1‑149 (MQ=255)
                                   aacacttCCTCGATCTGAGTCGGGAACACATTGAACCCGCGAATGATCAGCATGTCGTCACTTCGGCCGGTGATCTTGCCGATCAGC‑CGCAAAGGCCGA‑GCTGTGCCCTGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCAt                                                                                                                >  2:358293/6‑149 (MQ=255)
                                          ccTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAt                                                                                                         >  6:307469/1‑149 (MQ=255)
                                          ccTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGAGCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCaa                                                                                                          <  4:296491/148‑1 (MQ=255)
                                             cGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGcc                                                                                                      >  1:247544/1‑149 (MQ=255)
                                                tCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGcctcc                                                                                                   >  6:259638/1‑149 (MQ=255)
                                                    aGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGc                                                                                               >  7:127471/1‑149 (MQ=255)
                                                      tCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTc                                                                                             <  7:253012/149‑1 (MQ=255)
                                                         ggAACACATTGACCGCGCTTATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCTGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGc                                                                                          <  2:67752/149‑1 (MQ=255)
                                                            acacATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAg                                                                                       <  6:87165/149‑1 (MQ=255)
                                                             cacaTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAgg                                                                                      <  4:350085/149‑1 (MQ=255)
                                                              acaTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGt                                                                                     >  1:97461/1‑149 (MQ=255)
                                                              acaTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGt                                                                                     <  8:127471/149‑1 (MQ=255)
                                                               caTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGt                                                                                     >  5:26382/1‑148 (MQ=255)
                                                               caTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTg                                                                                    <  5:307469/149‑1 (MQ=255)
                                                                aTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTg                                                                                    >  3:186804/1‑148 (MQ=255)
                                                                 ttGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGaa                                                                                  >  8:74520/1‑149 (MQ=255)
                                                                    aCCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATAc                                                                               >  1:18347/1‑149 (MQ=255)
                                                                     ccccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATAcc                                                                              <  6:26382/149‑1 (MQ=255)
                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTTAATACCa                                                                             <  2:247544/149‑1 (MQ=255)
                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATAc                                                                               >  7:210751/1‑147 (MQ=255)
                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATAc                                                                               >  2:240265/1‑147 (MQ=255)
                                                                       ccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCACCCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCa                                                                             <  1:311727/149‑1 (MQ=255)
                                                                       ccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCgg                                                                                                                                             >  8:236283/1‑84 (MQ=255)
                                                                       ccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCgg                                                                                                                                             <  7:236283/84‑1 (MQ=255)
                                                                        cgcgAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGc                                                                           >  8:75504/1‑149 (MQ=255)
                                                                            aaTGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCTGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGc                                                                       <  5:18071/149‑1 (MQ=255)
                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGaa                                                                                  >  8:146471/1‑133 (MQ=255)
                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGaa                                                                                  <  7:146471/133‑1 (MQ=255)
                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGt                                                                  >  5:14358/1‑149 (MQ=255)
                                                                                   aGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTg                                                                <  7:74520/149‑1 (MQ=255)
                                                                                     cATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGc                                                              <  5:54534/149‑1 (MQ=255)
                                                                                            tcACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGAGGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCg                                                        >  1:263487/1‑149 (MQ=255)
                                                                                              aCTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGgc                                                      <  7:85698/148‑1 (MQ=255)
                                                                                               cTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGgcag                                                    <  8:324716/149‑1 (MQ=255)
                                                                                                  cGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCAc                                                 <  5:374725/149‑1 (MQ=255)
                                                                                                        gTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCGTCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGcc                                           >  2:310564/1‑149 (MQ=255)
                                                                                                         tGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCg                                          <  5:259638/149‑1 (MQ=255)
                                                                                                           aTCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCgg                                         <  8:210751/148‑1 (MQ=255)
                                                                                                                 ccGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAgg                                  <  8:236109/149‑1 (MQ=255)
                                                                                                                 ccGATCCGA‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAg                                   <  8:158062/148‑1 (MQ=255)
                                                                                                                    aTCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTc                               >  4:54546/1‑149 (MQ=255)
                                                                                                                       cgc‑cgcATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgat                            <  1:14288/149‑1 (MQ=255)
                                                                                                                       cgc‑cgcATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgat                            <  2:263487/149‑1 (MQ=255)
                                                                                                                               tCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCg                     >  2:193728/1‑149 (MQ=255)
                                                                                                                                 ggCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCggg                   <  6:14358/149‑1 (MQ=255)
                                                                                                                                 gaccGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCgg                    <  4:319667/146‑1 (MQ=255)
                                                                                                                                  gCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGa                  <  5:382994/149‑1 (MQ=255)
                                                                                                                                    cGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGa                  >  8:188953/1‑147 (MQ=255)
                                                                                                                                         cctcacccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAg            >  6:62065/6‑149 (MQ=255)
                                                                                                                                           gTGGCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGGCTCCTTGCTCAGCGAGGTGAATTCCAGCTCGCCCAGTTGGCGATGCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTg          <  3:325405/149‑1 (MQ=255)
                                                                                                                                           gTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTg          <  2:18347/149‑1 (MQ=255)
                                                                                                                                           gTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTg          <  1:193728/149‑1 (MQ=255)
                                                                                                                                              cccGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTc       >  3:133781/1‑149 (MQ=255)
                                                                                                                                                 ggCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTc    <  3:176643/149‑1 (MQ=255)
                                                                                                                                                   cagcagCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAAGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTccc  <  4:324995/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGC‑CGCATCGGCCGA‑GCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTGAATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCCC  >  NC_002947/3713424‑3713717

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