Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,630 C→G 7.6% 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 7.6% 199.7 / 4.8 63I255M (ATC→ATG) ‡benRBenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base G (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3580482‑3580777
                                                                                                                                                      |                                                                                                                                                   
gCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATc                                                                                                                                                     >  8:1050/1‑149 (MQ=255)
   tGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACg                                                                                                                                                  >  1:286465/1‑149 (MQ=255)
               aGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGCCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTAc                                                                                                                                      <  1:107595/149‑1 (MQ=255)
                   tGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                                  <  3:57956/149‑1 (MQ=255)
                        tAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCTGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAg                                                                                                                             >  6:268690/1‑149 (MQ=255)
                         agagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                            >  4:242432/1‑149 (MQ=255)
                          gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCTGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                            <  5:268690/148‑1 (MQ=255)
                          gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTTCCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                            <  1:236409/149‑1 (MQ=255)
                          gagCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGc                                                                                                                            <  6:52134/148‑1 (MQ=255)
                            gcgcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGa                                                                                                                         <  4:91684/149‑1 (MQ=255)
                              gcATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                       <  5:126187/149‑1 (MQ=255)
                                   cTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTgg                                                                                                                   <  5:238054/148‑1 (MQ=255)
                                    ttGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGaa                                                                                                                 >  8:257368/1‑149 (MQ=255)
                                        ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                              <  8:223843/148‑1 (MQ=255)
                                        ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCAGAAACTGGAACggg                                                                                                             >  6:73890/1‑149 (MQ=255)
                                        ttACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGATTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACggg                                                                                                             >  3:234906/1‑149 (MQ=255)
                                           cATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGc                                                                                                          <  2:286465/149‑1 (MQ=255)
                                            aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                         <  3:242432/149‑1 (MQ=255)
                                              cGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAt                                                                                                        <  2:35600/148‑1 (MQ=255)
                                              cGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAt                                                                                                        >  1:35600/1‑148 (MQ=255)
                                                 aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcc                                                                                                     <  1:119541/148‑1 (MQ=255)
                                                   gCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATTGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctg                                                                                                   >  1:275652/1‑149 (MQ=255)
                                                    cAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgcc                                                                                                 >  4:141315/1‑149 (MQ=255)
                                                     aaCCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                            >  6:315915/1‑105 (MQ=255)
                                                       ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGa                                                                                              <  7:74499/149‑1 (MQ=255)
                                                        cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                    <  6:116819/126‑1 (MQ=255)
                                                        cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                    >  5:116819/1‑126 (MQ=255)
                                                             aTTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGAcc                                                                                        <  7:1050/149‑1 (MQ=255)
                                                                   aTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                        >  8:17956/1‑91 (MQ=255)
                                                                   aTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCGTGATGCCCAAgtgct                                                                                                                                        <  7:17956/95‑5 (MQ=255)
                                                                      gTCGGCTGAAGTGCTTGCATTGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGTCCAGCACTTGGGCATGACGCCCAAGCTTTACGTGCGCCAGCGGAAACTGGTACGGGTGCATGCCTGCCTGAGCGACCCTACTtgt                                                                                <  3:130858/148‑1 (MQ=255)
                                                                       tCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                                  <  2:151909/129‑1 (MQ=255)
                                                                       tCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                                  >  1:151909/1‑129 (MQ=255)
                                                                          gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                           >  6:306812/1‑149 (MQ=255)
                                                                          gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                           >  4:320112/1‑149 (MQ=255)
                                                                                     ggCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGa                                                                >  7:234841/1‑149 (MQ=255)
                                                                                        agagCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACg                                                              >  2:115984/1‑148 (MQ=255)
                                                                                          agCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                           >  5:274567/1‑149 (MQ=255)
                                                                                           gCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAt                                                          <  3:320112/149‑1 (MQ=255)
                                                                                                    cATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCAGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTgg                                                 <  5:73890/149‑1 (MQ=255)
                                                                                                          tttGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGTGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACg                                           >  6:117563/1‑149 (MQ=255)
                                                                                                          tttGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGAGCAGTGTGACGGAATTGGCCCTGGATTACg                                           >  8:161554/1‑149 (MQ=255)
                                                                                                                ttCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACTGTTTcc                                     <  4:234906/149‑1 (MQ=255)
                                                                                                                  cGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCt                                    <  3:141315/148‑1 (MQ=255)
                                                                                                                          tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                            >  5:328644/1‑149 (MQ=255)
                                                                                                                          tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                            >  2:27771/1‑149 (MQ=255)
                                                                                                                          tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                            >  3:56128/1‑149 (MQ=255)
                                                                                                                          tGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                            >  1:154816/1‑149 (MQ=255)
                                                                                                                             gttgtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTc                                      >  1:181635/1‑136 (MQ=255)
                                                                                                                             gttgtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTc                                      <  2:181635/136‑1 (MQ=255)
                                                                                                                                gtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGt                      >  1:332380/1‑149 (MQ=255)
                                                                                                                                gtt‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGt                      >  5:112808/1‑149 (MQ=255)
                                                                                                                                     gACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGtttt                   <  7:257368/148‑1 (MQ=255)
                                                                                                                                      aCCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                 >  7:264837/1‑149 (MQ=255)
                                                                                                                                      aCCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc                 >  6:293812/1‑149 (MQ=255)
                                                                                                                                        cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCga               <  6:274567/149‑1 (MQ=255)
                                                                                                                                          gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              <  7:302741/148‑1 (MQ=255)
                                                                                                                                          gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              >  8:302741/1‑148 (MQ=255)
                                                                                                                                            gcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTa                                             <  1:295133/114‑1 (MQ=255)
                                                                                                                                            gcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTa                                             >  2:295133/2‑115 (MQ=255)
                                                                                                                                              ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              <  6:327306/144‑1 (MQ=255)
                                                                                                                                              ttGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgag              >  5:327306/1‑144 (MQ=255)
                                                                                                                                               tGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTACTTCATCTGGGGCGGTTTTCCGAGATTTAt        >  6:330946/1‑149 (MQ=255)
                                                                                                                                                     tCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGgcag  <  6:112808/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCTGAGTGCCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTA‑TGCGTTGTT‑CGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCGGCAG  >  NC_002947/3580482‑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: ‑