Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,715,396 C→G 8.3% G21G (GGC→GGG)  tatB → Sec‑independent protein translocase protein TatB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,715,3960CG 8.3% 186.1 / 6.3 60G21G (GGC→GGG) tatBSec‑independent protein translocase protein TatB
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (2/3);  total (29/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 = 3.35e-01

CACCGCGCAGCCGCAGCAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATCCTGCA  >  NC_002947/5715255‑5715542
                                                                                                                                             |                                                                                                                                                   
cacCGCGCAGCCGCAGCAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAgc                                                                                                                                              <  6:8430/149‑1 (MQ=255)
  ccGCGCAGCCGCAGCAAGGCCAAACAAACGAAGGCAAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAgcgc                                                                                                                                            >  6:111645/1‑149 (MQ=255)
       cagccgcagcAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCACCTGcc                                                                                                                                       >  6:287021/1‑149 (MQ=255)
        agccgcagcAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGGGCCGAGCACCAGcag                                                                                                                                      <  1:152125/149‑3 (MQ=255)
        acccgcagcAAGG‑CACCCCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCAGGCGCTCGGGGCCGAGCACCAGCAgg                                                                                                                                     <  1:292294/147‑1 (MQ=255)
          ccgcagcAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCggg                                                                                                                                    >  6:226473/1‑149 (MQ=255)
            gcagcaAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTg                                                                                                                                  >  7:121457/1‑149 (MQ=255)
              agcaAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGcc                                                                                                                                >  2:123840/1‑149 (MQ=255)
               gcaAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAGGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCcg                                                                                                                               >  8:34909/1‑149 (MQ=255)
                     ccACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑a                                                                                                                        <  6:291847/149‑1 (MQ=255)
                           cATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCGGCAGGCGCTCGGGGCCGAGCACCAGCAggg                                                                                                                                    <  8:90404/132‑1 (MQ=38)
                           cATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCGGCAGGCGCTCGGGGCCGAGCACCAGCAggg                                                                                                                                    >  7:90404/1‑132 (MQ=38)
                             tCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAggg                                                                                                                <  1:195770/149‑1 (MQ=255)
                                aaGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtg                                                                                                             >  6:293446/1‑149 (MQ=255)
                                 aggccaggccCACCCGGTCCAGGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAggcgctc                                                                                                            <  7:34909/149‑6 (MQ=255)
                                 aggccaggccCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            <  3:18279/149‑1 (MQ=255)
                                 aggccaggccCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            <  6:189014/149‑1 (MQ=255)
                                   gccaggccCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTgg                                                                                                          <  7:83762/149‑1 (MQ=255)
                                     caggccCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGAt                                                                                                        >  1:168418/1‑149 (MQ=255)
                                         cccACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATcggc                                                                                                    <  5:293446/149‑1 (MQ=255)
                                          ccACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCa                                                                                                                   >  2:253434/1‑133 (MQ=255)
                                          ccACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCa                                                                                                                   <  1:253434/133‑1 (MQ=255)
                                                    cAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACg                                                                                         >  3:82429/1‑149 (MQ=255)
                                                       gagCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            >  3:251009/1‑127 (MQ=255)
                                                       gagCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            <  4:251009/127‑1 (MQ=255)
                                                        agCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGc                                                                                     <  5:305418/149‑1 (MQ=255)
                                                                 ggAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAg                                                                             >  3:230160/1‑148 (MQ=255)
                                                                   aaaGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGAACGGCCGGCTCAAACGCAGCTTCAACAGCa                                                                           >  7:138108/1‑148 (MQ=255)
                                                                    aaGACTGACCCATGTTCGGCATGAGTTTCAGCGAGTTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATc                                                                         <  6:283731/149‑1 (MQ=255)
                                                                    aaGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATc                                                                         <  2:38157/149‑1 (MQ=255)
                                                                       aCTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATc                                                                         >  4:49955/1‑146 (MQ=255)
                                                                           aCCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            >  4:75933/1‑107 (MQ=255)
                                                                           aCCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCtgt                                                                                                            <  3:75933/107‑1 (MQ=255)
                                                                               aTGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGt                                                              >  8:302148/1‑149 (MQ=255)
                                                                                tGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTg                                                             <  7:302148/149‑1 (MQ=255)
                                                                                 gTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTgg                                                            >  3:102093/1‑149 (MQ=255)
                                                                                  ttCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGTAAGTg                                                             >  7:116263/1‑147 (MQ=255)
                                                                                    cGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAgc                                                         >  5:163334/1‑149 (MQ=255)
                                                                                             tttCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGa                                                           >  4:226973/1‑138 (MQ=255)
                                                                                             tttCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGa                                                           <  3:226973/138‑1 (MQ=255)
                                                                                                 aGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAgcgc                                                       <  3:38684/138‑1 (MQ=255)
                                                                                                 aGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAgcgc                                                       >  4:38684/1‑138 (MQ=255)
                                                                                                        tgctgctCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCgacgac                                     <  8:121457/149‑1 (MQ=255)
                                                                                                             ctCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCACCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATccg                                <  5:287021/149‑1 (MQ=255)
                                                                                                                 tcgGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGCAATAg                                                <  7:292785/129‑1 (MQ=255)
                                                                                                                 tcgGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAg                                                >  8:292785/1‑129 (MQ=255)
                                                                                                                 tcgGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCcgccgc                            <  1:214282/149‑1 (MQ=255)
                                                                                                                    gCCTGGTCGCCCTGCT‑GTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAAc                        >  6:54218/1‑149 (MQ=255)
                                                                                                                      cTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACt                       <  1:123840/149‑1 (MQ=255)
                                                                                                                        ggTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGCAACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTg                      >  1:80148/1‑149 (MQ=255)
                                                                                                                         gTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCTC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGc                     <  8:138108/148‑1 (MQ=255)
                                                                                                                            gCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATaa                 >  2:214247/1‑149 (MQ=255)
                                                                                                                            gCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATa                  <  6:163334/148‑1 (MQ=255)
                                                                                                                             cccGGCAGGCGCTCGGGGCCGAGCACCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAAc                >  3:293967/1‑149 (MQ=255)
                                                                                                                                      tGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATc       <  2:112740/149‑1 (MQ=255)
                                                                                                                                         tCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAAc                        <  1:77678/129‑1 (MQ=255)
                                                                                                                                         tCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAAc                        >  2:77678/1‑129 (MQ=255)
                                                                                                                                           ggCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATCCTg    >  5:218768/1‑147 (MQ=255)
                                                                                                                                           ggCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATCCTg    <  4:102093/147‑1 (MQ=255)
                                                                                                                                           ggCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATCCTGCa  <  3:49955/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CACCGCGCAGCCGCAGCAAGGCCACACCATCGAAGGCCAGGCCCACCCGGTCCAAGAGCCGCAGCGGAAAGACTGACCCATGTTCGGCATGAGTTTCAGCGAGCTGCTGCTCGTCGGCCTGGTCGCCCTGCTGGTGCTCGGCCCCGAGCGCCTGCCGGGTGCCGCGCGC‑ACGGCAGGGCTGTGGATCGGCCGGCTCAAACGCAGCTTCAACAGCATCAAGATGGAAGTGGAGCGCGAAATAGGCGCCGACGACATCCGCCGCCAACTGCATAACGAGCACATCCTGCA  >  NC_002947/5715255‑5715542

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