Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,559,926 C→G 5.9% R58R (CGG→CGC)  PP_3143 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,559,9260CG 5.9% 223.5 / 3.8 68R58R (CGG→CGC) PP_3143LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATGCAG  >  NC_002947/3559786‑3560063
                                                                                                                                            |                                                                                                                                         
gacgGGATCGAACGCACACGCAATGTGACGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCGCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAg                                                                                                                                   <  5:171686/149‑1 (MQ=255)
gacgGGAGCGAACGCCCACGCAAGTTGCCGCGGAATTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACAGCCGTCGGCCTCGTCCACCGTCTGGGCCAGGCGTAGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAg                                                                                                                                   <  6:268434/149‑1 (MQ=255)
 acgGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAg                                                                                                                                   >  3:374974/1‑148 (MQ=255)
  cgGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACTGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCa                                                                                                                                    <  3:175641/146‑1 (MQ=255)
  cgGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACTGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCa                                                                                                                                    >  4:175641/1‑146 (MQ=255)
         gAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCg                                                                                                                          <  7:196120/149‑1 (MQ=255)
          aaCGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCgg                                                                                                                         <  7:349566/149‑1 (MQ=255)
          aaCGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCgg                                                                                                                         <  7:251601/149‑1 (MQ=255)
               acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGcc                                                                                                                    >  8:172269/1‑149 (MQ=255)
               acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGcc                                                                                                                    >  4:110885/1‑149 (MQ=255)
               acacGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGcc                                                                                                                    >  4:105535/1‑149 (MQ=255)
                   gCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAgg                                                                                                                >  2:33164/1‑149 (MQ=255)
                        gtgCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGctc                                                                                                           >  2:220423/1‑149 (MQ=255)
                         tgCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGctc                                                                                                           <  1:369958/148‑1 (MQ=255)
                                gATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGtt                                                                                                   >  6:38953/1‑149 (MQ=255)
                                 aTTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTg                                                                                                  <  7:172269/149‑1 (MQ=255)
                                    tCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGgc                                                                                                <  3:105535/148‑1 (MQ=255)
                                    tCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGgc                                                                                                <  3:110885/148‑1 (MQ=255)
                                      ccccGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGgcgc                                                                                              <  3:172192/148‑1 (MQ=255)
                                      ccccGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCa                                                                                             >  6:72437/1‑149 (MQ=255)
                                        ccGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCCGGCGCTCTTCGAGTTGGCGCATc                                                                                           >  4:98432/1‑149 (MQ=255)
                                               cAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTgctggc                                                                                    >  4:53327/1‑149 (MQ=255)
                                                 gCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTgctggct                                                                                   >  5:187341/1‑148 (MQ=255)
                                                  ccGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCa                                                                                                                                    >  2:335043/1‑98 (MQ=255)
                                                  ccGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCa                                                                                                                                    <  1:335043/98‑1 (MQ=255)
                                                    gcagcaGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGAcc                                                                               >  5:228707/1‑149 (MQ=255)
                                                           cTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCg                                                                        >  3:324893/1‑149 (MQ=255)
                                                              ggCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCACTTGGCGCATCTGCTGGCTGACCGCGCCCgg                                                                       <  1:33164/147‑1 (MQ=255)
                                                               gCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCTGGCGCTCTTCCAGTTGGCGCATctgctg                                                                                      <  7:118107/131‑1 (MQ=255)
                                                               gCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATctgctg                                                                                      >  8:118107/1‑131 (MQ=255)
                                                                cTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGt                                                                    >  2:250391/1‑148 (MQ=255)
                                                                 ttCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGa                                                                  >  8:265901/1‑149 (MQ=255)
                                                                 ttCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGa                                                                  >  8:319285/1‑149 (MQ=255)
                                                                    gATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGAtg                                                                <  3:53327/148‑1 (MQ=255)
                                                                     aTACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGAtgtg                                                              <  7:265901/149‑1 (MQ=255)
                                                                     aTACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGAtgtg                                                              <  4:374974/149‑1 (MQ=255)
                                                                        cTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCa                                                            <  5:72437/148‑1 (MQ=255)
                                                                            cGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAgg                                                                                                                                  <  1:319890/74‑1 (MQ=255)
                                                                            cGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAgg                                                                                                                                  >  2:319890/1‑74 (MQ=255)
                                                                                ggCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCgg                                                   <  5:38953/149‑1 (MQ=255)
                                                                                    tCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGca                                                 >  6:269238/1‑147 (MQ=255)
                                                                                     cGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAg                                              >  2:255534/1‑149 (MQ=255)
                                                                                         cACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATctgctg                                                                                      <  2:343908/105‑1 (MQ=255)
                                                                                         cACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATctgctg                                                                                      >  1:343908/1‑105 (MQ=255)
                                                                                          aCCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGt                                         <  5:80154/149‑1 (MQ=255)
                                                                                           ccGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTg                                        <  7:319285/149‑1 (MQ=255)
                                                                                                 gggCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGa                                   >  4:332047/1‑148 (MQ=255)
                                                                                                      aGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTg                             >  4:321806/1‑149 (MQ=255)
                                                                                                         cGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCg                          <  3:321806/149‑1 (MQ=255)
                                                                                                           ttGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCggg                        >  1:48444/1‑149 (MQ=255)
                                                                                                              ccccTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTg                     >  8:4035/1‑149 (MQ=255)
                                                                                                              ccccTCGGCGCGCAGCTCCACTCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCt                      <  8:62387/148‑1 (MQ=255)
                                                                                                               cccTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgc                    >  1:93264/1‑149 (MQ=255)
                                                                                                                 cTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGAt                                                                 >  6:302589/1‑102 (MQ=255)
                                                                                                                 cTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGAt                                                                 <  5:302589/102‑1 (MQ=255)
                                                                                                                 cTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgcgc                  <  6:228707/149‑1 (MQ=255)
                                                                                                                   cGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAg                >  5:122788/1‑149 (MQ=255)
                                                                                                                    ggCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCa                 <  3:332047/147‑1 (MQ=255)
                                                                                                                        cTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCaaa           <  1:255534/149‑1 (MQ=255)
                                                                                                                             cTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACAt      >  5:14273/1‑149 (MQ=255)
                                                                                                                                cACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATg     >  4:226379/1‑147 (MQ=255)
                                                                                                                                cACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATGc    <  7:4035/148‑1 (MQ=255)
                                                                                                                                cACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATGc    <  2:93264/148‑1 (MQ=255)
                                                                                                                                 aCGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATGCAg  >  2:172639/1‑149 (MQ=255)
                                                                                                                                   gccggGCGCGCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgc                    <  4:116839/129‑1 (MQ=255)
                                                                                                                                   gccggGCGCGCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgc                    >  3:116839/1‑129 (MQ=255)
                                                                                                                                   gccggGCGCGCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgc                    >  3:111436/1‑129 (MQ=255)
                                                                                                                                   gccggGCGCGCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCGGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTgc                    <  4:111436/129‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GACGGGATCGAACGCACACGCAATGTGCCGCGGATTTCCCCCGCGTCCAGCCGCAGCAGCTCGGCTTCGATACTGCCGTAGGCCTCGGCCACCGTCTGGGCCAGGCGTTGCCCCTCGGCGCTCAGCTCCACGCCGCGCGCCCGGCGCAGGAACAGCCGGTAGCCCAGGCGCTCTTCCAGTTGGCGCATCTGCTGGCTGACCGCGCCCGGGGTGATGTGCAGCTCTTCGGCACAGCGGGTGAAGGACAGGTGCCGGGCTGCGCAGGCAAATACATGCAG  >  NC_002947/3559786‑3560063

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