Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,137,005 A→T 7.5% intergenic (‑255/‑136) fadE ← / → gstB medium‑long chain acyl‑CoA dehydrogenase/glutathione S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,137,0050AT 7.5% 193.8 / 5.3 67intergenic (‑255/‑136)fadE/gstBmedium‑long chain acyl‑CoA dehydrogenase/glutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (2/3);  total (31/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 = 4.58e-01

AGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGAAGA  >  NC_002947/2136858‑2137150
                                                                                                                                                    |                                                                                                                                                 
aGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTa                                                                                                                                                   >  2:186788/1‑148 (MQ=255)
   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAGAAc                                                                                                                                               >  7:70970/1‑149 (MQ=255)
   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                                                                                               >  7:76111/1‑149 (MQ=255)
   tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                                                                                               >  2:72643/1‑149 (MQ=255)
     aTTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                                                                                               <  6:193140/147‑1 (MQ=255)
         gTTCGCCTGTACCGGCCTATTCGCGGGTAACCCCCTCCTACAACGGAACGTGGCGATGCCCCCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACg                                                                                                                                          >  1:352376/1‑148 (MQ=255)
          ttCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTTTTACCCAcac                                                                                                                                         <  8:179100/148‑3 (MQ=255)
             gCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACTACGGTACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAgg                                                                                                                                      <  6:60010/148‑1 (MQ=255)
             gCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAgg                                                                                                                                      >  6:49099/1‑148 (MQ=255)
                    ccGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCg                                                                                                                               <  5:134139/148‑1 (MQ=255)
                    ccGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCg                                                                                                                               <  3:359432/148‑1 (MQ=255)
                       gCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGa                                                                                                                           >  7:309814/1‑149 (MQ=255)
                          tATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGt                                                                                                                        >  6:156052/1‑149 (MQ=255)
                                 gggTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACg                                                                                                                 <  6:130881/149‑1 (MQ=255)
                                    tAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACggtg                                                                                                              <  1:186788/149‑1 (MQ=255)
                                      aaCCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACggtggt                                                                                                            <  5:156052/149‑1 (MQ=255)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTTTTACCcacacata                                                                                                                                      <  6:321623/120‑3 (MQ=255)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTTTTACCcacacata                                                                                                                                      >  5:321623/1‑118 (MQ=255)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTTTTACCcacacat                                                                                                                                       <  2:258084/119‑2 (MQ=39)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTTTTACCcacacat                                                                                                                                       >  1:258084/1‑118 (MQ=255)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATAc                                                                                                                                 >  7:67374/1‑125 (MQ=255)
                                         cccTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATAc                                                                                                                                 <  8:67374/125‑1 (MQ=255)
                                         cccTCCTACAACGGAAAGTGGCGATGCCAACGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAg                                                                                                         >  3:320462/1‑149 (MQ=255)
                                              ctACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgt                                                                                                    <  8:76111/149‑1 (MQ=255)
                                              ctACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgt                                                                                                    <  8:70970/149‑1 (MQ=255)
                                                aCAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTgg                                                                                                  >  6:349851/1‑149 (MQ=255)
                                                  aaCGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTgg                                                                                                  >  7:111038/1‑147 (MQ=255)
                                                  aaCGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTgg                                                                                                  <  8:304687/147‑1 (MQ=255)
                                                  aaCGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTgg                                                                                                  >  7:304687/1‑147 (MQ=255)
                                                  aaCGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGAGTGGGTAAAACCCACGCAGGACTACCGGTTACGAGAAGACGGTGGTGAGCCTGTGGGa                                                                                                >  3:45561/1‑149 (MQ=255)
                                                    cGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGc                                                                                              <  7:38626/149‑1 (MQ=255)
                                                     ggAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgt                                                                                                    >  8:351557/1‑142 (MQ=255)
                                                     ggAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgt                                                                                                    <  7:351557/142‑1 (MQ=255)
                                                       accGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTATTACCCACGGAGGAATACCGGTGACGTGAAGTCGGTGGTGAGCCTGTGGGAGCggg                                                                                           <  2:298124/147‑1 (MQ=255)
                                                          gTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACTGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCtt                                                                                        >  5:306762/1‑149 (MQ=255)
                                                                aTGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACTGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCcgc                                                                                  <  6:306762/149‑1 (MQ=255)
                                                                 cgCCACCGCATTCATCAATCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCcgc                                                                                  >  2:196172/2‑148 (MQ=255)
                                                                   ccACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCcgcg                                                                                 >  1:327678/1‑147 (MQ=255)
                                                                    cACCGATTGATGAATGCGGTGGCGGCCGCGTCGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCcgcg                                                                                 >  1:31240/1‑146 (MQ=255)
                                                                     aCCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGACCGGGCTTGCCCGCGaa                                                                               >  7:271197/1‑147 (MQ=255)
                                                                      ccGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                            <  4:338109/149‑1 (MQ=255)
                                                                       cGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGc                                                                                          >  3:323194/1‑134 (MQ=255)
                                                                        ggcttGATGTATGGGGTGGCGGCGGCGACGGTGAATGACCGGGGGGTTGGCTGGGAGGGTCGGCCTATGTGTGGGTAAATCCGACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGc                                                                                          <  4:323194/131‑1 (MQ=255)
                                                                        gATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                          <  2:327678/149‑1 (MQ=255)
                                                                          aTTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACGCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCACGAAGAGGcc                                                                         >  6:198053/1‑149 (MQ=255)
                                                                          aTTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGcc                                                                         >  3:200463/1‑149 (MQ=255)
                                                                                  aTGCTGTGGCGGCAGCTACGGTTAATGACCGGGGGGTTGTCTTGGAGGCTCGGCCTATGTGTTTGTAAACCCCACGAAGGAATACCGGTGAGGTGAAGAGGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAg                                                                  <  4:45561/148‑1 (MQ=37)
                                                                                          gcggcCGCGAGGGTGAATGACCGGGGGGGTGGGTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGTGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGaa                                                          <  8:309814/148‑1 (MQ=255)
                                                                                          gcggcCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAg                                                         <  1:72643/149‑1 (MQ=255)
                                                                                               cgcgACGGTTAATGCCCGGGGGGTTGGGTGGGATGCTCGGCCTTTGTTGTTGTAAAACCCACGGATGAATACCGTTGACGTGTAGTCGGTGGGGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTg                                                    <  8:271197/149‑1 (MQ=25)
                                                                                               cgcgACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTg                                                    <  5:49099/149‑1 (MQ=255)
                                                                                                 cgACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACTCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCACGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGa                                                  <  5:198053/149‑1 (MQ=255)
                                                                                                       tGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGAc                                            >  3:140776/1‑149 (MQ=255)
                                                                                                            gACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGc                                       >  6:5155/1‑149 (MQ=255)
                                                                                                              ccGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGcc                                      <  2:31240/148‑1 (MQ=255)
                                                                                                               cccgggcTTGGCTGGCAGTCTCGGACTATGTGTGGGTAAAACACACGCAGAAATACCTGTCACGAGAAGACGCTGGTGAGCATGTCGGAGCGGGCTTGCCCGCGAAGAGGCACGAACATGCCGCGAAGCACTGGAGTAGACTTAGCCAc                                    >  8:179114/4‑149 (MQ=255)
                                                                                                               cGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCAc                                    <  1:196172/149‑1 (MQ=255)
                                                                                                                      ttGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCa                             <  4:200463/149‑1 (MQ=255)
                                                                                                                         gctggcAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGa                                                           <  7:53502/116‑1 (MQ=255)
                                                                                                                         gctggcAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGa                                                           >  8:53502/1‑116 (MQ=255)
                                                                                                                          ctggcAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTc                          >  5:184198/1‑148 (MQ=255)
                                                                                                                             gcAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCaaaa                      <  5:5155/149‑1 (MQ=255)
                                                                                                                                 gCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCGACCTCCGCTCTCAAAAGGAt                  >  3:325636/1‑149 (MQ=255)
                                                                                                                                 gCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGAt                  <  4:169869/149‑1 (MQ=255)
                                                                                                                                       ccTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGcc            <  4:140776/149‑1 (MQ=255)
                                                                                                                                               tgGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGaa    <  5:349851/149‑1 (MQ=255)
                                                                                                                                                 ggTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTg      >  4:261998/1‑145 (MQ=255)
                                                                                                                                                 ggTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGGGAGCCTGTGGGAGCTGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAATGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGAAGa  <  1:195714/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGCTCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGAAGA  >  NC_002947/2136858‑2137150

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: