Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,258,778 A→C 5.1% intergenic (‑107/+90) 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,7780AC 5.1% 886.4 / 25.9 254intergenic (‑107/+90)rocE/PP_3728amino acid permease RocE/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (111/130);  new base C (5/8);  total (116/138)
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.58e-01

Only 100 of 263 total aligned reads displayed.

TCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACGG  >  NC_002947/4258629‑4258927
                                                                                                                                                     |                                                                                                                                                     
 cAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCac                                                                                                                                                      <  2:2157624/150‑2 (MQ=255)
  aGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAac                                                                                                                                                     >  1:1051930/1‑150 (MQ=255)
    ggCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAacac                                                                                                                                                   >  1:4316639/1‑150 (MQ=255)
       gcgGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGt                                                                                                                                                <  1:5409203/150‑1 (MQ=255)
        cgGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGtt                                                                                                                                               <  1:4650688/150‑1 (MQ=255)
        cgGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGtt                                                                                                                                               <  1:295941/150‑1 (MQ=255)
         gggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtcc                                                                                                                                              <  2:4113277/150‑3 (MQ=255)
         gggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTg                                                                                                                                              <  2:525499/150‑1 (MQ=255)
          ggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtcca                                                                                                                                             <  1:2245675/150‑4 (MQ=255)
           gTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTTGCACACGGtccat                                                                                                                                            <  1:4075079/150‑5 (MQ=255)
               gCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATTAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATgg                                                                                                                                        >  2:4626934/1‑150 (MQ=255)
                 tGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAt                                                                                                                                      <  2:4979712‑M2/150‑2 (MQ=255)
                   gATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAtta                                                                                                                                    <  2:4650299‑M2/150‑4 (MQ=255)
                     ttttttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttatt                                                                                                                                  >  2:3168661‑M2/1‑145 (MQ=255)
                        tttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattgga                                                                                                                               <  1:5355775‑M2/150‑9 (MQ=255)
                          tCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGaaa                                                                                                                             >  2:875890/1‑150 (MQ=255)
                           cAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAAc                                                                                                                            <  1:4430437/150‑1 (MQ=255)
                            aaCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACg                                                                                                                           >  1:4899632/1‑150 (MQ=255)
                             aCGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACgg                                                                                                                          <  2:2992423/150‑1 (MQ=255)
                              cGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGa                                                                                                                         <  1:1275441/150‑1 (MQ=255)
                              cGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGa                                                                                                                         <  1:4129673/150‑1 (MQ=255)
                               gCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                        <  1:2303975/150‑1 (MQ=255)
                               gCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                        <  2:2903938/150‑1 (MQ=255)
                                 tGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGt                                                                                                                      >  2:1591297/1‑150 (MQ=255)
                                   tcagtcaTGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTaa                                                                                                                    <  1:3829782/150‑1 (MQ=255)
                                     agtcaTGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAgg                                                                                                                  <  2:4498089/150‑1 (MQ=255)
                                     agtcaTGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAgg                                                                                                                  <  1:2314074/150‑1 (MQ=255)
                                       tcaTGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGgcg                                                                                                                >  2:170114/1‑150 (MQ=255)
                                           ggCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGt                                                                                                            <  2:3331163/150‑1 (MQ=255)
                                            gCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGtt                                                                                                           >  1:2888014/1‑150 (MQ=255)
                                              acacCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTtat                                                                                                         <  1:3565755/150‑1 (MQ=255)
                                               cacCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTtata                                                                                                        <  2:4302742/150‑1 (MQ=255)
                                                  cTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGa                                                                                                     >  1:2438658/1‑150 (MQ=255)
                                                   tCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGaa                                                                                                    >  1:5108852/1‑150 (MQ=255)
                                                   tCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAGTATTGGAAACGGATGTAAGGCGCCGTTATACGaa                                                                                                    >  2:2386121/1‑150 (MQ=255)
                                                    cGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAAc                                                                                                   <  1:1777289/150‑1 (MQ=255)
                                                     gCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACg                                                                                                  >  2:2095700/1‑150 (MQ=255)
                                                     gCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACg                                                                                                  >  2:4233164/1‑150 (MQ=255)
                                                           gCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGa                                                                                            <  1:4967882/150‑1 (MQ=255)
                                                            cGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGaa                                                                                           <  1:4796408/150‑1 (MQ=255)
                                                            cGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGaa                                                                                           <  2:335107/150‑1 (MQ=255)
                                                             ggCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAg                                                                                          >  1:1033458/1‑150 (MQ=255)
                                                               ccGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGAt                                                                                        >  2:2810281/1‑150 (MQ=255)
                                                                        gggAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGt                                                                               <  1:4738958/150‑1 (MQ=255)
                                                                             ggTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTaa                                                                          >  2:541900/1‑150 (MQ=255)
                                                                               ttCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCGATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGc                                                                        >  1:93050/1‑150 (MQ=255)
                                                                               ttCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGc                                                                        >  1:93205/1‑150 (MQ=255)
                                                                                tCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCt                                                                       >  2:5487236/1‑150 (MQ=255)
                                                                                tCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGAACGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCt                                                                       >  1:3398781/1‑150 (MQ=255)
                                                                                  gAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttt                                                                     >  1:4856088/1‑150 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                   >  1:3494519/1‑150 (MQ=255)
                                                                                    aTGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCttttt                                                                   <  2:1212129/150‑1 (MQ=255)
                                                                                        ggggCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTAc                                                               >  1:2558215/1‑150 (MQ=255)
                                                                                           gCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTc                                                            >  2:2995779/1‑150 (MQ=255)
                                                                                           gCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTc                                                            <  1:5089167/150‑1 (MQ=255)
                                                                                            cTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCa                                                           <  1:402696/150‑1 (MQ=255)
                                                                                            cTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCa                                                           >  2:2336608/1‑150 (MQ=255)
                                                                                            cTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCa                                                           <  1:2054250/150‑1 (MQ=255)
                                                                                             tGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAg                                                          >  1:2118152/1‑150 (MQ=255)
                                                                                             tGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAg                                                          >  2:3092344/1‑150 (MQ=255)
                                                                                               aGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGaaa                                                                                                                             >  2:1929739/1‑81 (MQ=255)
                                                                                               aGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGct                                                        <  1:378853/150‑1 (MQ=255)
                                                                                                  gcTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTg                                                     <  2:187462/150‑1 (MQ=255)
                                                                                                   cTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGa                                                    >  1:1548730/1‑150 (MQ=255)
                                                                                                    tACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGaa                                                   <  2:3483716/150‑1 (MQ=255)
                                                                                                    tACGGCATGTCAGTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTACGGCGCCGGTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTGTGTACGTCAGCTCTGaa                                                   >  2:5074412/1‑150 (MQ=255)
                                                                                                         cATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGc                                              <  1:4014683/150‑1 (MQ=255)
                                                                                                          aTGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCt                                             <  2:800520/150‑1 (MQ=255)
                                                                                                          aTGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTGTGTACGTCAGCTCTGAATCAGCt                                             >  2:3880828/1‑150 (MQ=255)
                                                                                                             tCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGAt                                                                                                                        <  1:3346562/72‑1 (MQ=255)
                                                                                                              cATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAACCGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTACACAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTc                                         <  2:3589199/150‑1 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                        >  2:2257075/1‑150 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                        >  2:685048/1‑150 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                        <  2:4856088/150‑1 (MQ=255)
                                                                                                               attattCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTcc                                        >  2:2582104/1‑150 (MQ=255)
                                                                                                                  attCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGc                                     <  2:1095328/150‑1 (MQ=255)
                                                                                                                     cTTGTCTCCAATAAGCCTCATTGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTg                                  <  2:4922310/150‑1 (MQ=255)
                                                                                                                      ttGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattggagacaagaataatgacatgccgtagcgctcagcccctcattcgaaccctcccacgcaccggccgcaatggcga                                                       <  2:1752614‑M2/128‑81 (MQ=255)
                                                                                                                      ttGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTgg                                 >  1:5588313/1‑150 (MQ=255)
                                                                                                                      ttGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTgg                                 <  2:2664330/150‑1 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                            <  1:5021475/150‑1 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                            >  1:78734/1‑150 (MQ=255)
                                                                                                                           tcCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTg                            >  2:1147447/1‑150 (MQ=255)
                                                                                                                            ccAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTgg                           >  2:4404418/1‑150 (MQ=255)
                                                                                                                            ccAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTgg                           >  1:2705575/1‑150 (MQ=255)
                                                                                                                             cAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGAttattggagacaagaataatgacatgccgtagcgctcagcccctcattcgaaccctcccacgcaccggccgcaatggcgaggtgtgccatgactgacagcgttgaaaaa                          <  2:730966‑M2/150‑110 (MQ=255)
                                                                                                                               aTAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCt                        <  2:3665645/150‑1 (MQ=255)
                                                                                                                                     catggcaaCCGTGTTGCACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTa                  >  1:4618414/8‑150 (MQ=255)
                                                                                                                                      tCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAg                 >  2:2816670/1‑150 (MQ=255)
                                                                                                                                      tCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAg                 >  1:1423528/1‑150 (MQ=255)
                                                                                                                                      tCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAg                 >  2:2751547/1‑150 (MQ=255)
                                                                                                                                         gcaaCCGTGTTGCACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTg              >  1:3197201/4‑150 (MQ=255)
                                                                                                                                          caaCCGTGTTGCACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTgg             >  1:1648673/3‑150 (MQ=255)
                                                                                                                                             ccGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTg          >  2:1326446/1‑150 (MQ=255)
                                                                                                                                                  tgCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCaa     >  2:3059462/1‑150 (MQ=255)
                                                                                                                                                    cAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACg   >  2:5177160/1‑150 (MQ=255)
                                                                                                                                                    cAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACg   <  1:1391291/150‑1 (MQ=255)
                                                                                                                                                     aaCACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACgg  >  1:692073/1‑150 (MQ=255)
                                                                                                                                                     aaCACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACgg  >  1:862202/1‑150 (MQ=255)
                                                                                                                                                     aaCACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACgg  >  1:278587/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTTGCCATGGATTATTGGAAACGGATGTAAGGCGCCGTTATACGAACGGATCGAAGATGTAAAAAGTCGTAAGCTTTTTGTACGTCAGCTCTGAATCAGCTGTTCCAGCGTGGCACTGGCCTGACCTAGCTGGGTGGCCAACGG  >  NC_002947/4258629‑4258927

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: