Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,136,228 G→T 8.0% R146L (CGG→CTG)  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,2280GT 8.0% 226.5 / 8.8 75R146L (CGG→CTG) PP_3640AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base T (2/4);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

AGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTGGTCGCCC  >  NC_002947/4136084‑4136374
                                                                                                                                                 |                                                                                                                                                  
aGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGgcg                                                                                                                                                <  3:26426/149‑1 (MQ=255)
  gcgcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGTATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCt                                                                                                                                              <  6:243649/149‑1 (MQ=255)
    gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGt                                                                                                                                            >  3:168385/1‑149 (MQ=255)
    gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGt                                                                                                                                            >  3:183188/1‑149 (MQ=255)
    gcTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGGCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTgg                                                                                                                                            <  7:97611/149‑2 (MQ=255)
     cTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGtt                                                                                                                                           >  7:30636/1‑149 (MQ=255)
        gATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCgg                                                                                                                                        >  8:224878/1‑149 (MQ=255)
            gCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCt                                                                                                                                     <  8:38464/148‑1 (MQ=255)
             ccAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCAGGACGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCTGGCt                                                                                                                                     >  1:275003/1‑148 (MQ=255)
             ccAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGa                                                                                                                                   >  7:277868/1‑149 (MQ=255)
              cAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTggga                                                                                                                                          >  1:335696/1‑138 (MQ=255)
              cAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTggga                                                                                                                                          <  2:335696/141‑4 (MQ=255)
               aGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATc                                                                                                                                 >  6:174969/1‑149 (MQ=255)
                  ccgccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTg                                                                                                                              <  5:41918/149‑1 (MQ=255)
                   cgtcgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTgg                                                                                                                             <  5:289494/146‑1 (MQ=255)
                    gtcgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGga                                                                                                                            <  7:287802/147‑1 (MQ=255)
                    gccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGga                                                                                                                            <  5:307119/149‑1 (MQ=255)
                     ccgACGGCGTTGGCTTGGTCTGCGGGTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGTGGTTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTTTTCGGGCTGATCCTGgaa                                                                                                                           <  5:247714/149‑1 (MQ=255)
                     ccgACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGAATGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaa                                                                                                                           <  1:374447/149‑1 (MQ=255)
                       gACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaaga                                                                                                                         <  8:30636/149‑1 (MQ=255)
                       gACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGgaag                                                                                                                          >  6:296494/1‑148 (MQ=255)
                          ggCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATg                                                                                                                      <  4:210677/149‑1 (MQ=255)
                                 gCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGc                                                                                                                 >  4:211646/1‑147 (MQ=255)
                                 gCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCga                                                                                                               <  8:103395/149‑1 (MQ=255)
                                      gTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGc                                                                                                          <  5:296494/149‑1 (MQ=255)
                                       tCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGAGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCt                                                                                                         >  2:300253/1‑149 (MQ=255)
                                         tGCCGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTCAAGATCCAAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                                       >  6:49305/1‑149 (MQ=255)
                                           cGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                                       <  5:83064/147‑1 (MQ=255)
                                           cGGCTTCTTCCA‑TTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTgg                                                                                                       <  2:68413/146‑1 (MQ=255)
                                              cttcttCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCAt                                                                                                  <  4:168385/149‑1 (MQ=255)
                                               ttcttcCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCAt                                                                                                  >  8:222375/1‑148 (MQ=255)
                                                 ctgccATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATca                                                                                                <  5:42780/145‑1 (MQ=255)
                                                    ccATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGTCATCACAg                                                                                             <  3:131251/148‑1 (MQ=255)
                                                       ttttCA‑GTCGGGCCTCTGTTCGCTGATTGTCGAGGGTCGCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGc                                                                                         <  4:107840/149‑1 (MQ=255)
                                                          tCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGcc                                                                                                    <  6:299637/135‑1 (MQ=255)
                                                          tCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGcc                                                                                                    >  5:299637/1‑135 (MQ=255)
                                                           cA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTg                                                                                                                                             >  2:35705/1‑93 (MQ=38)
                                                           cA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGCCCGCGCTGCGCTGCCGCTg                                                                                                                                             <  1:35705/93‑1 (MQ=38)
                                                            a‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACt                                                                                     <  4:183188/148‑1 (MQ=255)
                                                                 gggCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAAc                                                                                <  6:106536/149‑1 (MQ=255)
                                                                 gggCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAAc                                                                                <  5:174969/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)
                                                                                                                                        gcgcagcgcGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTg         >  4:72479/1‑149 (MQ=255)
                                                                                                                                         cgcagcgcGGGCGCTGTTCGGGCTGATGCTGGAAGAATGCCGGCGAGATGCTGGGACATCACTGATGCTTCTGGAACGGGTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGGATGTGCACGCAGGGCAATCGCTGGGCGGGCTg         <  5:49305/148‑1 (MQ=255)
                                                                                                                                           cagcgcGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTGGtc      <  2:249042/149‑1 (MQ=255)
                                                                                                                                             gcgcGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTGGtcgc    <  6:9877/149‑1 (MQ=255)
                                                                                                                                               gcgGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTGGtcgccc  <  3:173991/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGGCGCTGGATAGCCAGGCCGCCGACGGCGTTGGCTTGGTCTGCGGCTTCTTCCATTTTCA‑GTCGGGCCTCTGCTCGCTGATTGTCGAGGGTCTCGCTGACTGGATCTTGTTACGCGCTGAAGATCCCAGCGGCAGCGCAGCGCGGGCGCTGTTCGGGCTGATCCTGGAAGAATGCCGGCGAGATGCTGGGCCATCACAGATGCTACTGGAACGGCTGACGCACTTGTTGTTCCTGTACGTGCTGCGCCAGCAGGTGCACGCAGGGCAATCGCTGGGCGGGCTGGTCGCCC  >  NC_002947/4136084‑4136374

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