Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,628 A→G 8.0% I255V (ATC→GTC) ‡ benR → BenABC operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,580,6280AG 8.0% 174.4 / 5.1 62I255V (ATC→GTC) ‡benRBenABC operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (34/23);  new base G (2/3);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCG  >  NC_002947/3580491‑3580773
                                                                                                                                         |                                                                                                                                                 
ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                        <  1:234311/149‑1 (MQ=255)
ccccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                        >  8:23839/1‑149 (MQ=255)
 cccGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAg                                                                                                                                        >  6:277431/1‑148 (MQ=255)
       ggCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGCCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCGTGATGCCCAAgtgctggt                                                                                                                                 <  2:112214/149‑8 (MQ=255)
          tGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgc                                                                                                                              >  7:273573/1‑149 (MQ=255)
            ggCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTgcgc                                                                                                                            <  3:99104/149‑1 (MQ=255)
                        ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCAAgtgctg                                                                                                                                   >  4:141386/1‑125 (MQ=255)
                        ttCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCGCTTGGGCGTGATGCCCAAgtgctg                                                                                                                                   <  3:141386/130‑6 (MQ=255)
                         tCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                >  4:308926/1‑148 (MQ=255)
                             gATTACATCGATAGCTACCTGATATTGGAATTGTCTGCTGAAGTGGTTGCAGTGCAGGCTTGCATGAGTTTGGGTTGGCTGTATGCGTTGTTCGACCTGGACTTTGTCATCAGGCCCATGCTTTACTTGCGCCAGCGGAAACTGTAAc                                                                                                            <  4:174923/148‑1 (MQ=255)
                             gATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAAc                                                                                                            <  8:273573/148‑1 (MQ=255)
                              aTTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACgg                                                                                                          <  7:23839/149‑1 (MQ=255)
                                tACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                   >  7:207150/1‑138 (MQ=255)
                                tACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGaaa                                                                                                                   <  8:207150/138‑1 (MQ=255)
                                   aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCa                                                                                                     >  2:3407/1‑149 (MQ=255)
                                   aTCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAACAGTACGTGCGCCAGCGGAAACTGGACCGGGTGCa                                                                                                     >  8:249297/1‑149 (MQ=255)
                                      gAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgc                                                                                                  >  6:64901/1‑149 (MQ=255)
                                        aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                <  3:10214/133‑1 (MQ=255)
                                        aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTg                                                                                                                >  4:10214/1‑133 (MQ=255)
                                        aCGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcct                                                                                                >  1:243297/1‑149 (MQ=255)
                                         cGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctg                                                                                               >  1:68303/1‑149 (MQ=255)
                                              ccTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCAAGAGTTTGCGTTCGCTGTATGCGATGTTCGACCAGCACTTGGGCAACACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGCTGAATgcctgcct                                                                                            >  8:30601/1‑147 (MQ=255)
                                               cTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAg                                                                                         >  2:39025/1‑149 (MQ=255)
                                                         aaTTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGTGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACt                                                                               >  6:59592/1‑149 (MQ=255)
                                                           ttGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGAGCATGCCTGCCTGAGCGACCCTACTtg                                                                             >  4:150021/1‑149 (MQ=255)
                                                                ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                         <  5:277431/148‑1 (MQ=255)
                                                                ggCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGc                                                                         <  3:308926/148‑1 (MQ=255)
                                                                 gCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGt                                                                       <  2:243297/149‑1 (MQ=255)
                                                                  cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTg                                                                      <  2:221451/149‑1 (MQ=255)
                                                                  cTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTg                                                                      >  3:257305/1‑149 (MQ=255)
                                                                          cTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAgtgt                                                              >  8:245058/1‑149 (MQ=255)
                                                                             gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                           <  7:245058/149‑1 (MQ=255)
                                                                             gCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGAc                                                           <  5:64901/149‑1 (MQ=255)
                                                                                   cAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAAtt                                                     >  6:103652/1‑149 (MQ=255)
                                                                                    aGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCTGAAACTGGAACGGGTGCATGCCTGCATGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTg                                                    >  4:301969/1‑149 (MQ=255)
                                                                                      gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                            >  5:27530/1‑107 (MQ=255)
                                                                                      gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATgcctgcct                                                                                            <  6:27530/107‑1 (MQ=255)
                                                                                      gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                       <  1:9642/144‑1 (MQ=255)
                                                                                      gCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGaa                                                       >  2:9642/1‑144 (MQ=255)
                                                                                        ttGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGccc                                                >  8:23362/1‑149 (MQ=255)
                                                                                           cATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTg                                              >  2:38328/1‑148 (MQ=255)
                                                                                            aTGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                            <  1:39025/149‑1 (MQ=255)
                                                                                            aTGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGa                                            >  8:41363/1‑149 (MQ=255)
                                                                                             tgagtttgagttCGCTGTATGCGTTGTTAGACCAGCACTTGGGCATCACGCACAAGCATTACGTGCTCCAGCGGAAACTGGAACGTGTGCATGCCTGCCTGAGCGACCCTACTTGTGGAGTGCGCAGTGTGACGGAATTGGCCCTGGAt                                           >  6:218497/1‑149 (MQ=255)
                                                                                             tGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGAt                                           >  5:134285/1‑149 (MQ=255)
                                                                                                   tGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGACAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGt                                     >  1:250026/1‑149 (MQ=255)
                                                                                                       ttCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                              >  7:104217/1‑88 (MQ=255)
                                                                                                       ttCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCAtgcctgc                                                                                              <  8:104217/88‑1 (MQ=255)
                                                                                                          gCTGTATGTGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTc                              <  5:59592/149‑1 (MQ=255)
                                                                                                            tGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGAGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCAt                            >  6:145104/1‑149 (MQ=255)
                                                                                                              tATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCt                          <  2:68303/149‑1 (MQ=255)
                                                                                                                tGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTgg                        >  3:98999/1‑149 (MQ=255)
                                                                                                                  cGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGAGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTggg                       <  3:150021/148‑1 (MQ=255)
                                                                                                                   gttgttCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGc                     >  2:9032/1‑149 (MQ=255)
                                                                                                                          gACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTc              >  1:316680/1‑149 (MQ=255)
                                                                                                                            ccAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc             >  2:15269/1‑148 (MQ=255)
                                                                                                                             cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGACCGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTcc             >  3:289524/1‑147 (MQ=255)
                                                                                                                             cAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCACTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCga           <  7:187512/149‑1 (MQ=255)
                                                                                                                               gCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga         <  7:23362/149‑1 (MQ=255)
                                                                                                                                tgcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCGGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga         <  2:195520/146‑1 (MQ=255)
                                                                                                                                tgcTTGGGCGTGATGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCGGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCgaga         >  1:195520/3‑148 (MQ=255)
                                                                                                                                      ggCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCg  <  3:191820/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
CCCCGCAGGCTGGGCTAGAGCGCATTCTTGATTACATCGAACGCAACCTGAAATTGGAATTGTCGGCTGAAGTGCTGGCAGAGCAGGCTTGCATGAGTTTGCGTTCGCTGTATGCGTTGTTCGACCAGCACTTGGGCATCACGCCCAAGCATTACGTGCGCCAGCGGAAACTGGAACGGGTGCATGCCTGCCTGAGCGACCCTACTTGTGGCGTGCGCAGTGTGACGGAATTGGCCCTGGATTACGGTTTCCTTCATCTGGGGCGGTTTTCCGAGATTTATCG  >  NC_002947/3580491‑3580773

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