Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,309 G→C 8.1% Q37E (CAG→GAG) ‡ 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,3090GC 8.1% 186.9 / 6.5 62Q37E (CAG→GAG) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base C (2/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

CCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTC  >  NC_002947/1812164‑1812449
                                                                                                                                                    |                                                                                                                                            
cccATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCg                                                                                                                                           >  4:255843/1‑149 (MQ=255)
    tGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTg                                                                                                                                       <  3:111631/149‑1 (MQ=255)
     ggCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGa                                                                                                                                      >  7:289894/1‑149 (MQ=255)
      gcgcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TACAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTgga                                                                                                                                     <  2:39496/149‑3 (MQ=255)
       cgcgcAGCGCGTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTggag                                                                                                                                    <  4:248544/149‑4 (MQ=255)
        gcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                        <  5:367244/144‑1 (MQ=255)
        gcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                        >  6:367244/1‑144 (MQ=255)
               gcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGGAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGAt                                                                                                                              >  6:277301/1‑147 (MQ=255)
                cTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAg                                                                                                                           <  3:189622/149‑1 (MQ=255)
                    cTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTGGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCATGCGCGTTCTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGa                                                                                                                         <  1:280001/148‑1 (MQ=255)
                       tcGGCCATCCTGCAAGCGTTGGCAATGC‑TGCATTAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGTGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCAAGCCGTGCGATCTGCCGGCTGATGTGCGATGATGAAAc                                                                                                                      <  1:161695/146‑1 (MQ=255)
                           ccATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAg                                                                                                                >  7:316851/1‑149 (MQ=255)
                            cttCCTGC‑AGCATTGCCAACGCTTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAgg                                                                                                               >  3:361747/3‑149 (MQ=255)
                              tCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCg                                                                                                             >  2:217694/1‑149 (MQ=255)
                               ccTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGt                                                                                                            >  7:292052/1‑149 (MQ=255)
                                   cAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTcgg                                                                                                        >  2:144457/1‑149 (MQ=255)
                                   cAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTcgg                                                                                                        <  5:37813/149‑1 (MQ=255)
                                            gCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCa                                                                                                <  7:271731/148‑1 (MQ=255)
                                                 gc‑tgcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGAGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAAc                                                                                          >  2:40245/1‑149 (MQ=255)
                                                    tgcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTg                                                                                        >  5:339325/1‑149 (MQ=255)
                                                    gaccattaGGTCTGCCGGGCTTCGCTCAGGGTTACCGGGCGGGGGCTGG‑GATACTGCAGGCGGGTT‑TGCCGGCGGGCCTCAAGCCGGGCGAGCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACt                                                                                         <  1:91770/141‑1 (MQ=255)
                                                     gcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTgg                                                                                       <  8:289894/149‑1 (MQ=255)
                                                      cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGACACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGc                                                                                      >  4:324371/1‑149 (MQ=255)
                                                      cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGc                                                                                      >  7:139426/1‑149 (MQ=255)
                                                       aGGAAGGTCTGACCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGGCGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAATCCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGcc                                                                                     >  8:80564/1‑149 (MQ=255)
                                                        ggAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGccc                                                                                    <  1:217694/149‑1 (MQ=255)
                                                           aGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGAt                                                                                 >  1:84681/1‑149 (MQ=255)
                                                                  cccGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCAGATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGc                                                                            >  1:244916/1‑148 (MQ=255)
                                                                   ccGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCAc                                                                         <  8:316851/149‑1 (MQ=255)
                                                                    cGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAACCCCCCAGGCGTGCGGAGGCGGCAGTGACCTGGCCCGATTCGGCCACg                                                                        >  8:271902/1‑149 (MQ=255)
                                                                          cGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGa                                                                    <  3:255843/147‑1 (MQ=255)
                                                                             tCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAAt                                                               >  7:323923/1‑149 (MQ=255)
                                                                             tCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAAt                                                               >  3:258442/1‑149 (MQ=255)
                                                                               aGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGGAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtc                                                             <  5:277301/149‑1 (MQ=255)
                                                                                 ggTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACg                                                                        >  1:158354/1‑136 (MQ=255)
                                                                                 ggTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACg                                                                        <  2:158354/136‑1 (MQ=255)
                                                                                  gTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtcgtc                                                          <  6:339325/149‑1 (MQ=255)
                                                                                  gTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtcgtc                                                          <  2:84681/149‑1 (MQ=255)
                                                                                               tgctgc‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCAt                                             >  3:255014/1‑149 (MQ=255)
                                                                                               tgctgc‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCAt                                             >  3:250142/1‑149 (MQ=255)
                                                                                                 ctgc‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCg                                           >  8:94249/1‑149 (MQ=255)
                                                                                                  tgc‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                          >  3:359496/1‑149 (MQ=255)
                                                                                                    c‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                          <  1:192034/147‑1 (MQ=255)
                                                                                                    c‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgct                                        <  8:139426/149‑1 (MQ=255)
                                                                                                    c‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgct                                        >  4:327262/1‑149 (MQ=255)
                                                                                                      gATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctg                                       >  8:105237/1‑149 (MQ=255)
                                                                                                       aTACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctgc                                      <  7:94249/149‑1 (MQ=255)
                                                                                                             cagGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCaac                                 >  4:116748/1‑146 (MQ=255)
                                                                                                             cagGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGGTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCg                                <  8:323923/149‑1 (MQ=255)
                                                                                                               ggCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGAt                              <  4:255014/149‑1 (MQ=255)
                                                                                                                gcgcGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGAt                              <  4:250142/148‑1 (MQ=255)
                                                                                                                  gcgTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTAt                           >  1:149950/1‑149 (MQ=255)
                                                                                                                          cAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgt                     <  1:144457/148‑1 (MQ=255)
                                                                                                                           aGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATc                                 <  4:71144/135‑1 (MQ=255)
                                                                                                                           aGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgtgt                   <  6:42244/149‑1 (MQ=255)
                                                                                                                           aGGCGGTCCTCCAGCCGGGCGATCTGCCGCCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATc                                 >  3:71144/1‑135 (MQ=255)
                                                                                                                            ggCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgtgtg                  >  2:153620/1‑149 (MQ=255)
                                                                                                                                  acatcAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaat             >  7:96937/5‑148 (MQ=255)
                                                                                                                                    tccAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat          >  7:53884/1‑149 (MQ=255)
                                                                                                                                    atcAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat          >  1:599/3‑149 (MQ=255)
                                                                                                                                       aGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGtt       >  3:262441/1‑149 (MQ=255)
                                                                                                                                            ggCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTc  >  2:34111/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGC‑TGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGC‑GATACAGCAGGCGCGTT‑TGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTC  >  NC_002947/1812164‑1812449

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑