Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,050,071 G→T 8.4% R522S (CGC→AGC)  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,050,0710GT 8.4% 1138.7 / 51.5 308R522S (CGC→AGC) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base G (130/152);  new base T (14/12);  total (144/164)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 316 total aligned reads displayed.

GCTCAACAGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGG  >  CP009273/1049922‑1050218
                                                                                                                                                     |                                                                                                                                                   
gCTCAACAGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGgcg                                                                                                                                                     >  1:1939786/1‑150 (MQ=255)
    aaCAGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGt                                                                                                                                                 >  2:3887232/1‑150 (MQ=255)
     aCAGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGtt                                                                                                                                                >  1:4072685/1‑150 (MQ=255)
       aGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCa                                                                                                                                              <  2:1406383/150‑1 (MQ=255)
       aGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCa                                                                                                                                              <  1:3400567/150‑1 (MQ=255)
           tggttATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAgg                                                                                                                                          >  2:54506/1‑150 (MQ=255)
            ggttATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAggg                                                                                                                                         <  2:1702538/150‑1 (MQ=255)
              ttATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCt                                                                                                                                       >  2:2566625/1‑150 (MQ=255)
               tataACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTc                                                                                                                                      >  2:910967/1‑150 (MQ=255)
               tataACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTc                                                                                                                                      >  2:1645414/1‑150 (MQ=255)
               tataACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTc                                                                                                                                      >  1:1970850/1‑150 (MQ=255)
                ataACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCa                                                                                                                                     >  1:4349040/1‑150 (MQ=255)
                   aCCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCg                                                                                                                                  >  1:2285905/1‑150 (MQ=255)
                    ccTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGc                                                                                                                                 >  1:3018507/1‑150 (MQ=255)
                    ccTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCAACCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGc                                                                                                                                 >  2:2522476/1‑150 (MQ=255)
                      tGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTg                                                                                                                               >  1:2444185/1‑150 (MQ=255)
                      tGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTg                                                                                                                               <  2:3234837/150‑1 (MQ=255)
                         cGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGAcc                                                                                                                            <  2:1053008/150‑1 (MQ=255)
                          gAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCg                                                                                                                           >  2:2103159/1‑150 (MQ=255)
                          gAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCg                                                                                                                           >  2:63182/1‑150 (MQ=255)
                          gAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCg                                                                                                                           <  2:374856/150‑1 (MQ=255)
                          gAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCg                                                                                                                           <  2:2048383/150‑1 (MQ=255)
                          gAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGGTGCCAGGCGAATCTGGCGACCTATCACCCGTCCGCTCAGTAATTGCAGGGTACTTACCAGCAGCGGGCTCGGTTCAAAGTGCTCATCGCTGACCg                                                                                                                           >  2:4348979/1‑150 (MQ=255)
                            aTACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGaa                                                                                                                         >  1:2012437/1‑150 (MQ=255)
                            aTACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGaa                                                                                                                         >  1:161518/1‑150 (MQ=255)
                                   gtGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCt                                                                                                                  >  1:1899328/1‑150 (MQ=255)
                                   gtGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCt                                                                                                                  <  2:1981084/150‑1 (MQ=255)
                                   gtGGATCGCCCATTAACGCGCACGGCATATCATCGGAAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCt                                                                                                                  >  2:1648517/1‑150 (MQ=255)
                                       aTCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGcc                                                                                                              <  1:3183994/150‑1 (MQ=255)
                                       aTCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGcc                                                                                                              <  1:3229533/150‑1 (MQ=255)
                                       aTCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGcc                                                                                                              <  1:951032/150‑1 (MQ=255)
                                            ccATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTg                                                                                                         <  1:4014053/150‑1 (MQ=255)
                                             cATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGc                                                                                                        >  1:3070238/1‑150 (MQ=255)
                                              aTTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCt                                                                                                       >  1:1549662/1‑150 (MQ=255)
                                                tAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTc                                                                                                     <  2:4479196/150‑1 (MQ=255)
                                                tAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTc                                                                                                     >  1:1670993/1‑150 (MQ=255)
                                                   cgcgCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGAt                                                                                                  <  2:3494199/150‑1 (MQ=255)
                                                   cgcgCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGAt                                                                                                  <  2:3363403/150‑1 (MQ=255)
                                                    gcgcACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATa                                                                                                 <  1:4431636/150‑1 (MQ=255)
                                                    gcgcACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATa                                                                                                 <  2:3058840/150‑1 (MQ=255)
                                                     cgcACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAg                                                                                                >  2:3435408/1‑150 (MQ=255)
                                                       cACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCg                                                                                              >  2:2393219/1‑150 (MQ=255)
                                                            cATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGaata                                                                                         >  1:2239338/1‑150 (MQ=255)
                                                              tatCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGaataat                                                                                       >  2:2242811/1‑150 (MQ=255)
                                                              tatCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGaataat                                                                                       <  1:1474753/150‑1 (MQ=255)
                                                               atcatcGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATc                                                                                      <  1:2660651/150‑1 (MQ=255)
                                                                  tttggCATTTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCgag                                                                                   <  2:2942363/147‑1 (MQ=255)
                                                                  atcGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCgag                                                                                   <  1:3711482/150‑1 (MQ=255)
                                                                   tcGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCgaga                                                                                  <  1:3644643/150‑1 (MQ=255)
                                                                    cGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGaa                                                                                 <  2:1241866/150‑1 (MQ=255)
                                                                     ggCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAat                                                                                <  2:1779094/150‑1 (MQ=255)
                                                                       cAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAatat                                                                              >  2:2216777/1‑150 (MQ=255)
                                                                        aaTTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATc                                                                             <  1:1112703/150‑1 (MQ=255)
                                                                          ttGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGt                                                                           >  2:452239/1‑150 (MQ=255)
                                                                            gCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTg                                                                         <  1:1354560/150‑1 (MQ=255)
                                                                              cGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAg                                                                       >  2:2782687/1‑150 (MQ=255)
                                                                               gTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAgg                                                                      >  2:514376/1‑150 (MQ=255)
                                                                               gTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAgg                                                                      >  2:26014/1‑150 (MQ=255)
                                                                                 tGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGAt                                                                    <  2:3599874/150‑1 (MQ=255)
                                                                                  gCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATg                                                                   <  1:4272928/150‑1 (MQ=255)
                                                                                    cAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGt                                                                 >  2:4426628/1‑150 (MQ=255)
                                                                                    cAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGt                                                                 <  1:2527992/150‑1 (MQ=255)
                                                                                      ggCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCa                                                               <  2:2194882/150‑1 (MQ=255)
                                                                                       gCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAg                                                              <  2:592325/150‑1 (MQ=255)
                                                                                        cGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGc                                                             >  1:3278534/1‑150 (MQ=255)
                                                                                          aaTCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCaa                                                           <  1:4034033/150‑1 (MQ=255)
                                                                                          aaTCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCaa                                                           <  2:4300364/150‑1 (MQ=255)
                                                                                          aaTCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCaa                                                           <  2:4300351/150‑1 (MQ=255)
                                                                                                   aCCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCAcc                                                  <  2:2876804/150‑1 (MQ=255)
                                                                                                    ccTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCa                                                 >  1:96984/1‑150 (MQ=255)
                                                                                                     cTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCag                                                >  2:3429572/1‑150 (MQ=255)
                                                                                                     cTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCag                                                >  2:1156563/1‑150 (MQ=255)
                                                                                                      tttCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCaga                                               >  2:3241363/1‑150 (MQ=255)
                                                                                                         cACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTc                                            >  2:4159118/1‑150 (MQ=255)
                                                                                                         cACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTc                                            >  2:1932102/1‑150 (MQ=255)
                                                                                                          aCCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCa                                           >  2:2713652/1‑150 (MQ=255)
                                                                                                          aCCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCa                                           <  1:141937/150‑1 (MQ=255)
                                                                                                          aCCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCa                                           <  1:4138538/150‑1 (MQ=255)
                                                                                                                ccGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAAtt                                     >  1:3429455/1‑150 (MQ=255)
                                                                                                                  gCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGc                                   >  2:2466839/1‑150 (MQ=255)
                                                                                                                    tCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGccc                                 <  1:2855861/150‑1 (MQ=255)
                                                                                                                    tCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGccc                                 <  1:2695777/150‑1 (MQ=255)
                                                                                                                      aTTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGc                               <  1:565607/150‑1 (MQ=255)
                                                                                                                      aTTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGc                               <  1:2972140/150‑1 (MQ=255)
                                                                                                                        tAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCaa                             >  2:1764607/1‑150 (MQ=255)
                                                                                                                        tAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCaa                             >  2:2642086/1‑150 (MQ=255)
                                                                                                                        tAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCaa                             >  2:4112891/1‑150 (MQ=255)
                                                                                                                           ttGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAatc                          <  2:32160/150‑1 (MQ=255)
                                                                                                                            tGCAGGGTACTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAatca                         <  1:4354077/150‑1 (MQ=255)
                                                                                                                                gggTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCAc                     <  2:2772492/150‑1 (MQ=255)
                                                                                                                                    aCTTTCCAGCAGCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTg                 <  1:2626144/150‑1 (MQ=255)
                                                                                                                                    aCTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTg                 >  2:2663351/1‑150 (MQ=255)
                                                                                                                                      tttCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTggg               <  1:1754743/150‑1 (MQ=255)
                                                                                                                                       ttCCATCAGCGGGCGCGGTTCAAAGGGATCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGc              <  1:1014454/150‑1 (MQ=255)
                                                                                                                                       ttCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGc              <  1:895278/150‑1 (MQ=255)
                                                                                                                                         ccAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGt            <  1:4250361/150‑1 (MQ=255)
                                                                                                                                         ccAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGt            <  1:4250668/150‑1 (MQ=255)
                                                                                                                                          cagcagCGGGCTCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGtt           <  1:2534012/150‑1 (MQ=255)
                                                                                                                                          cagcagCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGtt           >  2:3092575/1‑150 (MQ=255)
                                                                                                                                                   gcgcGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCggg  <  1:2582850/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTCAACAGGTTGGTTATAACCTGACGAATACGTCGTGGATCGCCCATTAACGCGCACGGCATATCATCGGCAATTGCCGTTGCCAGGCGAATCGGGCGACCTTTCACCCGTCCGCTCATTAATTGCAGGGTACTTTCCAGCAGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGG  >  CP009273/1049922‑1050218

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: