Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,302 G→C 6.9% A39G (GCC→GGC) ‡ PP_1615 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,812,3020GC 6.9% 180.8 / 3.5 58A39G (GCC→GGC) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812155‑1812445
                                                                                                                                                   |                                                                                                                                                
gCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCggg                                                                                                                                                  <  5:148054/148‑1 (MQ=255)
     gTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCt                                                                                                                                            <  4:224703/149‑1 (MQ=255)
      tCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCATCCAGCCGGGCGATCTg                                                                                                                                           >  1:157394/1‑149 (MQ=255)
        gCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGcc                                                                                                                                         >  7:187919/1‑149 (MQ=255)
             tGGAGCGCATCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTTCAGGAAGGTCTGGCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTg                                                                                                                                    <  4:184675/149‑1 (MQ=255)
              ggCGCGCAGCGCTGCCTCCCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCAGGCTGCGCTCAGGGGTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTgg                                                                                                                                   <  5:217020/149‑2 (MQ=255)
               gcgcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGAt                                                                                                                                  >  2:208155/1‑149 (MQ=255)
                    cAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCCGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCg                                                                                                                             <  7:118705/149‑1 (MQ=255)
                     aGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGa                                                                                                                            <  1:125975/149‑1 (MQ=255)
                      gcgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGAt                                                                                                                           <  5:228265/149‑1 (MQ=255)
                      gcgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGAt                                                                                                                           <  6:221048/149‑1 (MQ=255)
                       cgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGa                                                                                                                            >  7:215753/1‑147 (MQ=255)
                       cgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGa                                                                                                                            <  8:215753/147‑1 (MQ=255)
                         cTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAg                                                                                                                        <  3:59641/149‑1 (MQ=255)
                              tCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                   >  3:244757/1‑149 (MQ=255)
                              tCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                   >  5:230183/1‑149 (MQ=255)
                                aCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAccc                                                                                                                 >  3:17569/1‑149 (MQ=255)
                                aCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAccc                                                                                                                 <  8:154533/149‑1 (MQ=255)
                                                 gTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTcggcggcg                                                                                               <  4:244757/149‑1 (MQ=255)
                                                     gCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAg                                                                                           >  6:203444/1‑149 (MQ=255)
                                                      cAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGt                                                                                          <  1:208155/149‑1 (MQ=255)
                                                              cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGc                                                                                  <  3:100528/149‑1 (MQ=255)
                                                              cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGc                                                                                  >  5:301995/1‑149 (MQ=255)
                                                                ggAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGccc                                                                                >  7:172344/1‑149 (MQ=255)
                                                                ggAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTgg                                                                                                                                   >  7:71825/1‑98 (MQ=255)
                                                                ggAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTgg                                                                                                                                   <  8:71825/99‑2 (MQ=255)
                                                                      tCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCg                                                                          >  8:79426/1‑149 (MQ=255)
                                                                       cTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCgg                                                                         <  8:264079/149‑1 (MQ=255)
                                                                          cccGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGc                                                                        <  1:236656/147‑1 (MQ=255)
                                                                               cTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGc                                                                  >  8:60739/1‑148 (MQ=255)
                                                                                ttCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGa                                                                >  6:172913/1‑149 (MQ=255)
                                                                                     tCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCGCGATTCGGCCACGGCGCCAAAt                                                           >  3:34056/1‑149 (MQ=255)
                                                                                     tCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAAt                                                           >  6:44547/1‑149 (MQ=255)
                                                                                       aGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtc                                                         <  4:68585/149‑1 (MQ=255)
                                                                                          gTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGACCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtcgtc                                                      >  1:116380/1‑149 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGa                                                    <  5:172913/148‑1 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGa                                                    <  7:79426/148‑1 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGa                                                    <  8:172344/148‑1 (MQ=255)
                                                                                               ccGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGc                                                 >  2:168679/1‑149 (MQ=255)
                                                                                                 ggcggGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCt                                               <  5:44547/149‑1 (MQ=255)
                                                                                                       tgctgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCAt                                         <  6:230183/149‑1 (MQ=255)
                                                                                                         ctgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCg                                       <  2:116380/149‑1 (MQ=255)
                                                                                                            cGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgc                                     <  5:203444/148‑1 (MQ=255)
                                                                                                             gATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctg                                   >  2:79747/1‑149 (MQ=255)
                                                                                                              aTACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctgc                                  <  6:301995/149‑1 (MQ=255)
                                                                                                                  agcagGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCAt                              <  7:60739/149‑1 (MQ=255)
                                                                                                                            tttGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATccc                    >  6:51054/1‑149 (MQ=255)
                                                                                                                               gCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATcc                     <  5:174288/145‑1 (MQ=255)
                                                                                                                               gCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTAATCGATTATcc                     >  6:174288/1‑145 (MQ=255)
                                                                                                                                cAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgtg                >  2:264224/1‑149 (MQ=255)
                                                                                                                                cAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgt                 <  1:79747/148‑1 (MQ=255)
                                                                                                                                   gCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTgg             >  8:186118/1‑149 (MQ=255)
                                                                                                                                     ggTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGc            >  2:246709/1‑148 (MQ=255)
                                                                                                                                     ggTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGAAAGCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGc            >  1:257171/1‑148 (MQ=255)
                                                                                                                                        acatcAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGTGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGCCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCTTCGATTATCCCTGTGGGGCaa          >  1:226247/5‑148 (MQ=255)
                                                                                                                                          tccAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat      >  5:45902/1‑149 (MQ=255)
                                                                                                                                          tccAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat      >  8:292298/1‑149 (MQ=255)
                                                                                                                                            cAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGt    <  7:189962/149‑1 (MQ=255)
                                                                                                                                              gCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGttt  >  3:120946/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812155‑1812445

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