Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,313,924 C→A 30.8% R135S (CGC→AGC)  atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,313,9240CA30.8% 661.8 / 281.5 341R135S (CGC→AGC) atoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base C (104/131);  new base A (53/53);  total (157/184)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

TTGCCCCCATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGATATTCGCCGCCA  >  CP009273/2313775‑2314072
                                                                                                                                                     |                                                                                                                                                    
ttGCCCCCATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTc                                                                                                                                                      <  2:2240897/150‑1 (MQ=255)
ttGCCCCCATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTc                                                                                                                                                      <  2:3013516/150‑1 (MQ=255)
   cccccATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTcgcg                                                                                                                                                   >  2:334923/1‑150 (MQ=255)
    ccccATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGa                                                                                                                                                  >  1:521362/1‑150 (MQ=255)
     cccATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGaa                                                                                                                                                 >  1:2707687/1‑150 (MQ=255)
       cATTACCGAAAAGATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGt                                                                                                                                               <  1:4089257/150‑1 (MQ=255)
           aCCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATg                                                                                                                                           <  1:3893396/150‑1 (MQ=255)
             cGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCg                                                                                                                                         <  2:4202306/150‑1 (MQ=255)
              gAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGt                                                                                                                                        <  1:1412149/150‑1 (MQ=255)
               aaaaTATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTa                                                                                                                                       <  2:506159/150‑1 (MQ=255)
                aaaTATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTAc                                                                                                                                      <  1:3703007/150‑1 (MQ=255)
                   tatCACTCACGCCTTCCCTGGAATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCATGGTCGCGAAGTGATGCGTACaaa                                                                                                                                   <  1:4755561/150‑1 (MQ=255)
                    atCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAAt                                                                                                                                  <  1:3853155/150‑1 (MQ=255)
                     tcactcacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATa                                                                                                                                 <  2:1288668/150‑1 (MQ=255)
                      cactcacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATAc                                                                                                                                >  1:1822617/1‑150 (MQ=255)
                        ctcacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATAccc                                                                                                                              <  1:163904/150‑1 (MQ=255)
                         tcacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATAcccc                                                                                                                             <  1:1246990/150‑1 (MQ=255)
                         tcacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATAcccc                                                                                                                             >  1:3459342/1‑150 (MQ=255)
                          cacGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCt                                                                                                                            >  2:559615/1‑150 (MQ=255)
                           acGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCtt                                                                                                                           <  2:4491727/150‑1 (MQ=255)
                           acGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCtt                                                                                                                           <  2:4973446/150‑1 (MQ=255)
                            cGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCttt                                                                                                                          <  2:2409919/150‑1 (MQ=255)
                             gCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTa                                                                                                                         <  1:1922885/150‑1 (MQ=255)
                              ccTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAg                                                                                                                        <  1:1266540/150‑1 (MQ=255)
                              ccTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAg                                                                                                                        <  2:2751920/150‑1 (MQ=255)
                               cTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGt                                                                                                                       <  1:4533046/150‑1 (MQ=255)
                                  cccTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTa                                                                                                                    <  1:5055787/150‑1 (MQ=255)
                                          tCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggca                                                                                                            <  2:4591964/150‑1 (MQ=255)
                                          tCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggca                                                                                                            >  2:1002849/1‑150 (MQ=255)
                                           cGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcag                                                                                                           <  1:3052316/150‑1 (MQ=255)
                                           cGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcag                                                                                                           >  1:3122695/1‑150 (MQ=255)
                                            ggTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcagg                                                                                                          >  2:3943191/1‑150 (MQ=255)
                                              tGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCTCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcaggca                                                                                                        <  2:185185/150‑1 (MQ=255)
                                              tGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcaggca                                                                                                        <  2:4933875/150‑1 (MQ=255)
                                                cTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTcaggcaggcagg                                                                                                      >  1:4646963/1‑150 (MQ=255)
                                                 tGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGt                                                                                                     <  2:1962190/150‑1 (MQ=255)
                                                 tGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGt                                                                                                     <  1:1554674/150‑1 (MQ=255)
                                                  ggTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTg                                                                                                    >  2:775359/1‑150 (MQ=255)
                                                   gTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTgc                                                                                                   >  2:4010551/1‑150 (MQ=255)
                                                    ttattaCAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTgcg                                                                                                  >  2:3099237/1‑150 (MQ=255)
                                                     tattaCAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTgcgc                                                                                                 >  2:4695201/1‑150 (MQ=255)
                                                      attaCAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGcgcg                                                                                                >  1:1731900/1‑150 (MQ=255)
                                                         acaacaAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCgcggcg                                                                                             <  1:4323362/150‑1 (MQ=255)
                                                         acaacaAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCgcggcg                                                                                             >  2:593975/1‑150 (MQ=255)
                                                            acaaAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGata                                                                                          >  2:1094328/1‑150 (MQ=255)
                                                                aaTGCTGGATGCGATAATCACCTACGCTCCTTCTGCTCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTTTTCAGGCAGGCAGGTGCGCGGCGATAtttt                                                                                      <  1:1335995/150‑1 (MQ=255)
                                                                aaTGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATAtttt                                                                                      <  1:1873831/150‑1 (MQ=255)
                                                                      ggATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTc                                                                                >  2:2719656/1‑150 (MQ=255)
                                                                      ggATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTc                                                                                <  1:1990415/150‑1 (MQ=255)
                                                                        aTGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCaa                                                                              <  1:3592788/150‑1 (MQ=255)
                                                                        aTGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCaa                                                                              <  1:2123831/150‑1 (MQ=255)
                                                                         tGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAAt                                                                             <  1:1000941/150‑1 (MQ=255)
                                                                              tAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTcc                                                                        <  1:1523040/150‑1 (MQ=255)
                                                                                aTCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTcccc                                                                      <  1:3881132/150‑1 (MQ=255)
                                                                                  cACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCAt                                                                    <  1:583976/150‑1 (MQ=255)
                                                                                   aCCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCAtt                                                                   >  2:4125390/1‑150 (MQ=255)
                                                                                   aCCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCAtt                                                                   <  2:650778/150‑1 (MQ=255)
                                                                                     cTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGa                                                                 >  1:1138394/1‑150 (MQ=255)
                                                                                       aCGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGc                                                               >  1:3968160/1‑150 (MQ=255)
                                                                                             cTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATg                                                         <  2:951453/150‑1 (MQ=255)
                                                                                              ttCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATgg                                                        >  2:2332920/1‑150 (MQ=255)
                                                                                              ttCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATgg                                                        <  2:2129101/150‑1 (MQ=255)
                                                                                                cAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTg                                                      <  1:3354835/150‑1 (MQ=255)
                                                                                                cAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTg                                                      >  2:2916853/1‑150 (MQ=255)
                                                                                                    gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAAt                                                  <  2:1921197/150‑1 (MQ=255)
                                                                                                    gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAAt                                                  <  2:2473939/150‑1 (MQ=255)
                                                                                                    gcTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAAt                                                  <  2:3402994/150‑1 (MQ=255)
                                                                                                       atatCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCt                                               <  2:16870/150‑1 (MQ=255)
                                                                                                       atatCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCt                                               >  2:5045072/1‑150 (MQ=255)
                                                                                                        tatCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTc                                              <  1:2569163/150‑1 (MQ=255)
                                                                                                         atCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCg                                             >  1:2441768/1‑150 (MQ=255)
                                                                                                             gAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCta                                         <  1:322168/150‑1 (MQ=255)
                                                                                                             gAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCta                                         >  2:3587783/1‑150 (MQ=255)
                                                                                                              aataatGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCtat                                        <  2:3201900/150‑1 (MQ=255)
                                                                                                               ataatGTCGGCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCtata                                       <  2:1199964/150‑1 (MQ=255)
                                                                                                                  atGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCt                                    >  2:597023/1‑150 (MQ=255)
                                                                                                                  atGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCt                                    >  2:3240273/1‑150 (MQ=255)
                                                                                                                   tGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTg                                   <  1:2658361/150‑1 (MQ=255)
                                                                                                                       ggCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGcc                               <  1:875770/150‑1 (MQ=255)
                                                                                                                       ggCGTTACCATTGCCGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGcc                               >  2:2163704/1‑150 (MQ=255)
                                                                                                                           ttACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATg                           <  1:1041875/150‑1 (MQ=255)
                                                                                                                              ccATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAAt                        <  2:3421286/150‑1 (MQ=255)
                                                                                                                              ccATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAAt                        >  1:2888399/1‑150 (MQ=255)
                                                                                                                              ccATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAAt                        <  1:4446953/150‑1 (MQ=255)
                                                                                                                                aTTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTa                      <  2:3678168/150‑1 (MQ=255)
                                                                                                                                 ttGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTaa                     <  2:1455738/150‑1 (MQ=255)
                                                                                                                                  tGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAAc                    >  1:2472092/1‑150 (MQ=255)
                                                                                                                                    ccGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCACTGAATTAACCg                  >  2:3757815/1‑150 (MQ=255)
                                                                                                                                    ccGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCg                  <  1:4972027/150‑1 (MQ=255)
                                                                                                                                    ccGCAGATCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCg                  <  2:2534510/150‑1 (MQ=255)
                                                                                                                                      gCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGaa                >  2:4151509/1‑150 (MQ=255)
                                                                                                                                        aGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGa              >  1:946001/1‑150 (MQ=255)
                                                                                                                                        aGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGa              >  2:5086338/1‑150 (MQ=255)
                                                                                                                                         gATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGat             <  1:3760284/150‑1 (MQ=255)
                                                                                                                                          aTCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGata            >  1:2965437/1‑150 (MQ=255)
                                                                                                                                           tCACCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGatat           <  2:3831287/150‑1 (MQ=255)
                                                                                                                                             aCCCTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGATATTc         >  1:1305475/1‑150 (MQ=255)
                                                                                                                                              accTGGTCGCGAAGTGATGCGTACAAATACACCTTTAGTTTATTCAGGCAGGCAGGTTCGCGGCGAGATTTTGAATTCAATGCTCCACATTGAGCGTAATGGTGAAATCCTCGTCTATATCTGGGCCAAAGACTAAACCGAAGATATTcg        >  2:4438271/2‑150 (MQ=255)
                                                                                                                                                cTGGTAGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGATATTcgcc      >  2:1261915/1‑150 (MQ=255)
                                                                                                                                                   gTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGATATTCgccgcc   <  2:4216996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGCCCCCATTACCGAAAATATCACTCACGCCTTCCCTGGCATCGGTGCTGGTTATTACAACAAAATGCTGGATGCGATAATCACCTACGCGCCTTCAGCGCTATATCAGAATAATGTCGGCGTTACCATTGCCGCAGATCACCCTGGTCGCGAAGTGATGCGTACAAATACCCCTTTAGTTTATTCAGGCAGGCAGGTGCGCGGCGATATTTTGAATTCAATGCTCCCCATTGAGCGTAATGGTGAAATCCTCGGCTATATCTGGGCCAATGAATTAACCGAAGATATTCGCCGCCA  >  CP009273/2313775‑2314072

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: