Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,258,777 C→G 5.1% intergenic (‑106/+91) rocE ← / ← PP_3728 amino acid permease RocE/two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,258,7770CG 5.1% 874.2 / 26.5 255intergenic (‑106/+91)rocE/PP_3728amino acid permease RocE/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (110/132);  new base G (5/8);  total (115/140)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 258 total aligned reads displayed.

TCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACG  >  NC_002947/4258629‑4258926
                                                                                                                                                    |                                                                                                                                                     
 cAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCaa                                                                                                                                                     <  1:3315359/150‑1 (MQ=255)
 cAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCaa                                                                                                                                                     <  1:3544940/150‑1 (MQ=255)
  aGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAac                                                                                                                                                    >  1:1051930/1‑150 (MQ=255)
  aGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAac                                                                                                                                                    >  1:4372991/1‑150 (MQ=255)
  aGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAac                                                                                                                                                    >  1:693193/1‑150 (MQ=255)
        cgGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGtt                                                                                                                                              <  1:295941/150‑1 (MQ=255)
         gggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtcc                                                                                                                                             <  2:4113277/150‑3 (MQ=255)
         gggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTg                                                                                                                                             <  2:525499/150‑1 (MQ=255)
          ggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGc                                                                                                                                            <  2:1276825/150‑1 (MQ=255)
               gCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATTAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATgg                                                                                                                                       >  2:4626934/1‑150 (MQ=255)
                 tGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAt                                                                                                                                     <  1:2754150‑M2/150‑2 (MQ=255)
                   gATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAtta                                                                                                                                   >  2:1942416‑M2/1‑147 (MQ=255)
                     ttttttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttatt                                                                                                                                 >  2:1722162‑M2/1‑145 (MQ=255)
                     ttttgtCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttatt                                                                                                                                 >  2:1847372‑M2/1‑145 (MQ=255)
                      tttttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtccatga                                                                                                                                         <  2:2276055/141‑7 (MQ=255)
                        tttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattgga                                                                                                                              <  1:5355775‑M2/150‑9 (MQ=255)
                          tCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGaaa                                                                                                                            >  2:875890/1‑150 (MQ=255)
                           cAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAAc                                                                                                                           <  1:4430437/150‑1 (MQ=255)
                           cAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAAc                                                                                                                           <  2:308765/150‑1 (MQ=255)
                            aaCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACg                                                                                                                          >  1:4899632/1‑150 (MQ=255)
                            aaCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACg                                                                                                                          >  2:2143071/1‑150 (MQ=255)
                            aaCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACg                                                                                                                          >  2:3683655/1‑150 (MQ=255)
                             aCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACgg                                                                                                                         <  2:2992275/150‑1 (MQ=255)
                              cGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGa                                                                                                                        <  1:1275441/150‑1 (MQ=255)
                              cGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGa                                                                                                                        <  1:4129673/150‑1 (MQ=255)
                               gCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       <  1:1968765/150‑1 (MQ=255)
                               gCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       <  1:2303975/150‑1 (MQ=255)
                               gCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       >  1:5005444/1‑150 (MQ=255)
                               gCTGTCAGTCAGGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       <  2:4663289/150‑1 (MQ=255)
                                 tGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGGGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCAGGTCAGTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCAGGGAGTATTGGAAACGGATGt                                                                                                                     >  1:823400/1‑150 (MQ=255)
                                     agtcaTGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAgg                                                                                                                 >  2:3618049/1‑150 (MQ=255)
                                           ggCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGt                                                                                                           >  2:4873765/1‑150 (MQ=255)
                                            gCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGtt                                                                                                          >  1:2888014/1‑150 (MQ=255)
                                                  cTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGa                                                                                                    >  2:2618150/1‑150 (MQ=255)
                                                  cTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGa                                                                                                    >  2:1037277/1‑150 (MQ=255)
                                                   tCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGaa                                                                                                   >  1:5108852/1‑150 (MQ=255)
                                                    cGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAAc                                                                                                  <  1:1777289/150‑1 (MQ=255)
                                                            cGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGaa                                                                                          <  1:4796408/150‑1 (MQ=255)
                                                            cGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGaa                                                                                          <  2:266546/150‑1 (MQ=255)
                                                             ggCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAg                                                                                         >  1:3866703/1‑150 (MQ=255)
                                                             ggCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAg                                                                                         >  1:1033458/1‑150 (MQ=255)
                                                              gCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGa                                                                                        <  2:4237101/150‑1 (MQ=255)
                                                               ccGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGAt                                                                                       <  1:4676796/150‑1 (MQ=255)
                                                               ccGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGAt                                                                                       >  2:2810281/1‑150 (MQ=255)
                                                                   tGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTaa                                                                                   <  2:135849/150‑1 (MQ=255)
                                                                        gggAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGt                                                                              <  1:4738958/150‑1 (MQ=255)
                                                                         ggAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTc                                                                             <  2:4996861/150‑1 (MQ=255)
                                                                          gAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtcc                                                                                                                                             >  1:4098012/1‑83 (MQ=255)
                                                                          gAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCg                                                                            <  2:4539065/150‑1 (MQ=255)
                                                                             ggTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTaa                                                                         >  2:541900/1‑150 (MQ=255)
                                                                                tCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGAACGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCt                                                                      >  1:3398781/1‑150 (MQ=255)
                                                                                  gAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttt                                                                    >  1:932044/1‑150 (MQ=255)
                                                                                  gAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttt                                                                    >  1:4856088/1‑150 (MQ=255)
                                                                                  gAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttt                                                                    <  2:5254106/150‑1 (MQ=255)
                                                                                   aaTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCtttt                                                                   >  2:3435072/1‑150 (MQ=255)
                                                                                    aTGAGGGGCTGATCGCTTCGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGTATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  1:2203896/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  2:1120653/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  2:2396652/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  2:1901825/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  2:1165787/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCatttt                                                                  <  2:2228729/150‑1 (MQ=255)
                                                                                    aTGAGGGGCTGAGAGCTACGTCATGTCATTATTCTTGTCTCCAATAAGCCTCATTGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                  <  1:2386768/150‑1 (MQ=255)
                                                                                        ggggCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTAc                                                              >  1:2558215/1‑150 (MQ=255)
                                                                                          ggCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGt                                                            <  2:5545502/150‑1 (MQ=255)
                                                                                           gCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTc                                                           <  1:5089167/150‑1 (MQ=255)
                                                                                           gCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTc                                                           >  2:2995457/1‑150 (MQ=255)
                                                                                            cTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCa                                                          >  2:2336608/1‑150 (MQ=255)
                                                                                            cTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCa                                                          <  1:402696/150‑1 (MQ=255)
                                                                                             tGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAg                                                         >  1:2118152/1‑150 (MQ=255)
                                                                                               aGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGct                                                       <  1:378853/150‑1 (MQ=255)
                                                                                                 cgcTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGctct                                                     <  2:1621228/150‑1 (MQ=255)
                                                                                                  gcTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTg                                                    <  2:4338225/150‑1 (MQ=255)
                                                                                                    tACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCTTAAGCTTTTTGTACGTCAGCTCTGaa                                                  >  2:3008182/1‑150 (MQ=255)
                                                                                                    tACGGCATGTCAGTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTACGGCGCCGGTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTGTGTACGTCAGCTCTGaa                                                  >  2:5074412/1‑150 (MQ=255)
                                                                                                        gCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAg                                              <  2:3104453/150‑1 (MQ=255)
                                                                                                         cATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGc                                             <  1:4014798/150‑1 (MQ=255)
                                                                                                          aTGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCt                                            >  2:3959555/1‑150 (MQ=255)
                                                                                                          aTGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCt                                            <  2:800520/150‑1 (MQ=255)
                                                                                                             tCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       >  2:3346562/1‑72 (MQ=255)
                                                                                                             tCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       <  1:3346562/72‑1 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattggagacaagaataatgacatgccgtagcgctcagcccctcattcgaaccctcccacgcaccggccgcaatggcgaggtgtgccatgactg                                       <  1:919373‑M2/150‑96 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                       <  2:2828828/70‑1 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                       >  2:685048/1‑150 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                       >  2:4405535/1‑150 (MQ=255)
                                                                                                                     cTTGTCTCCAATAAGCCTCATTGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTg                                 <  2:4922310/150‑1 (MQ=255)
                                                                                                                      ttGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTgg                                >  1:269498/1‑150 (MQ=255)
                                                                                                                      ttGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTgg                                >  1:5588313/1‑150 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                           >  2:3191796/1‑150 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                           >  1:1469616/1‑150 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                           <  1:5021475/150‑1 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                           >  1:2308831/1‑150 (MQ=255)
                                                                                                                            ccAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTgg                          <  1:3634192/150‑1 (MQ=255)
                                                                                                                             ctctacgCCTCATGGACCGTGTGCAAGAGGGTTGCCTTGGATTTTTGGCAACGGATGGAAGGCGCCGTTATACGACCGGATCGAAGATGTAAACAGTCGTAAGCTTTTTGTACGTCAGCTCTGACTCTGCTGTTCGAGCGTGGCACTGGc                         <  2:3495864/144‑1 (MQ=255)
                                                                                                                             cAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattggagacaagaataatgacatgccgtagcgctcagcccctcattcgaaccctcccacgcaccggccgcaatggcgaggtgtgccatgactgacagcgttgaaaaa                         <  2:730966‑M2/150‑110 (MQ=255)
                                                                                                                                 aaGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGa                     <  2:2787635/150‑1 (MQ=255)
                                                                                                                                     catggcaaCCGTGTTGCACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTa                 >  1:4618414/8‑150 (MQ=255)
                                                                                                                                      tCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAg                >  1:1423528/1‑150 (MQ=255)
                                                                                                                                          caaCCGTGTTGCACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTgg            >  1:1648673/3‑150 (MQ=255)
                                                                                                                                             ccGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTg         <  1:540139/150‑1 (MQ=255)
                                                                                                                                                    cAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACg  <  1:1391291/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACG  >  NC_002947/4258629‑4258926

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