Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,630 C→G 6.2% I255M (ATC→ATG)  benR → BenABC operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,580,6300CG 6.2% 208.4 / 3.6 65I255M (ATC→ATG) benRBenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base G (2/2);  total (35/30)
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.30e-01

TGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3580485‑3580777
                                                                                                                                                  |                                                                                                                                                   
tGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACg                                                                                                                                                  >  8:173521/1‑149 (MQ=255)
  aGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGc                                                                                                                                                 <  3:100341/148‑1 (MQ=255)
      ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAg                                                                                                                                            >  7:339944/1‑149 (MQ=255)
       cccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                            >  3:138555/1‑148 (MQ=255)
        ccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATc                                                                                                                                                     >  4:207862/1‑138 (MQ=255)
            aGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgctgg                                                                                                                                      <  8:339944/149‑7 (MQ=255)
              gCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGt                                                                                                                                    <  5:217920/149‑1 (MQ=255)
                  gggtAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgcgc                                                                                                                                >  3:39969/4‑149 (MQ=255)
                     tAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAg                                                                                                                             >  8:112684/1‑149 (MQ=255)
                         gcgcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCgg                                                                                                                          >  3:67421/1‑148 (MQ=255)
                          cgcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGa                                                                                                                         >  8:42983/1‑148 (MQ=255)
                             aTTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                       >  8:266595/1‑147 (MQ=255)
                             aTTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                       <  7:266595/147‑1 (MQ=255)
                                 ttGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                                 >  4:316622/1‑149 (MQ=255)
                                    aTTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                              <  6:176808/149‑1 (MQ=255)
                                       aCATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGt                                                                                                            <  7:42983/148‑1 (MQ=255)
                                        cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                        <  6:83709/119‑5 (MQ=255)
                                        cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                        >  5:83709/1‑115 (MQ=255)
                                        cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTg                                                                                                           >  5:143538/1‑148 (MQ=255)
                                            gAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgc                                                                                                      <  6:143538/149‑1 (MQ=255)
                                              aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcct                                                                                                    >  3:78568/1‑149 (MQ=255)
                                               cGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctg                                                                                                   >  3:197169/1‑149 (MQ=255)
                                                gCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACt                                                                                                                     >  6:57777/1‑130 (MQ=255)
                                                gCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACt                                                                                                                     <  5:57777/130‑1 (MQ=255)
                                                  aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCGTTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                                >  2:89846/1‑149 (MQ=255)
                                                  aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                                >  7:113610/1‑149 (MQ=255)
                                                  aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                                <  4:67421/149‑1 (MQ=255)
                                                     cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                             >  6:21756/1‑149 (MQ=255)
                                                     cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                             <  5:21756/149‑1 (MQ=255)
                                                      tGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCaa                                                                                                                                             <  6:108332/100‑1 (MQ=255)
                                                      tGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCaa                                                                                                                                             >  5:108332/1‑100 (MQ=255)
                                                        aaaTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCGAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGa                                                                                          >  2:304789/1‑149 (MQ=255)
                                                         aaTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAc                                                                                         <  4:78568/149‑1 (MQ=255)
                                                          aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                        <  3:328774/149‑1 (MQ=255)
                                                          aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                        <  6:11998/149‑1 (MQ=255)
                                                           ttGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAccc                                                                                       <  7:173521/149‑1 (MQ=255)
                                                                 ttGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGATCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACtt                                                                                  <  4:138555/148‑1 (MQ=255)
                                                                    tCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTtgtg                                                                               <  8:277258/148‑1 (MQ=255)
                                                                           aaGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTgc                                                                         <  4:197169/147‑1 (MQ=255)
                                                                                cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgtgt                                                                  <  3:316622/149‑1 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTCGGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                             >  1:34480/1‑136 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTg                                                                                                           <  1:60335/105‑1 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTg                                                                                                           >  2:60335/1‑105 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                             <  2:34480/135‑1 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                               >  1:358527/1‑149 (MQ=255)
                                                                                   gCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                               >  8:289538/1‑149 (MQ=255)
                                                                                         cAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAtt                                                         >  5:94656/1‑149 (MQ=255)
                                                                                            gCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGc                                                      >  8:2303/1‑149 (MQ=255)
                                                                                              ttGCATGAGTTTGCGTTC‑GCGGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGccc                                                    >  8:135343/1‑149 (MQ=255)
                                                                                                  aTGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                                >  3:178204/1‑149 (MQ=255)
                                                                                                                     tATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCt                              >  5:349255/1‑149 (MQ=255)
                                                                                                                          gttgttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGGCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGc                         >  1:269072/1‑149 (MQ=255)
                                                                                                                          gttgttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGc                         >  6:336652/1‑149 (MQ=255)
                                                                                                                             gttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCgg                       <  8:113610/148‑1 (MQ=255)
                                                                                                                                  acgaGCACTTGGGCATCACGGCCAATCATTTCGTGCGCCATGGGAAACTGGAACGGGTGCATACCTGCGTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGATTTGGCCCTTGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                 <  5:138828/146‑1 (MQ=255)
                                                                                                                                  aCCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                 >  4:129598/1‑149 (MQ=255)
                                                                                                                                   ccAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGGACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                                >  7:322208/1‑85 (MQ=255)
                                                                                                                                   ccAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGGACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                                <  8:322208/85‑1 (MQ=255)
                                                                                                                                    cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCga               <  7:112684/149‑1 (MQ=255)
                                                                                                                                        aCTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTTGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGAtt           >  2:341531/1‑149 (MQ=255)
                                                                                                                                        aCTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGAtt           >  3:372861/1‑149 (MQ=255)
                                                                                                                                          ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTa         >  3:220494/1‑149 (MQ=255)
                                                                                                                                             ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTGCGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATc       <  5:341919/148‑1 (MQ=255)
                                                                                                                                             ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCg      <  7:29354/149‑1 (MQ=255)
                                                                                                                                              gCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCgg     <  3:129598/149‑1 (MQ=255)
                                                                                                                                               cATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgc    <  2:269072/149‑1 (MQ=255)
                                                                                                                                                  cACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgcag  >  8:300977/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTC‑GCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3580485‑3580777

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