Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,249,903 T→C 6.9% intergenic (‑82/+95) PP_4624 ← / ← cidB alpha/beta family hydrolase/holin CidA anti‑holin regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,249,9030TC 6.9% 172.2 / 2.2 58intergenic (‑82/+95)PP_4624/cidBalpha/beta family hydrolase/holin CidA anti‑holin regulator
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base C (2/2);  total (27/31)
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.97e-01

CCCGGAACATGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTAAGCAGCCCGGTCAG  >  NC_002947/5249758‑5250051
                                                                                                                                                 |                                                                                                                                                    
cccGGAACATGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGctctc                                                                                                                                                   >  7:353362/1‑149 (MQ=255)
 ccGGAACATGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGctctc                                                                                                                                                   >  1:155172/1‑148 (MQ=255)
     aaCATGAAACTGTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTATTGGTATTGGCGGAATATGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGt                                                                                                                                              <  7:264267/149‑1 (MQ=255)
        aTGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGa                                                                                                                                           <  1:160561/149‑1 (MQ=255)
         tGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAg                                                                                                                                          <  7:296860/149‑1 (MQ=255)
         tGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAg                                                                                                                                          <  8:353362/149‑1 (MQ=255)
         tGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAg                                                                                                                                          <  5:159415/149‑1 (MQ=255)
          gAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAg                                                                                                                                          >  6:261026/1‑148 (MQ=255)
          gAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACATCGGCACTTACATGTGCTCTCTATGTCGAgg                                                                                                                                         >  3:324044/1‑149 (MQ=255)
           aaaCTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAgg                                                                                                                                         >  3:180508/1‑148 (MQ=255)
              ctctGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGTCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCaaa                                                                                                                                     <  6:86269/149‑1 (MQ=255)
               tctGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCCCTCTATGTCGAGGCAAAg                                                                                                                                    >  3:314341/1‑149 (MQ=255)
               tctGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCCCTCTATGTCGAGGCAAAg                                                                                                                                    >  7:249814/1‑149 (MQ=255)
                  gtgCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCAGGTGTAGTGGTATTGGCGGATTAGGGCAAGGCTGCGGTCATTTCATGTAACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTc                                                                                                                                  >  7:20128/1‑148 (MQ=255)
                   tgCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCCCTCTATGTCGAGGCAAAGTCTa                                                                                                                                <  8:249814/149‑1 (MQ=255)
                   tgCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCCCTCTATGTCGAGGCAAAGTCTa                                                                                                                                <  4:314341/149‑1 (MQ=255)
                    gCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAg                                                                                                                               >  1:2216/1‑149 (MQ=255)
                     cccGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAg                                                                                                                               <  1:183668/148‑1 (MQ=255)
                      ccGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGc                                                                                                                              >  8:278920/1‑148 (MQ=255)
                      ccGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGc                                                                                                                              >  6:64104/1‑148 (MQ=255)
                                         aaTGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGaa                                                                                                          <  5:64104/149‑1 (MQ=255)
                                           tGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACt                                                                                                        >  7:2968/1‑149 (MQ=255)
                                               aTCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGGATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCt                                                                                                              <  1:299363/139‑1 (MQ=255)
                                               aTCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCt                                                                                                              >  2:299363/1‑139 (MQ=255)
                                                tctcGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGt                                                                                                   >  1:309706/1‑149 (MQ=255)
                                                  tcGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCAAGCCAGCAGCACATTGGCTGGAACTTTCgat                                                                                                  >  8:331783/1‑146 (MQ=255)
                                                  tcGGCCATAGGCATCAGATGGGTCCCGGTGTAGTAGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTg                                                                                                 >  5:40902/1‑149 (MQ=255)
                                                    ggCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATATGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTGTGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGt                                                                                                   <  4:296245/145‑1 (MQ=255)
                                                    ggCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTgc                                                                                                <  8:334902/148‑1 (MQ=255)
                                                     gCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTgcca                                                                                              <  2:155172/149‑1 (MQ=255)
                                                      ccATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTgccag                                                                                             >  5:370219/1‑149 (MQ=255)
                                                                 aGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTtgat                                                                                   <  7:278920/148‑1 (MQ=255)
                                                                  gATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTtgatga                                                                                 <  6:135083/149‑1 (MQ=255)
                                                                          ccGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGa                                                                         <  2:312452/149‑1 (MQ=255)
                                                                          ccGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGa                                                                         <  4:180508/149‑1 (MQ=255)
                                                                           cGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGc                                                                           >  3:81904/1‑146 (MQ=255)
                                                                           cGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGa                                                                         <  4:81904/148‑1 (MQ=255)
                                                                                  gTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTg                                                                  <  6:370219/148‑1 (MQ=255)
                                                                                    ggTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCCCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGAc                                                                                >  8:83108/1‑132 (MQ=255)
                                                                                    ggTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCCCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGAc                                                                                <  7:83108/132‑1 (MQ=255)
                                                                                        ttGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGaa                                                           >  5:25334/1‑149 (MQ=255)
                                                                                        ttGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGaa                                                           <  8:2968/149‑1 (MQ=255)
                                                                                           gCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACgg                                                        >  7:35945/1‑149 (MQ=255)
                                                                                             ggAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCAAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTc                                                      <  5:261026/149‑1 (MQ=255)
                                                                                              gAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGt                                                       >  7:2689/1‑147 (MQ=255)
                                                                                                   gggCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCagag                                                  >  2:163014/1‑147 (MQ=255)
                                                                                                       aaGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCTGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCaa                                            <  3:325719/149‑1 (MQ=255)
                                                                                                                   aTTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAgg                                >  5:37660/1‑149 (MQ=255)
                                                                                                                          gTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGcgccgc                          >  6:18947/1‑148 (MQ=255)
                                                                                                                             aCCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAAc                      <  5:18947/149‑1 (MQ=255)
                                                                                                                              ccTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCTTTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACa                     <  4:324044/149‑1 (MQ=255)
                                                                                                                                tCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGGACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAgg                   >  4:98263/1‑149 (MQ=255)
                                                                                                                                  ggCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAgg                   <  6:37660/147‑1 (MQ=255)
                                                                                                                                   gCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTa                <  2:2216/149‑1 (MQ=255)
                                                                                                                                       ttACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTAag              >  6:128703/1‑147 (MQ=255)
                                                                                                                                         aCATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCagag                                                  >  5:122726/1‑109 (MQ=255)
                                                                                                                                         aCATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCagag                                                  <  6:122726/109‑1 (MQ=255)
                                                                                                                                         aCATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAgg                   <  3:44780/140‑1 (MQ=255)
                                                                                                                                         aCATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAgg                   >  4:44780/1‑140 (MQ=255)
                                                                                                                                         aCATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTAagca            <  8:2689/147‑1 (MQ=255)
                                                                                                                                                 tctctATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTAAGCAGCCCGGTCAg  >  3:131510/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCGGAACATGAAACTCTGTGCCCGGGTACGCCAGACGCACAATGGAATCTCGGCCATAGGCATCAGATGGGTCCCGGTGTAGTGGTATTGGCGGAATAGGGCAAGGCTGCGGTCATTGCATGTCACCTCGGCACTTACATGTGCTCTCTATGTCGAGGCAAAGTCTAGCCAGCAGCCCATTGGCTGGAACTTGCGTTGCCAGCCAGCTTGATGACTGAGCGAGTCAGTGGCACGAACGGTCAGAGCAGCAAGGCCAACAAAGGCGCCGCGAACAGGTTAAGCAGCCCGGTCAG  >  NC_002947/5249758‑5250051

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