Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,192 G→C 5.6% R57R (CGC→CGG)  TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,1920GC 5.6% 295.6 / 6.6 89R57R (CGC→CGG) TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (37/47);  new base C (3/2);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGC  >  NC_002947/5297046‑5297334
                                                                                                                                                   |                                                                                                                                              
cgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCg                                                                                                                                              <  3:140703/149‑1 (MQ=255)
cgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACAGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCg                                                                                                                                              >  6:353302/1‑149 (MQ=255)
  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                                                                                                                                     <  1:166996/149‑1 (MQ=255)
         tGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCaa                                                                                                                                     <  7:217250/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                                                                                                                 <  4:274022/149‑1 (MQ=255)
                             gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCAt                                                                                                                 <  7:80498/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                                                                                                                                  >  3:324641/1‑115 (MQ=255)
                                             gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGCCCGCTCGATCAACta                                                                                                                                  <  4:324641/116‑2 (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                                                                              <  5:212924/149‑1 (MQ=255)
                                                                tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              <  6:231645/149‑1 (MQ=255)
                                                                  cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                            <  5:353302/149‑1 (MQ=255)
                                                                   tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGac                                                                                                             <  4:317585/115‑1 (MQ=255)
                                                                   tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGac                                                                                                             >  3:317585/1‑115 (MQ=255)
                                                                   tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  2:68908/1‑146 (MQ=255)
                                                                   tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  4:55235/1‑146 (MQ=255)
                                                                   tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                              >  5:195802/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                                                >  1:47529/1‑106 (MQ=255)
                                                                                                                                          gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                <  2:47529/106‑1 (MQ=255)
                                                                                                                                           ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                               >  5:22409/1‑74 (MQ=38)
                                                                                                                                           ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                               <  6:22409/74‑1 (MQ=38)
                                                                                                                                              aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc  <  6:83368/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGC  >  NC_002947/5297046‑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: ‑