Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,137,007 A→T 8.5% intergenic (‑257/‑134) 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,0070AT 8.5% 190.0 / 5.8 67intergenic (‑257/‑134)fadE/gstBmedium‑long chain acyl‑CoA dehydrogenase/glutathione S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGAAGA  >  NC_002947/2136861‑2137150
                                                                                                                                                   |                                                                                                                                               
tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAGAAc                                                                                                                                               >  7:70970/1‑149 (MQ=255)
tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                                                                                               >  2:72643/1‑149 (MQ=255)
tCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAAc                                                                                                                                               >  7:76111/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:70970/149‑1 (MQ=255)
                                           ctACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCtgt                                                                                                    <  8:76111/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)
                                                                                                                                                   |                                                                                                                                               
TCATTTGTTCGCCTGTACCGGCCTATTCGCGGGTAAACCCCTCCTACAACGGAACGTGGCGATGCCACCGATTGATGAATGCGGTGGCGGCCGCGACGGTGAATGACCGGGGGGTTGGCTGGCAGGCTCGGCCTATGTGTGGGTAAAACCCACGCAGGAATACCGGTGACGTGAAGACGGTGGTGAGCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCGGCGAAGGTCTGGAGTAGACTGTGCCACCTCCGCACTCAAAAGGATGTTGCCATGCTGAAGA  >  NC_002947/2136861‑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: