Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,303 C→A 6.9% A39S (GCC→TCC) ‡ 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,3030CA 6.9% 182.7 / 3.7 58A39S (GCC→TCC) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base A (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 = 3.24e-01

GTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812160‑1812445
                                                                                                                                               |                                                                                                                                               
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                                                                                                                           <  6:221048/149‑1 (MQ=255)
                 gcgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGAt                                                                                                                           <  5:228265/149‑1 (MQ=255)
                  cgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGa                                                                                                                            <  8:215753/147‑1 (MQ=255)
                  cgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGa                                                                                                                            >  7:215753/1‑147 (MQ=255)
                    cTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAg                                                                                                                        <  3:59641/149‑1 (MQ=255)
                         tCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                   >  5:230183/1‑149 (MQ=255)
                         tCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                   >  3:244757/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                                                                                                                                   <  8:71825/99‑2 (MQ=255)
                                                           ggAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTgg                                                                                                                                   >  7:71825/1‑98 (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                                                    <  8:172344/148‑1 (MQ=255)
                                                                                        aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGa                                                    <  7:79426/148‑1 (MQ=255)
                                                                                        aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGa                                                    <  5:172913/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)
                                                                                                                                               |                                                                                                                                               
GTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAG‑GCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812160‑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: ‑