Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,313,672 A→T 6.2% L139Q (CTG→CAG)  benE‑I ← benzoate transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,313,6720AT 6.2% 301.4 / 6.1 96L139Q (CTG→CAG) benE‑Ibenzoate transport protein
Reads supporting (aligned to +/- strand):  ref base A (51/39);  new base T (4/2);  total (55/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCCAGCAGCAGTATGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGACCAAGA  >  NC_002947/2313528‑2313819
                                                                                                                                                |                                                                                                                                                    
cccAGCAGCAGTATGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCcagcag                                                                                                                                                  <  5:23161/149‑1 (MQ=255)
 ccAGCAGCAGTATGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTAGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAgg                                                                                                                                                 >  8:235989/1‑149 (MQ=255)
        cagTATGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCcggcgg                                                                                                                                          <  3:317754/149‑1 (MQ=255)
             tGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCTGCTGGCcggcatcctgc                                                                                                                                     <  7:158274/149‑8 (MQ=255)
                cGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATgccg                                                                                                                                  >  1:139602/1‑149 (MQ=255)
                 gACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATgccgc                                                                                                                                 >  8:108827/1‑149 (MQ=255)
                    cACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATGAGCAACCCGGCGAGGCTGGTGTCGACCGACAGGCTGGCACAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCATTGCCGCCgg                                                                                                                              <  3:62398/149‑1 (MQ=255)
                      cGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGc                                                                                                                            <  7:108827/149‑1 (MQ=255)
                       ggCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGCCGCCAATGCCGCCGGGCa                                                                                                                           >  5:18645/1‑149 (MQ=255)
                            aaCGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTAcg                                                                                                                      >  6:98077/1‑149 (MQ=255)
                             aCGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTCcgc                                                                                                                     >  3:81201/1‑149 (MQ=255)
                                      gAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGcc                                                                                                            >  3:38692/1‑149 (MQ=255)
                                            tttGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAg                                                                                                              >  1:316659/1‑141 (MQ=255)
                                            tttGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAg                                                                                                              <  2:316659/141‑1 (MQ=255)
                                               gAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTAcgc                                                                                                                     >  6:6113/1‑131 (MQ=255)
                                               gAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTAcgc                                                                                                                     <  5:6113/131‑1 (MQ=255)
                                                aaaGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGaa                                                                                                     <  7:130886/146‑1 (MQ=255)
                                                aaaGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGaa                                                                                                     >  8:130886/1‑146 (MQ=255)
                                                aaaGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCa                                                                                                  >  2:167393/1‑149 (MQ=255)
                                                 aaGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCAc                                                                                                 >  6:7504/1‑149 (MQ=255)
                                                     cacaTACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGaa                                                                                                     <  6:329068/141‑1 (MQ=255)
                                                     cacaTACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGaa                                                                                                     >  5:329068/1‑141 (MQ=255)
                                                     cacaTACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCCAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATTCCGCCGGGCACTACGCGCACCACCCGCTCGAAGCACCCg                                                                                              >  5:210175/1‑148 (MQ=255)
                                                      acaTACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCt                                                                                            <  8:169024/149‑1 (MQ=255)
                                                          aCGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGc                                                                                        >  1:55750/1‑149 (MQ=255)
                                                           cGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGc                                                                                        <  5:7504/148‑1 (MQ=255)
                                                           cGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGc                                                                                        <  7:102967/148‑1 (MQ=255)
                                                            gCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCa                                                                                          <  6:184329/145‑1 (MQ=255)
                                                            gCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCa                                                                                          >  5:184329/1‑145 (MQ=255)
                                                             cggCAATCA‑CCACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCa                                                                                    <  8:181874/147‑1 (MQ=255)
                                                             cTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCa                                                                                                                           <  6:211080/111‑1 (MQ=255)
                                                             cTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCa                                                                                                                           >  5:211080/1‑111 (MQ=255)
                                                                 aaTCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGca                                                                                 >  5:135488/1‑149 (MQ=255)
                                                                  aTCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGcac                                                                                <  2:139602/149‑1 (MQ=255)
                                                                     aGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGcaccac                                                                             <  7:235989/149‑1 (MQ=255)
                                                                      gCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGcaccac                                                                             >  1:298389/1‑148 (MQ=255)
                                                                      gCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGcaccac                                                                             >  2:86028/1‑148 (MQ=255)
                                                                             ggCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCGGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCcg‑g                                                                    >  1:251605/1‑149 (MQ=255)
                                                                             ggCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCcg‑g                                                                    <  5:199155/149‑1 (MQ=255)
                                                                                gAGGCTGGGGTCGACAGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCcg‑gcgg                                                                 >  2:226494/1‑149 (MQ=255)
                                                                                   gCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGccc                                                               >  6:240221/1‑148 (MQ=255)
                                                                                             accgacAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTa                                                    <  6:135488/149‑1 (MQ=255)
                                                                                               cgacAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAgg                                                  <  5:98077/149‑1 (MQ=255)
                                                                                                gacAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAgg                                                  >  3:79185/1‑148 (MQ=255)
                                                                                                gacAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGc                                                 <  1:86028/149‑1 (MQ=255)
                                                                                                     gCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTa                                            >  3:209023/1‑149 (MQ=255)
                                                                                                      cTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTAt                                           <  1:167393/149‑1 (MQ=255)
                                                                                                      cTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTAt                                           >  4:100268/1‑149 (MQ=255)
                                                                                                       tGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACt                                                                                                                         >  8:78153/1‑71 (MQ=255)
                                                                                                       tGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACt                                                                                                                         <  7:78153/71‑1 (MQ=255)
                                                                                                       tGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGcc                                                               >  5:41161/1‑128 (MQ=255)
                                                                                                       tGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGcc                                                               <  6:41161/128‑1 (MQ=255)
                                                                                                       tGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTACGCACCTaag                                          >  1:309285/1‑147 (MQ=255)
                                                                                                        ggCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCcg‑gcg                                                                  >  6:154631/1‑124 (MQ=255)
                                                                                                        ggCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCcg‑gcg                                                                  <  5:154631/124‑1 (MQ=255)
                                                                                                         gctaaGGCACCGATGCCGAACTGCAGGAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCACTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGc                                        <  6:18645/146‑1 (MQ=255)
                                                                                                         gCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGc                                        >  5:148447/1‑149 (MQ=255)
                                                                                                            aaGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTc                                     <  4:38692/149‑1 (MQ=255)
                                                                                                             aGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCg                                    <  2:298389/149‑1 (MQ=255)
                                                                                                               gCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGaa                                                                          >  4:86213/1‑110 (MQ=255)
                                                                                                               gCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGaa                                                                          <  3:86213/110‑1 (MQ=255)
                                                                                                                  ccGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTa                               <  4:79185/149‑1 (MQ=255)
                                                                                                                   cGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAg                              >  5:290059/1‑149 (MQ=255)
                                                                                                                   cGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAg                              >  7:45600/1‑149 (MQ=255)
                                                                                                                    gATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTa                                                    <  2:75907/126‑1 (MQ=255)
                                                                                                                    gATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTa                                                    >  1:75907/1‑126 (MQ=255)
                                                                                                                    gATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAgg                             <  2:309285/149‑1 (MQ=255)
                                                                                                                     aTGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAggg                            >  7:51120/1‑149 (MQ=255)
                                                                                                                            aCTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTg                     <  6:290059/149‑1 (MQ=255)
                                                                                                                              tGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTg                     <  8:225386/147‑1 (MQ=255)
                                                                                                                              tGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTg                     >  7:225386/1‑147 (MQ=255)
                                                                                                                              tGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGc                   >  3:13399/1‑149 (MQ=255)
                                                                                                                                cagcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCGCGCGGCCGACAGCAGGTAGGCACCTATGGCCt                                      <  1:34909/129‑1 (MQ=255)
                                                                                                                                cagcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCt                                      >  2:34909/1‑128 (MQ=255)
                                                                                                                                cagcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGc                   >  2:111521/1‑147 (MQ=255)
                                                                                                                                 agcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTAt                                           <  8:74781/122‑1 (MQ=255)
                                                                                                                                 agcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTAt                                           >  7:74781/1‑122 (MQ=255)
                                                                                                                                 agcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAg                >  5:210857/1‑149 (MQ=255)
                                                                                                                                  gcagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAgg               >  7:103368/1‑149 (MQ=255)
                                                                                                                                   cagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAgg                             >  8:286137/1‑134 (MQ=255)
                                                                                                                                   cagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAgg                             <  7:286137/134‑1 (MQ=255)
                                                                                                                                   cagGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAggg              >  2:196730/1‑149 (MQ=255)
                                                                                                                                    agGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCatct              >  1:178488/1‑145 (MQ=255)
                                                                                                                                    agGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGc             <  3:100268/149‑1 (MQ=255)
                                                                                                                                     tgccgCCGGCCTGCTGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCg            >  5:336434/5‑149 (MQ=255)
                                                                                                                                     ggATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCg            <  2:55750/149‑1 (MQ=255)
                                                                                                                                      gATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCgg           >  3:81028/1‑149 (MQ=255)
                                                                                                                                      gATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCgg           >  1:80087/1‑149 (MQ=255)
                                                                                                                                       aTGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGt          <  8:45600/149‑1 (MQ=255)
                                                                                                                                        cgCCGGCCTGCTGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAgg                                                      <  7:284358/103‑1 (MQ=255)
                                                                                                                                        cgCCGGCCTGCTGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAgg                                                      >  8:284358/2‑104 (MQ=255)
                                                                                                                                        cgCCGGCCTGCTGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGt          >  1:220430/2‑148 (MQ=255)
                                                                                                                                        cgCCGGCCTGCTGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGt          >  1:215833/2‑148 (MQ=255)
                                                                                                                                         gCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGa        >  5:264232/1‑149 (MQ=255)
                                                                                                                                           cGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGAc       >  6:170025/1‑148 (MQ=255)
                                                                                                                                            ggCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGAcc      >  4:48497/1‑148 (MQ=255)
                                                                                                                                             gCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGACCaa    <  5:240221/149‑1 (MQ=255)
                                                                                                                                               cagcagGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGACCAAGa  <  6:264232/149‑1 (MQ=255)
                                                                                                                                                agcagGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAggg              >  8:241548/1‑136 (MQ=255)
                                                                                                                                                agcagGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAggg              <  7:241548/136‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCCAGCAGCAGTATGGCGACCACGGCGTAACGGGCGGTGAAGCGTTTGAAAGCCACATACGCTGCAATCAGCAACCCGGCGAGGCTGGGGTCGACCGACAGGCTGGCAAAGGCACCGATGCCGAACTGCAGCAGGATGCCGGCCAGCAGGCCGGCGGCAATGCCGCCGGGCACTACGCGCACCAGCCGCTCGAAGCACCCGCTCAGCCCCAGCACCACGAAGCCG‑GCGGCCGACAGCAGGTAGGCACCTATGGCCTCGCTGTAGGGCACGCTGGCCAGGGCGGTGACCAAGA  >  NC_002947/2313528‑2313819

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