Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,195 G→C 5.9% A56A (GCC→GCG) ‡ TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,1950GC 5.9% 276.5 / 6.5 85A56A (GCC→GCG) ‡TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (35/45);  new base C (3/2);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

TTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGC  >  NC_002947/5297048‑5297334
                                                                                                                                                    |                                                                                                                                           
ttaccAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCggg                                                                                                                                            <  3:87077/146‑1 (MQ=255)
ttGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCggg                                                                                                                                            >  3:118246/1‑149 (MQ=255)
    cAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                        >  1:225693/1‑149 (MQ=255)
    cAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                        >  1:25498/1‑149 (MQ=255)
    cAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                        <  7:258895/149‑1 (MQ=255)
     aGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                                                                        >  5:5870/1‑148 (MQ=255)
     aGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                       <  3:118713/149‑1 (MQ=255)
      gTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                                                                                       >  4:178885/1‑148 (MQ=255)
       tGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCaa                                                                                                                                     <  7:217250/149‑1 (MQ=255)
       tGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCaa                                                                                                                                     <  1:166996/149‑1 (MQ=255)
         cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACt                                                                                                                                   <  6:207513/149‑1 (MQ=255)
          gcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTg                                                                                                                                  >  6:136248/1‑149 (MQ=255)
            cAGTTGCTCGCCGTCGTTGGTGAAGAGCTTGGCCTGGGCCACTTCTTGGGTTCTGGCTATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTATTTGATCGAGCGGGCGATCAACTGCa                                                                                                                                <  5:187551/149‑1 (MQ=255)
            cAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCa                                                                                                                                <  6:145051/149‑1 (MQ=255)
            cAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCa                                                                                                                                <  4:118246/149‑1 (MQ=255)
                     gCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGt                                                                                                                       >  3:274022/1‑149 (MQ=255)
                      cccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGtt                                                                                                                      <  6:5870/149‑1 (MQ=255)
                       ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCg                                                                                                                                          <  4:303208/128‑1 (MQ=255)
                       ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                     <  5:136248/149‑1 (MQ=255)
                        cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                    <  8:198746/149‑1 (MQ=255)
                        cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                    >  6:369865/1‑149 (MQ=255)
                         tCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                    >  6:204703/1‑148 (MQ=255)
                           gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCAt                                                                                                                 <  7:80498/149‑1 (MQ=255)
                           gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCAt                                                                                                                 <  4:274022/149‑1 (MQ=255)
                            ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                                                                              >  4:250258/1‑135 (MQ=255)
                            ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                                                                              <  3:250258/135‑1 (MQ=255)
                             tgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCa                                                                                                                  >  5:190376/1‑146 (MQ=255)
                              gttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                                                                               <  8:322516/148‑1 (MQ=255)
                              gttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGa                                                                                                              >  8:51579/1‑149 (MQ=255)
                              gttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGATATGGTTGGCATAGa                                                                                                              >  5:251282/1‑149 (MQ=255)
                                 gAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGaca                                                                                                            <  3:80304/148‑1 (MQ=255)
                                  aaGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACg                                                                                                          <  2:225693/149‑1 (MQ=255)
                                    gACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCa                                                                                                        <  2:25498/149‑1 (MQ=255)
                                      ccTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTCTCCTCGAACACCACAGCAG‑CTTCTACGCTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAgg                                                                                                      >  5:10691/1‑149 (MQ=255)
                                       cTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                     >  6:384456/1‑149 (MQ=255)
                                        gatctcTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                                                                    <  3:337656/144‑1 (MQ=255)
                                          ggCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                                                                    <  5:384456/147‑1 (MQ=255)
                                           gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGCCCGCTCGATCAACta                                                                                                                                  <  4:324641/116‑2 (MQ=255)
                                           gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGCCCGCTCGATCAACta                                                                                                                                  >  3:324641/1‑115 (MQ=255)
                                             cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                               >  2:243145/1‑149 (MQ=255)
                                             cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                               >  6:212924/1‑149 (MQ=255)
                                              tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCTGGGCGTGctc                                                                                              >  2:2815/1‑147 (MQ=255)
                                              tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                                                                    >  4:337656/1‑143 (MQ=255)
                                                 gCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                                                                           >  5:83368/1‑149 (MQ=255)
                                                   cACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                                                                                         >  2:451/1‑149 (MQ=255)
                                                     cttcttCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGTGCTCGCCctgc                                                                                       <  7:51579/149‑1 (MQ=255)
                                                        cttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagca                                                                                      <  5:204703/147‑1 (MQ=255)
                                                         ttcGATACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAggg                                                                                   >  2:285885/1‑149 (MQ=255)
                                                            ggTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCTGGGCGTGCTCGCCctgcgtgtct                                                                                 <  1:2815/148‑7 (MQ=255)
                                                             gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                               >  3:253742/1‑149 (MQ=255)
                                                              tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              <  6:231645/149‑1 (MQ=255)
                                                              tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              <  5:212924/149‑1 (MQ=255)
                                                                cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                            <  5:353302/149‑1 (MQ=255)
                                                                 tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGac                                                                                                             >  3:317585/1‑115 (MQ=255)
                                                                 tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGac                                                                                                             <  4:317585/115‑1 (MQ=255)
                                                                 tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  4:55235/1‑146 (MQ=255)
                                                                 tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  5:195802/1‑146 (MQ=255)
                                                                 tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  2:68908/1‑146 (MQ=255)
                                                                   gCAATCACCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCAG‑CTTCTACGGTGTAGTTGATCGCGCGGGCGATCAACTGCAGGTGATGGTTGGCAGAGACACGCAGGGCGCGCAGGCGCAGCAGGGTGAGGgaga                                                                          >  8:167397/1‑148 (MQ=255)
                                                                     aaTCAACGCCAGCGATTCGTTGGCGCCGGCGCAGTCATCGAACACCACAGCCGCCTTCTACGGTCTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGATCACGACCAGCAGGGAGAGGGAGATa                                                                        >  1:290243/1‑149 (MQ=255)
                                                                           cGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCa                                                                 <  1:68908/149‑1 (MQ=255)
                                                                            gCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg                                                                <  5:369865/149‑1 (MQ=255)
                                                                              cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                              <  3:55235/149‑1 (MQ=255)
                                                                                  gATTCGTTGGCGGCGTCGCTGGCATCGAACACGACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAAc                                                          <  6:251282/149‑1 (MQ=255)
                                                                                   aTTCGTTGGCGGCGGCGCTGGCCGCGAACCCCACAGCCG‑CGTCTACGGGGTAGGTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                         <  6:131407/149‑1 (MQ=255)
                                                                                    ttCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACgg                                                        >  1:241944/1‑149 (MQ=255)
                                                                                      cGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgc                                                       <  4:253742/148‑1 (MQ=255)
                                                                                           gcggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc                                                                   <  2:228084/131‑1 (MQ=255)
                                                                                           gcggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc                                                                   >  1:228084/1‑131 (MQ=255)
                                                                                             ggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAgg                                               <  1:285885/149‑1 (MQ=255)
                                                                                                     tGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCt                                        >  3:2208/1‑148 (MQ=255)
                                                                                                       tCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                     >  4:208647/1‑149 (MQ=255)
                                                                                                       tCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                     <  1:243145/149‑1 (MQ=255)
                                                                                                         aTCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                   <  3:178885/149‑1 (MQ=255)
                                                                                                          tCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTg                                  <  4:2208/149‑1 (MQ=255)
                                                                                                            gAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTgg                                 <  1:451/148‑1 (MQ=255)
                                                                                                                  cacaGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGCTACCCGCCAGGc                                                              >  3:136984/1‑113 (MQ=255)
                                                                                                                    caGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                    <  6:266218/89‑1 (MQ=255)
                                                                                                                    caGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAAGCCCAGCAgg                                                                                    >  5:266218/1‑89 (MQ=255)
                                                                                                                              cTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGa                <  6:195802/148‑1 (MQ=255)
                                                                                                                                    tgtAGTTGATCGAGCGGGGGATCAACTGCTTGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCGCAGCAGGGTGAGGGAGATACCCGCCTGGCCAACGGCAAGCAAGGCGACGCTGAGTTGGCCACGGCCCTGCACCGAGCGCAgg         <  2:290243/149‑1 (MQ=255)
                                                                                                                                    tgcaGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg         >  6:153701/4‑149 (MQ=255)
                                                                                                                                        gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                <  2:47529/106‑1 (MQ=255)
                                                                                                                                        gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                >  1:47529/1‑106 (MQ=255)
                                                                                                                                         ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                               <  6:22409/74‑1 (MQ=38)
                                                                                                                                         ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                               >  5:22409/1‑74 (MQ=38)
                                                                                                                                            aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc  <  6:83368/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
TTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGC  >  NC_002947/5297048‑5297334

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