Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,136,218 A→C 8.7% S143R (AGC→CGC)  PP_3640 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,136,2180AC 8.7% 210.6 / 9.1 69S143R (AGC→CGC) PP_3640AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base C (2/4);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCAGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGC  >  NC_002947/4136073‑4136365
                                                                                                                                                  |                                                                                                                                                   
ccAGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGgcagcg                                                                                                                                                  >  3:210677/1‑149 (MQ=255)
ccAGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGgcagcg                                                                                                                                                  >  1:68413/1‑148 (MQ=255)
  aGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGgcagcgca                                                                                                                                                <  8:227649/149‑1 (MQ=255)
  aGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGgcagcgca                                                                                                                                                <  3:136307/149‑1 (MQ=255)
   ggCCATGCAGGCGATGGATAGCCAGGCCGCCGACGGCGTTGGCATGGTCTGCGGCTTCTTCCATTTTCAGTGGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGGTACGCGATGAAGATCACAGCGgcagctcag                                                                                                                                               >  7:353880/1‑149 (MQ=255)
    gCCATGCAGGCGCTGGATAGCCAGCCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGgcagcacagc                                                                                                                                              >  2:272625/1‑149 (MQ=255)
      cATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCagcgcagcgc                                                                                                                                            >  1:249042/1‑149 (MQ=255)
      cATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCagcgcagcgc                                                                                                                                            >  6:83064/1‑149 (MQ=255)
           aGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGgcg                                                                                                                                       <  3:26426/149‑1 (MQ=255)
             gcgcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGTATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCt                                                                                                                                     <  6:243649/149‑1 (MQ=255)
               gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGt                                                                                                                                   >  3:183188/1‑149 (MQ=255)
               gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGt                                                                                                                                   >  3:168385/1‑149 (MQ=255)
               gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGGCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTgg                                                                                                                                   <  7:97611/149‑2 (MQ=255)
                cTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGtt                                                                                                                                  >  7:30636/1‑149 (MQ=255)
                   gATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCgg                                                                                                                               >  8:224878/1‑149 (MQ=255)
                       gCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCt                                                                                                                            <  8:38464/148‑1 (MQ=255)
                        ccAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGGACGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCTGGCt                                                                                                                            >  1:275003/1‑148 (MQ=255)
                        ccAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGa                                                                                                                          >  7:277868/1‑149 (MQ=255)
                         cAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTggga                                                                                                                                 >  1:335696/1‑138 (MQ=255)
                         cAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTggga                                                                                                                                 <  2:335696/141‑4 (MQ=255)
                          aGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATc                                                                                                                        >  6:174969/1‑149 (MQ=255)
                             ccgccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTg                                                                                                                     <  5:41918/149‑1 (MQ=255)
                              cgtcgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTgg                                                                                                                    <  5:289494/146‑1 (MQ=255)
                               gtcgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGga                                                                                                                   <  7:287802/147‑1 (MQ=255)
                               gccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGga                                                                                                                   <  5:307119/149‑1 (MQ=255)
                                ccgACGGCGTTGGCTTGGTCTGCGGGTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGTGGTTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTTTTCGGGCTGATCCTGgaa                                                                                                                  <  5:247714/149‑1 (MQ=255)
                                ccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGAATGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaa                                                                                                                  <  1:374447/149‑1 (MQ=255)
                                  gACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaaga                                                                                                                <  8:30636/149‑1 (MQ=255)
                                  gACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaag                                                                                                                 >  6:296494/1‑148 (MQ=255)
                                     ggCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATg                                                                                                             <  4:210677/149‑1 (MQ=255)
                                            gCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGc                                                                                                        >  4:211646/1‑147 (MQ=255)
                                            gCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCga                                                                                                      <  8:103395/149‑1 (MQ=255)
                                                 gTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGc                                                                                                 <  5:296494/149‑1 (MQ=255)
                                                  tCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGAGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCt                                                                                                >  2:300253/1‑149 (MQ=255)
                                                    tGCCGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTCAAGATCCAAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                              >  6:49305/1‑149 (MQ=255)
                                                      cGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                              <  5:83064/147‑1 (MQ=255)
                                                      cGGCTTCTTCCATTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                              <  2:68413/146‑1 (MQ=255)
                                                         cttcttCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCAt                                                                                         <  4:168385/149‑1 (MQ=255)
                                                          ttcttcCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCAt                                                                                         >  8:222375/1‑148 (MQ=255)
                                                            ctgccATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATca                                                                                       <  5:42780/145‑1 (MQ=255)
                                                               ccATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGTCATCACAg                                                                                    <  3:131251/148‑1 (MQ=255)
                                                                  ttttCAGTCGGGCCTCTGTTCGCTGATTGTCGAGGGTCGCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGc                                                                                <  4:107840/149‑1 (MQ=255)
                                                                     tCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGcc                                                                                           >  5:299637/1‑135 (MQ=255)
                                                                     tCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGcc                                                                                           <  6:299637/135‑1 (MQ=255)
                                                                      cAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTg                                                                                                                                    >  2:35705/1‑93 (MQ=38)
                                                                      cAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTg                                                                                                                                    <  1:35705/93‑1 (MQ=38)
                                                                       aGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACt                                                                            <  4:183188/148‑1 (MQ=255)
                                                                            gggCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAAc                                                                       <  5:174969/149‑1 (MQ=255)
                                                                            gggCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAAc                                                                       <  6:106536/149‑1 (MQ=255)
                                                                              gCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACg                                                                      >  3:268191/1‑148 (MQ=255)
                                                                               ccTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGc                                                                    <  7:164381/149‑1 (MQ=255)
                                                                                ctctGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCt                                                                   >  5:254700/1‑149 (MQ=255)
                                                                                        gCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACt                                                           >  5:258874/1‑149 (MQ=255)
                                                                                        gCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCACAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCAc                                                            <  1:272625/148‑1 (MQ=255)
                                                                                          tGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACttg                                                         >  3:273523/1‑149 (MQ=255)
                                                                                           gATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACttgt                                                        >  2:352421/1‑149 (MQ=255)
                                                                                              tGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTtgttgt                                                     >  3:61265/1‑149 (MQ=255)
                                                                                               gTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTtgttgt                                                     <  8:277868/148‑1 (MQ=255)
                                                                                                tCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTc                                                   <  4:268191/149‑1 (MQ=255)
                                                                                                  gAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCt                                                 <  1:300253/149‑1 (MQ=255)
                                                                                                    gggTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGTGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTg                                                >  2:177446/1‑148 (MQ=255)
                                                                                                              gACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTgcg                                     >  7:147591/1‑149 (MQ=255)
                                                                                                                 tGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTgcgc                                    >  7:35614/1‑147 (MQ=255)
                                                                                                                 tGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTgcgc                                    >  8:292921/1‑147 (MQ=255)
                                                                                                                   gATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCTGAGCGCGGGCGCTGTTCGGGCTGATCCTGGATTAATGGCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCca                                  <  8:353880/147‑1 (MQ=255)
                                                                                                                     tCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCcagcag                              >  6:342879/1‑149 (MQ=255)
                                                                                                                         gTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGTGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCcagcag                              <  1:177446/145‑1 (MQ=255)
                                                                                                                             cgcgCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACACACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGc                      <  7:224878/149‑1 (MQ=255)
                                                                                                                              gcgcTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTAc                                                                             >  4:348562/1‑93 (MQ=255)
                                                                                                                              gcgcTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTAc                                                                             <  3:348562/93‑1 (MQ=255)
                                                                                                                               cgcTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAg                    <  4:273523/149‑1 (MQ=255)
                                                                                                                                 cTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAggg                  <  8:317017/149‑1 (MQ=255)
                                                                                                                                  tGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAggg                  >  5:9877/1‑148 (MQ=255)
                                                                                                                                      gATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATc             >  4:183032/1‑149 (MQ=255)
                                                                                                                                      gATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATc             >  4:173991/1‑149 (MQ=255)
                                                                                                                                          ccAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTg         >  3:89941/1‑149 (MQ=255)
                                                                                                                                           cAGCGGCAGCGCAGCGCGGGCGCTGTTCGCGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTgg        <  5:334482/149‑1 (MQ=255)
                                                                                                                                                 cagcgcagcgCGGGCGCTGTTCGGGCTGATCCTGGCAGAATGCCGGCGAGATGCGGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTgggcggtc  <  2:275003/149‑3 (MQ=255)
                                                                                                                                                 cagcgcagcgCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTgggcgggc  >  3:161815/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCAGGCCATGCAGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGC  >  NC_002947/4136073‑4136365

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: