Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,928 C→A 89.8% 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,9280CA89.8% 856.4 / 87.0 302E903* (GAA→TAA) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base C (18/13);  new base A (136/135);  total (154/148)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTT  >  CP009273/1048781‑1049074
                                                                                                                                                   |                                                                                                                                                  
ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACa                                                                                                                                                  <  1:3109228/150‑1 (MQ=255)
 gccgccGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACag                                                                                                                                                 <  1:348757/150‑1 (MQ=255)
  ccgccgCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACaga                                                                                                                                                <  1:3723085/150‑1 (MQ=255)
  ccgccgCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACaga                                                                                                                                                >  2:4695270/1‑150 (MQ=255)
   cgccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACagag                                                                                                                                               >  2:447594/1‑150 (MQ=255)
     ccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAGCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCg                                                                                                                                             >  1:2349622/1‑150 (MQ=255)
     ccgcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCg                                                                                                                                             <  2:1982062/150‑1 (MQ=255)
       gcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGc                                                                                                                                           <  1:1833697/150‑1 (MQ=255)
        cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCa                                                                                                                                          <  1:2046100/150‑1 (MQ=255)
        cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGACTTTCTTATGCAACCATGCTTACAGAGCGGCa                                                                                                                                          <  2:2476245/150‑1 (MQ=255)
          gCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAac                                                                                                                                        <  1:4271230/150‑1 (MQ=255)
           ccAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAaca                                                                                                                                       >  1:3423905/1‑150 (MQ=255)
           ccAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAaca                                                                                                                                       >  2:1608565/1‑150 (MQ=255)
               cGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGc                                                                                                                                   >  1:4018155/1‑150 (MQ=255)
                  gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtg                                                                                                                                >  1:950203/1‑150 (MQ=255)
                      ttttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAt                                                                                                                            >  2:4140610/1‑150 (MQ=255)
                        ttATTCTCTGCCATGATTTTGCCAGTCCGTTCAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAttt                                                                                                                          <  2:2887160/150‑1 (MQ=255)
                        ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAttt                                                                                                                          <  2:3323398/150‑1 (MQ=255)
                        ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAAttt                                                                                                                          >  2:3002911/1‑150 (MQ=255)
                           ttCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATttctt                                                                                                                       <  1:3637200/150‑1 (MQ=255)
                                 gCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGtgtg                                                                                                                 >  1:477547/1‑150 (MQ=255)
                                        ttttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggg                                                                                                          >  2:4913405/1‑150 (MQ=255)
                                          ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggggg                                                                                                        >  1:518044/1‑150 (MQ=255)
                                          ttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAgggggg                                                                                                        >  1:4907741/1‑150 (MQ=255)
                                           tGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                       <  2:276839/150‑1 (MQ=255)
                                           tGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                       >  2:5130069/1‑150 (MQ=255)
                                            gCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCa                                                                                                      <  1:141893/150‑1 (MQ=255)
                                             ccAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAg                                                                                                     <  1:3192347/150‑1 (MQ=255)
                                              cAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGa                                                                                                    <  2:92672/150‑1 (MQ=255)
                                               aGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGAt                                                                                                   >  1:2946651/1‑150 (MQ=255)
                                                  ccGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAAt                                                                                                <  2:3209774/150‑1 (MQ=255)
                                                  ccGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAATGGGGGCAGATAAt                                                                                                <  1:2104416/150‑1 (MQ=255)
                                                  ccGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAAt                                                                                                <  2:2303371/150‑1 (MQ=255)
                                                   cGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                               >  1:2567847/1‑150 (MQ=255)
                                                   cGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                               <  2:4683484/150‑1 (MQ=255)
                                                    gTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATgg                                                                                              >  2:3423987/1‑150 (MQ=255)
                                                         gTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgct                                                                                         <  1:124579/150‑1 (MQ=255)
                                                            gTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctgct                                                                                      <  2:362214/150‑1 (MQ=255)
                                                               taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAACTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                   <  2:208318/150‑1 (MQ=255)
                                                               taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                   >  1:4334958/1‑150 (MQ=255)
                                                               taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                   >  1:634566/1‑150 (MQ=255)
                                                               taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACAATGCTTACAGAGCTGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                   >  1:3214046/1‑150 (MQ=255)
                                                                   aaaTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTATGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                               >  2:4025091/1‑150 (MQ=255)
                                                                    aaTGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                              <  2:1368688/150‑1 (MQ=255)
                                                                      tGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                            >  2:4239766/1‑150 (MQ=255)
                                                                        ccACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                                          >  2:3456426/1‑150 (MQ=255)
                                                                         cACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                                         >  2:2416835/1‑150 (MQ=255)
                                                                         cACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                                         >  1:82163/1‑150 (MQ=255)
                                                                            tATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAg                                                                      >  2:124152/1‑150 (MQ=255)
                                                                             aTTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgc                                                                     >  1:2314205/1‑150 (MQ=255)
                                                                             aTTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgc                                                                     <  2:2458090/150‑1 (MQ=255)
                                                                              ttCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgct                                                                    <  2:4232574/150‑1 (MQ=255)
                                                                              ttCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgct                                                                    <  2:5020364/150‑1 (MQ=255)
                                                                                 tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                                 >  1:489127/1‑150 (MQ=255)
                                                                                  aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                                <  1:1940900/150‑1 (MQ=255)
                                                                                  aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                                <  1:5077975/150‑1 (MQ=255)
                                                                                           tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                                       >  1:5022249/1‑150 (MQ=255)
                                                                                            aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGc                                                      >  2:3842708/1‑150 (MQ=255)
                                                                                              cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                                    >  1:1005845/1‑150 (MQ=255)
                                                                                              cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                                    <  2:603046/150‑1 (MQ=255)
                                                                                               tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                                   <  1:4598530/150‑1 (MQ=255)
                                                                                                 gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                 >  1:2630392/1‑150 (MQ=255)
                                                                                                  cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                                >  2:462811/1‑150 (MQ=255)
                                                                                                   tGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                               >  1:723373/1‑150 (MQ=255)
                                                                                                   tGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                               <  1:269390/150‑1 (MQ=255)
                                                                                                     atatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGa                                             <  2:1821645/150‑1 (MQ=255)
                                                                                                      tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                            >  1:1631601/1‑150 (MQ=255)
                                                                                                       atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                                           >  1:4657128/1‑150 (MQ=255)
                                                                                                        ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTc                                          <  2:1041407/150‑1 (MQ=255)
                                                                                                          gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCAATCCTCGa                                        >  2:896250/1‑150 (MQ=255)
                                                                                                           aTTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAg                                       <  2:2834876/150‑1 (MQ=255)
                                                                                                                 aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                                 >  1:862976/1‑150 (MQ=255)
                                                                                                                   cgcgTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                               >  2:4263381/1‑150 (MQ=255)
                                                                                                                   cgcgTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                               >  1:4943177/1‑150 (MQ=255)
                                                                                                                   cgcgTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCt                               <  2:1590430/150‑1 (MQ=255)
                                                                                                                    gcgTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                              <  2:1393404/150‑1 (MQ=255)
                                                                                                                    gcgTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtt                              >  1:2609286/1‑150 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCTACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                            <  1:2059249/150‑1 (MQ=255)
                                                                                                                      gTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                            >  1:3881739/1‑150 (MQ=255)
                                                                                                                          aGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                        >  1:837118/1‑150 (MQ=255)
                                                                                                                           ggTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGcn                       <  2:3809135/150‑2 (MQ=255)
                                                                                                                            gTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg                      >  2:3984574/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    >  1:3011774/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    <  1:4782738/150‑1 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    >  2:1884759/1‑150 (MQ=255)
                                                                                                                              ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                    >  1:1770713/1‑150 (MQ=255)
                                                                                                                               tttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                   <  1:4392097/150‑1 (MQ=255)
                                                                                                                                ttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgc                  <  2:3524861/150‑1 (MQ=255)
                                                                                                                                 tcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg                 <  1:35700/150‑1 (MQ=255)
                                                                                                                                 tcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg                 >  2:2436888/1‑150 (MQ=255)
                                                                                                                                    tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg              >  1:2481108/1‑150 (MQ=255)
                                                                                                                                     aTGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCAATATGCATTCCCAGACTCGAGCACCTGCTTTTTAGCTGATGCGCGGc             >  2:1779716/1‑150 (MQ=255)
                                                                                                                                       gCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCAc           <  2:2587017/150‑1 (MQ=255)
                                                                                                                                        cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg          <  2:3649439/150‑1 (MQ=255)
                                                                                                                                         aaCCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGc         >  2:2539809/1‑150 (MQ=255)
                                                                                                                                           ccATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCCGGTGTGCGCCTAGCTGGCTGGCGATATGCATTCCCAGACGCGAGCACCTGCTGTTTCGCTGATGCGCGGCACGCtg       >  2:3130951/1‑149 (MQ=255)
                                                                                                                                             aTGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTa     >  2:3044884/1‑150 (MQ=255)
                                                                                                                                               gCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtt   <  1:1506710/150‑1 (MQ=255)
                                                                                                                                                cTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAttt  <  2:3560291/150‑1 (MQ=255)
                                                                                                                                                cTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAttt  >  2:2217699/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: