Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,427,829 T→G 13.0% E33D (GAA→GAC) ‡ PP_2126 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,427,8290TG13.0% 498.8 / 21.8 199E33D (GAA→GAC) ‡PP_2126LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/106);  new base G (29/2);  total (91/108)
Fisher's exact test for biased strand distribution p-value = 3.10e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GACAGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACA  >  NC_002947/2427680‑2427978
                                                                                                                                                     |                                                                                                                                                     
 aCAGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCtt                                                                                                                                                      >  2:2021750/1‑150 (MQ=255)
  cAGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTc                                                                                                                                                     <  1:3055328/150‑1 (MQ=255)
   aGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCg                                                                                                                                                    <  2:724913/150‑1 (MQ=255)
   aGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCg                                                                                                                                                    >  2:4016186/1‑150 (MQ=255)
     gATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGcc                                                                                                                                                  <  1:3288129/150‑1 (MQ=255)
     gATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGcc                                                                                                                                                  <  1:3770369/150‑1 (MQ=255)
       tCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGa                                                                                                                                                >  1:3292817/1‑150 (MQ=255)
        cAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGAc                                                                                                                                               <  2:120796/150‑1 (MQ=255)
        cAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGAc                                                                                                                                               <  1:2446125/150‑1 (MQ=255)
          gCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTACGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGAAGCCGAcga                                                                                                                                             <  1:1776844/150‑3 (MQ=255)
          gCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCa                                                                                                                                             >  1:2997426/1‑150 (MQ=255)
           cACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGACCAc                                                                                                                                            >  2:624861/1‑150 (MQ=255)
              tttGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCCACCCCGTc                                                                                                                                         >  2:1362914/1‑150 (MQ=255)
              tttGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCATAGCCCAacaccac                                                                                                                                         >  1:3810185/1‑147 (MQ=255)
                 gCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTGGCCGCCCACCTCATa                                                                                                                                      >  1:1216053/1‑150 (MQ=255)
                  cAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAg                                                                                                                                     <  1:1022306/150‑1 (MQ=255)
                   aGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGc                                                                                                                                    >  2:368461/1‑150 (MQ=255)
                      tACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCgg                                                                                                                                 <  1:3739201/150‑1 (MQ=255)
                        cTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCAAAGCCcggg                                                                                                                               >  1:1088543/1‑150 (MQ=255)
                         tGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTAGCCGACCACCGAATAACCGGGGa                                                                                                                              >  1:2433987/1‑150 (MQ=255)
                          gCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGat                                                                                                                             >  2:3022374/1‑150 (MQ=255)
                             gCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTAGCCCACCACCACAAAGACGGGGatata                                                                                                                          >  1:4299139/1‑149 (MQ=255)
                                ggCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCa                                                                                                                       <  1:4399706/150‑1 (MQ=255)
                                ggCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCa                                                                                                                       <  1:3221105/150‑1 (MQ=255)
                                     gCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCa                                                                                                                  <  2:1076940/150‑1 (MQ=255)
                                       aGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCATAGCCCACCACCTCATAGCCGGGGATATAGCAGACCact                                                                                                                >  1:2571157/1‑148 (MQ=255)
                                           gggTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGccc                                                                                                            <  1:97387/150‑1 (MQ=255)
                                           gggTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGccc                                                                                                            >  2:2925050/1‑150 (MQ=255)
                                           gggTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGccc                                                                                                            <  1:2351920/150‑1 (MQ=255)
                                             gTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGGCATAGCCGGGGATATCGCAGACCAGTGCCCGc                                                                                                          >  2:872767/1‑150 (MQ=255)
                                              tAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCCGCATCGCCGACCAAGTAATAGCCCGGGATATCGCACACCAGTTCCCGCa                                                                                                         >  2:3137809/1‑150 (MQ=255)
                                               aaGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGAGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTAGCCGACCACGTCAAAGCCGGGGAAAACTCAGAACAATGCCCGCAc                                                                                                        >  2:1884508/1‑150 (MQ=255)
                                                 gggCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCAccc                                                                                                      <  2:1797598/150‑1 (MQ=255)
                                                  ggCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTAGCCGGACACGTCATATGCGGGGTTATAGCAGAACAGTTGCCGCAcccc                                                                                                     >  2:5495/1‑150 (MQ=255)
                                                  ggCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGCCCCCGGCATAGCCGGGGCTATCGCCGCCCACGGCCCGCAcccc                                                                                                     >  2:761331/1‑150 (MQ=255)
                                                   gCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGCCCCCGTCCGAGCCGGGGCTATCGCCGACCCGGGCCCGCCCCCgg                                                                                                   >  2:770193/1‑150 (MQ=255)
                                                       ccAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGgcgc                                                                                                <  2:3259496/150‑1 (MQ=255)
                                                       ccAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGgcgc                                                                                                <  1:509113/150‑1 (MQ=255)
                                                        cAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTGGGCGTGGCCGCCCACCGCAAATCCGGGGGTATAGCCGAACCCTGCCCCCACCCCcgcgcg                                                                                               >  1:1545324/1‑150 (MQ=255)
                                                           gCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGAt                                                                                            <  1:368718/150‑1 (MQ=255)
                                                           gCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGACCCCGGCAAATCCGGGGGTATAGCCGACAAGTGGCCGCACCCCCGCGCggga                                                                                            >  1:1893917/1‑148 (MQ=255)
                                                            ccGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATc                                                                                           <  2:4248578/150‑1 (MQ=255)
                                                            ccGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATc                                                                                           <  1:3022374/150‑1 (MQ=255)
                                                               cTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGACCCCCTGATAGGCGGGGGTATAGGCTACCAGTGCCCGCACCCCCGCGCGGATCAGc                                                                                        >  1:1746371/1‑150 (MQ=255)
                                                                 ggTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGa                                                                                      >  2:1814280/1‑150 (MQ=255)
                                                                 ggTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGa                                                                                      <  2:1689802/150‑1 (MQ=255)
                                                                     gatcGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTgg                                                                                  <  2:4483642/147‑1 (MQ=255)
                                                                     gAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGGACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTgg                                                                                  <  2:1509309/150‑1 (MQ=255)
                                                                       aCGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTc                                                                                <  2:3852069/150‑1 (MQ=255)
                                                                       aCGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTc                                                                                <  1:3137809/150‑1 (MQ=255)
                                                                           aTGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCt                                                                            <  2:2434943/150‑1 (MQ=255)
                                                                            tGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTa                                                                           >  2:4283837/1‑150 (MQ=255)
                                                                              gagaTGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTAcc                                                                         <  1:4441255/150‑1 (MQ=255)
                                                                               agaTGTCCAGCAAGATGATTTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACca                                                                        <  1:4160708/150‑1 (MQ=255)
                                                                                 aTGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACcagc                                                                      <  1:2925050/150‑1 (MQ=255)
                                                                                   gTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCATCGCCGAACACCTCAAATCCGGGGATATAGCAGAACAGTGCCCGCAACCCGGCCCGGATAAGCGAGTGGTAAAATACcagcag                                                                    >  2:3798185/1‑150 (MQ=255)
                                                                                    tCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTAGCCGACCACGTAATAGCCGGGGATATCGCAGACCAGTTCCCCCACCCCGGCGCGGATAAGCGAGTGGTCAACTACCAGCagt                                                                   >  1:4490304/1‑150 (MQ=255)
                                                                                    tCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCCGCATCGCCCACCACGTAAAAGCCGGGGATAAAGCAGACCAGTGCCCGCACCCCCGCGCGGAACAGCGAGTGGTCAACAACcaccagg                                                                   >  1:1071949/1‑149 (MQ=255)
                                                                                      cAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCagtag                                                                 >  2:4120597/1‑150 (MQ=255)
                                                                                      cAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCATCGCCGACCACCTCATAGCCGGGGCTATCGCAGACCAGTGCCCCCACCCCGGCGCCGATCAGCGAGTGGTCAACTAACACCagaag                                                                 >  2:1759282/1‑150 (MQ=255)
                                                                                       aGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTAGCCCGCCACCTCAAATCCGGGGATATCGCAGACCAGAGCCCGCACCCCCGCGCGGGTCAGCGAGTGGTCATCAACCAGCagttgt                                                                >  1:2386229/1‑150 (MQ=255)
                                                                                         cAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCt                                                              >  2:3645010/1‑150 (MQ=255)
                                                                                         cAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCCCACCCCGGCGCGGATCAGCGAGTGGTCATATAACAGCAGTAGTCt                                                              >  2:3510432/1‑150 (MQ=255)
                                                                                               gatgTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATATCCGGGGATATAGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGt                                                        >  2:2794536/1‑150 (MQ=255)
                                                                                                  gTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAt                                                     <  2:4344085/150‑1 (MQ=255)
                                                                                                   tCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAtg                                                    <  2:3726285/150‑1 (MQ=255)
                                                                                                      ggCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAtggtt                                                 <  1:272674/150‑1 (MQ=255)
                                                                                                        cgcCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAtggttgg                                               <  1:85201/150‑1 (MQ=255)
                                                                                                         gcCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAtggttggt                                              <  1:1778692/150‑1 (MQ=255)
                                                                                                           cAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTaa                                            <  2:583610/150‑1 (MQ=255)
                                                                                                            aGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCACACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCAGGGTTGGTAAt                                           >  2:3188948/1‑150 (MQ=255)
                                                                                                             gtgtGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATa                                          <  1:4016186/150‑1 (MQ=255)
                                                                                                                    tACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTgg                                   <  2:2081481/150‑1 (MQ=255)
                                                                                                                     aCCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCGACCACGTCCTAGCCGGGGATATGGCAGACCACTGCCCGCACCCCGGCGCGGAACAGCGAGTGGTCCTCTACCAGCAGTAGTGTACAGGTCCAGGGTGGTAAAACTCCTGgc                                  >  2:1146683/1‑150 (MQ=255)
                                                                                                                      ccTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTAGCCGACCACGTCATATCCGGGGATATCGCAGACCACTGCCCCCACCCCCGCGCGGATCAGCGAGTGGTCATCTACCACCACTAGTATACAGGTCATGGTTGGTAATAATCCAGgcg                                 >  2:472847/1‑150 (MQ=255)
                                                                                                                         ggTCGAGTAGCTGGTTGCCGTCGTCGGCGTCGCCCACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATAAACCTGGCGCGt                              >  2:4292611/1‑150 (MQ=255)
                                                                                                                              aGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTg                         <  1:2253002/150‑1 (MQ=255)
                                                                                                                               gTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTgg                        <  1:3575701/150‑1 (MQ=255)
                                                                                                                               gTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTgg                        <  2:1036792/150‑1 (MQ=255)
                                                                                                                                tAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTggg                       <  1:2059587/150‑1 (MQ=255)
                                                                                                                                 aGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGt                      <  1:3104194/150‑1 (MQ=255)
                                                                                                                                        tGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCTGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCtgg               <  2:3534050/150‑1 (MQ=255)
                                                                                                                                        tGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCggg               <  2:81766/150‑1 (MQ=255)
                                                                                                                                          acgtggtCGGCGAAGCCGACGACGGCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGa             >  1:1399621/2‑150 (MQ=255)
                                                                                                                                          acgtggtCGGCGAAGCCGACGACGGCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGa             >  1:815844/2‑150 (MQ=255)
                                                                                                                                            gtcgtcGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCCGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGAcc           >  2:1019760/1‑150 (MQ=255)
                                                                                                                                             tcgtcgGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTACTACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCa          >  2:2554103/1‑150 (MQ=255)
                                                                                                                                               gtcgGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACg        >  2:2500780/1‑150 (MQ=255)
                                                                                                                                                   gCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGaa    <  1:1271662/150‑1 (MQ=255)
                                                                                                                                                   gCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGaa    <  1:2689058/150‑1 (MQ=255)
                                                                                                                                                   gCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGaa    <  1:2816034/150‑1 (MQ=255)
                                                                                                                                                   gCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGaa    >  1:2856913/1‑150 (MQ=255)
                                                                                                                                                   gCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGaa    <  1:722159/150‑1 (MQ=255)
                                                                                                                                                    cTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAAc   <  1:2021750/150‑1 (MQ=255)
                                                                                                                                                    cTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAAc   >  1:1876926/1‑150 (MQ=255)
                                                                                                                                                     ttCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACa  >  2:460148/1‑150 (MQ=255)
                                                                                                                                                     ttCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACa  >  2:528126/1‑150 (MQ=255)
                                                                                                                                                     ttCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACa  >  2:3751411/1‑150 (MQ=255)
                                                                                                                                                     ttCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACa  >  2:2175884/1‑150 (MQ=255)
                                                                                                                                                     ttCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTCGGTGCGCGGGGACCACGGGAACa  >  2:3325447/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACAGGATCAGCACTTTGCAGGTACTGCCGCTGGCCCGCAGTTGGGTAAGGGCGTCCAGGCCGCTGGTCGAACGCATGGAGATGTCCAGCAAGATGATGTCGGGCGCCAGTGTGCGTACCTGGTCGAGTAGCTGGTTGCCGTCGTCGGCTTCGCCGACCACGTCATAGCCGGGGATATCGCAGACCAGTGCCCGCACCCCGGCGCGGATCAGCGAGTGGTCATCTACCAGCAGTAGTCTACAGGTCATGGTTGGTAATACTCCTGGCGCGTTGCTGGGTGCGCGGGGACCACGGGAACA  >  NC_002947/2427680‑2427978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: