Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,928 C→A 91.5% E903* (GAA→TAA)  torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,048,9280CA91.5% 815.5 / 69.7 282E903* (GAA→TAA) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base C (16/8);  new base A (146/112);  total (162/120)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 282 total aligned reads displayed.

GGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTT  >  CP009273/1048781‑1049074
                                                                                                                                                   |                                                                                                                                                  
ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACa                                                                                                                                                  <  1:2107229/150‑1 (MQ=255)
  ccgccgATGCCAGCGGGTTCTTTTATTATCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCTGGTCTTTCTTATGCAACCATGCTTACaga                                                                                                                                                <  1:3448443/150‑1 (MQ=255)
   cgccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                               <  2:3506806/150‑1 (MQ=255)
      cgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCgg                                                                                                                                            <  2:376946/150‑1 (MQ=255)
      cgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCgg                                                                                                                                            <  2:2689715/150‑1 (MQ=255)
        cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCa                                                                                                                                          >  1:2343252/1‑150 (MQ=255)
          gCCAGCGGGTTCTTTTATTATCTGCCATGATTTTGCCAGTCCGTTAAGATGTATACCAAATGCCACTATTATAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGACTTTCTTATTCAACCATGCTTACAGAGCGGCAac                                                                                                                                        >  2:2245331/1‑150 (MQ=255)
             aGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACt                                                                                                                                     >  1:1714913/1‑150 (MQ=255)
               cGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGc                                                                                                                                   <  1:1471529/150‑1 (MQ=255)
               cGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGc                                                                                                                                   >  1:2686410/1‑150 (MQ=255)
                  gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtg                                                                                                                                >  2:1346280/1‑150 (MQ=255)
                       tttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAtt                                                                                                                           >  1:3896994/1‑150 (MQ=255)
                         tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttc                                                                                                                         <  1:2442301/150‑1 (MQ=255)
                         tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttc                                                                                                                         >  2:4137404/1‑150 (MQ=255)
                           ttCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttctt                                                                                                                       <  2:2211848/150‑1 (MQ=255)
                            tctcTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                      <  2:2481368/150‑1 (MQ=255)
                            tctcTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                      <  2:1037622/150‑1 (MQ=255)
                            tctcGGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                      <  2:168834/150‑1 (MQ=255)
                             ctctGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCg                                                                                                                     <  2:1075856/150‑1 (MQ=255)
                               ctGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCgtg                                                                                                                   <  1:4473219/150‑1 (MQ=255)
                                tGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCgtgt                                                                                                                  >  2:4601331/1‑150 (MQ=255)
                                  ccATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCACATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTgg                                                                                                                <  2:1567674/150‑1 (MQ=255)
                                    aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                              <  1:2470060/150‑1 (MQ=255)
                                     tGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAg                                                                                                             >  1:342433/1‑150 (MQ=255)
                                     tGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAg                                                                                                             >  1:3231664/1‑150 (MQ=255)
                                       aTTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggg                                                                                                           >  1:4002566/1‑150 (MQ=255)
                                          ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggggg                                                                                                        <  2:510443/150‑1 (MQ=255)
                                                 tCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTTATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATaa                                                                                                 >  1:2293001/1‑150 (MQ=255)
                                                      tAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgct                                                                                            >  1:1150766/1‑150 (MQ=255)
                                                        aGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgc                                                                                          <  1:2906966/150‑1 (MQ=255)
                                                            gTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctgct                                                                                      >  1:1186223/1‑150 (MQ=255)
                                                              ataCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTcc                                                                                    >  1:1494108/1‑150 (MQ=255)
                                                               taCCAAATGCCACTAGTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                   <  1:636106/150‑1 (MQ=255)
                                                                aCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAg                                                                                  >  2:1813396/1‑150 (MQ=255)
                                                                  cAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGtt                                                                                <  1:1584384/150‑1 (MQ=255)
                                                                     aTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTtgtg                                                                             <  1:1787148/150‑1 (MQ=255)
                                                                      tGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCaagatgac                                                                            >  2:4594286/1‑143 (MQ=255)
                                                                      tGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                            <  1:4009950/150‑1 (MQ=255)
                                                                      tGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                            <  1:1685306/150‑1 (MQ=255)
                                                                       gCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                           >  1:2316150/1‑150 (MQ=255)
                                                                       gCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                           >  2:4429812/1‑150 (MQ=255)
                                                                          aCTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                        <  2:3346239/150‑1 (MQ=255)
                                                                          aCTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                        >  2:393373/1‑150 (MQ=255)
                                                                           cTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                                       <  1:1373135/150‑1 (MQ=255)
                                                                            tATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAg                                                                      <  2:2622301/150‑1 (MQ=255)
                                                                            tATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAg                                                                      <  1:3116646/150‑1 (MQ=255)
                                                                              ttCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgct                                                                    <  1:1842609/150‑1 (MQ=255)
                                                                               tCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctg                                                                   <  2:4511872/150‑1 (MQ=255)
                                                                                 tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATCATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                                 >  2:207776/1‑150 (MQ=255)
                                                                                 tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                                 >  2:2769476/1‑150 (MQ=255)
                                                                                  aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                                <  2:3006666/150‑1 (MQ=255)
                                                                                  aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                                >  1:3605946/1‑150 (MQ=255)
                                                                                      gttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCg                                                            >  1:4111034/1‑150 (MQ=255)
                                                                                       ttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                                           <  1:50467/150‑1 (MQ=255)
                                                                                        tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                                          <  1:4309466/150‑1 (MQ=255)
                                                                                         cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCTGCAACACTGCGTGTAATTTCTTCGTGTGGAATGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGata                                                         >  2:4095356/1‑150 (MQ=255)
                                                                                         cTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGata                                                         >  1:1706107/1‑150 (MQ=255)
                                                                                             aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                                     >  1:4024818/1‑150 (MQ=255)
                                                                                             aCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCa                                                     >  2:2392846/1‑150 (MQ=255)
                                                                                                 gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                 <  2:1393864/150‑1 (MQ=255)
                                                                                                  cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCCTGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                                <  1:2028707/150‑1 (MQ=255)
                                                                                                    gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                              <  1:837104/150‑1 (MQ=255)
                                                                                                     atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGATTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTTGCTGGCGATATGCATTCCCAGa                                             >  2:1202989/1‑150 (MQ=255)
                                                                                                      tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                            >  1:913993/1‑150 (MQ=255)
                                                                                                       atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                                           >  2:3857472/1‑150 (MQ=255)
                                                                                                       atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                                           >  1:2977850/1‑150 (MQ=255)
                                                                                                        ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTc                                          <  1:1340009/150‑1 (MQ=255)
                                                                                                            ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                                      >  1:4244948/1‑150 (MQ=255)
                                                                                                             tCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                                     <  1:3614104/150‑1 (MQ=255)
                                                                                                                aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                                  >  2:3098049/1‑150 (MQ=255)
                                                                                                                aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                                  >  2:668753/1‑150 (MQ=255)
                                                                                                                aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                                  <  2:840526/150‑1 (MQ=255)
                                                                                                                aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                                  <  1:1087490/150‑1 (MQ=255)
                                                                                                                 aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                                 <  1:4571259/150‑1 (MQ=255)
                                                                                                                  tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                                >  2:356449/1‑150 (MQ=255)
                                                                                                                   cgcgTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                               >  2:1153993/1‑150 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                            >  1:3050591/1‑150 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                            >  2:2165358/1‑150 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                            >  2:3655688/1‑150 (MQ=255)
                                                                                                                       ttCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt                           >  1:2272476/1‑150 (MQ=255)
                                                                                                                         cAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAg                         >  1:3951963/1‑150 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGTTg                      >  2:1128775/1‑150 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                      <  1:1005346/150‑1 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                      >  2:1612183/1‑150 (MQ=255)
                                                                                                                             tctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGa                     >  2:1725402/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGTCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    >  1:4209266/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    >  2:4510878/1‑150 (MQ=255)
                                                                                                                               tttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                   <  2:2304069/150‑1 (MQ=255)
                                                                                                                               tttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                   >  1:2519532/1‑150 (MQ=255)
                                                                                                                                ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc                  <  2:4473458/150‑1 (MQ=255)
                                                                                                                                 tcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg                 >  1:1242150/1‑150 (MQ=255)
                                                                                                                                    tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg              >  2:4236050/1‑150 (MQ=255)
                                                                                                                                    tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg              >  1:2123782/1‑150 (MQ=255)
                                                                                                                                    tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg              <  2:1257439/150‑1 (MQ=255)
                                                                                                                                     aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGc             >  2:1453541/1‑150 (MQ=255)
                                                                                                                                        cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg          >  2:749970/1‑150 (MQ=255)
                                                                                                                                         aaCCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGc         >  2:4097683/1‑150 (MQ=255)
                                                                                                                                           ccATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCtt       <  2:2881464/150‑1 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt   <  2:3896994/150‑1 (MQ=255)
                                                                                                                                                cTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAttt  >  2:4293085/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTT  >  CP009273/1048781‑1049074

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: