Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 590,565 A→T 5.7% V112E (GTG→GAG)  cusR ← response regulator in two‑component regulatory system with CusS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273590,5650AT 5.7% 1064.9 / 36.4 341V112E (GTG→GAG) cusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base A (165/156);  new base T (5/15);  total (170/171)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

AGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCAC  >  CP009273/590416‑590713
                                                                                                                                                     |                                                                                                                                                    
 gaaCTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCAc                                                                                                                                                     >  1:2709419/1‑150 (MQ=255)
  aaCTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCAcc                                                                                                                                                    <  2:2034190/150‑1 (MQ=255)
    cTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCcg                                                                                                                                                  >  1:4758401/1‑150 (MQ=255)
        agcagAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGcc                                                                                                                                              >  2:2524500/1‑150 (MQ=255)
           agagTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCcagc                                                                                                                                           <  2:2020275/150‑1 (MQ=255)
           agagTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCcagc                                                                                                                                           >  1:2480596/1‑150 (MQ=255)
            gagTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCTCCCGCGCcagca                                                                                                                                          >  2:2135723/1‑150 (MQ=255)
             agTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCcagcag                                                                                                                                         >  1:1278366/1‑150 (MQ=255)
              gTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGt                                                                                                                                        >  1:2842243/1‑150 (MQ=255)
               tAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGtt                                                                                                                                       <  1:101683/150‑1 (MQ=255)
               tAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGtt                                                                                                                                       >  1:3808693/1‑150 (MQ=255)
               tAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGtt                                                                                                                                       >  1:3861445/1‑150 (MQ=255)
               tAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACAAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGtt                                                                                                                                       >  1:908723/1‑150 (MQ=255)
                  aCTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAg                                                                                                                                    <  1:2094526/150‑1 (MQ=255)
                     ctTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCaa                                                                                                                                 >  2:3286937/1‑150 (MQ=255)
                      tttACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCaaa                                                                                                                                <  2:699009/150‑1 (MQ=255)
                      tttACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCaaa                                                                                                                                >  2:1640149/1‑150 (MQ=255)
                           tGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGa                                                                                                                           <  2:4352777/150‑1 (MQ=255)
                            ggTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGaa                                                                                                                          >  1:685771/1‑150 (MQ=255)
                             gTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAAt                                                                                                                         <  1:4326810/150‑1 (MQ=255)
                               cAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATgg                                                                                                                       >  1:1090125/1‑150 (MQ=255)
                                  aGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAGCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCTCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCtt                                                                                                                    >  1:821228/1‑150 (MQ=255)
                                      aTGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGACAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCAcc                                                                                                                <  1:1488571/150‑1 (MQ=255)
                                       tGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCa                                                                                                               <  1:954579/150‑1 (MQ=255)
                                        gcgcGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAg                                                                                                              >  2:403208/1‑150 (MQ=255)
                                        gcgcGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAg                                                                                                              <  1:4338674/150‑1 (MQ=255)
                                        gcgcGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAg                                                                                                              >  2:5018067/1‑150 (MQ=255)
                                         cgcgTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAgg                                                                                                             <  2:5237716/150‑1 (MQ=255)
                                           cgTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTa                                                                                                           <  2:2864708/150‑1 (MQ=255)
                                           cgTGCCGCTGCGGGTGACTTTGCGGCAGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGAGAATCACCGCCGCCCCGCGCCGCACTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTACACCAGGTa                                                                                                           >  2:2321297/1‑150 (MQ=255)
                                            gTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAg                                                                                                          >  1:310095/1‑150 (MQ=255)
                                            gTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAg                                                                                                          >  2:3978006/1‑150 (MQ=255)
                                                cGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGACTGGCTTCACCAGGTAGTCAt                                                                                                      >  1:4632139/1‑150 (MQ=255)
                                                 gCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATc                                                                                                     >  2:3538617/1‑150 (MQ=255)
                                                     cGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGcc                                                                                                 >  2:4528738/1‑150 (MQ=255)
                                                       ggTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGcccc                                                                                               <  1:3413516/150‑1 (MQ=255)
                                                       ggTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGcccc                                                                                               <  2:5092080/150‑1 (MQ=255)
                                                         tGACTTTGCGGCTGACGAGATCGACCATCAAAGCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCa                                                                                             <  1:1932711/150‑1 (MQ=255)
                                                            cTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACt                                                                                          >  2:3165303/1‑150 (MQ=255)
                                                            cTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACt                                                                                          <  2:3432130/150‑1 (MQ=255)
                                                              ttGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCGGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTcc                                                                                        >  2:1523829/1‑150 (MQ=255)
                                                              ttGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTcc                                                                                        <  1:4440031/150‑1 (MQ=255)
                                                                   gCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGccc                                                                                   >  2:1676588/1‑150 (MQ=255)
                                                                    cTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGcccc                                                                                  <  2:1318635/150‑1 (MQ=255)
                                                                     tGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCt                                                                                 <  2:2571230/150‑1 (MQ=255)
                                                                      gACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCtt                                                                                <  1:5176165/150‑1 (MQ=255)
                                                                      gACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCtt                                                                                <  2:4318710/150‑1 (MQ=255)
                                                                         gagaTCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCTCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAc                                                                             <  2:795593/150‑1 (MQ=255)
                                                                         gagaTCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAc                                                                             >  1:2183443/1‑150 (MQ=255)
                                                                         gagaTCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAc                                                                             <  2:3980629/150‑1 (MQ=255)
                                                                         gagaTCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAc                                                                             <  2:1583435/150‑1 (MQ=255)
                                                                          agaTCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAcg                                                                            >  1:3633559/1‑150 (MQ=255)
                                                                                 ccATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGtt                                                                     <  2:4490850/150‑1 (MQ=255)
                                                                                   aTCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCa                                                                   <  1:4478736/150‑1 (MQ=255)
                                                                                   aTCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCa                                                                   >  2:190572/1‑150 (MQ=255)
                                                                                    tCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCaa                                                                  >  2:2100884/1‑150 (MQ=255)
                                                                                    tCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCaa                                                                  <  1:1060585/150‑1 (MQ=255)
                                                                                      aaaTCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATg                                                                >  2:1072248/1‑150 (MQ=255)
                                                                                        aTCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGt                                                              >  2:4057468/1‑150 (MQ=255)
                                                                                         tCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTg                                                             >  1:604531/1‑150 (MQ=255)
                                                                                         tCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTg                                                             <  1:4497042/150‑1 (MQ=255)
                                                                                            gCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGGTCAATGGTGCCa                                                          >  2:1810834/1‑150 (MQ=255)
                                                                                                    aaaCTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGt                                                  >  2:1482962/1‑150 (MQ=255)
                                                                                                      aCTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTaa                                                <  1:1543605/150‑1 (MQ=255)
                                                                                                         gACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCa                                             <  1:368212/150‑1 (MQ=255)
                                                                                                         gACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCa                                             <  1:5063557/150‑1 (MQ=255)
                                                                                                             ttCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAg                                         <  2:1986277/150‑1 (MQ=255)
                                                                                                                gATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAt                                      <  2:3342225/150‑1 (MQ=255)
                                                                                                                gATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAt                                      <  2:4904527/150‑1 (MQ=255)
                                                                                                                gATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAt                                      <  1:4087243/150‑1 (MQ=255)
                                                                                                                gATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAt                                      >  1:1987034/1‑150 (MQ=255)
                                                                                                                 aTAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCTCCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATc                                     >  2:3895703/1‑150 (MQ=255)
                                                                                                                  tAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCg                                    <  1:2797815/150‑1 (MQ=255)
                                                                                                                  tAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCg                                    <  1:3206365/150‑1 (MQ=255)
                                                                                                                    aTCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGc                                  >  2:1287334/1‑150 (MQ=255)
                                                                                                                     tCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCa                                 >  2:23694/1‑150 (MQ=255)
                                                                                                                      cACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGTAt                                <  2:3989217/150‑1 (MQ=255)
                                                                                                                         cgccgcCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATccc                             >  2:4154375/1‑150 (MQ=255)
                                                                                                                          gccgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCt                            <  2:3263441/150‑1 (MQ=255)
                                                                                                                          gccgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCt                            <  1:2529506/150‑1 (MQ=255)
                                                                                                                           ccgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCtt                           <  1:1805812/150‑1 (MQ=255)
                                                                                                                            cgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCttt                          <  2:4729407/150‑1 (MQ=255)
                                                                                                                            cgccCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCttt                          <  2:3618426/150‑1 (MQ=255)
                                                                                                                              ccccGCGCCGCAGTAATGTGCGCTCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTAt                        <  2:4750824/150‑1 (MQ=255)
                                                                                                                              ccccGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTAt                        <  2:1090125/150‑1 (MQ=255)
                                                                                                                              ccccGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTAt                        >  1:540803/1‑150 (MQ=255)
                                                                                                                                 cgcgCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTgg                     <  1:4537068/150‑1 (MQ=255)
                                                                                                                                   cgccgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCg                   <  2:2481666/150‑1 (MQ=255)
                                                                                                                                   cgccgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCg                   <  2:2709419/150‑1 (MQ=255)
                                                                                                                                    gccgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCgg                  >  1:746124/1‑150 (MQ=255)
                                                                                                                                    gccgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCgg                  <  2:2842243/150‑1 (MQ=255)
                                                                                                                                      cgcAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAg                >  2:3764636/1‑150 (MQ=255)
                                                                                                                                        cAGTAATGTGCGCTCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCg              >  1:4183220/1‑150 (MQ=255)
                                                                                                                                         aGTAATGTGCGCTCCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGt             <  1:461932/150‑1 (MQ=255)
                                                                                                                                         aGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGt             >  1:4052807/1‑150 (MQ=255)
                                                                                                                                          gTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTa            >  1:1145576/1‑150 (MQ=255)
                                                                                                                                           tAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTGTATTGGCGGAGCGTaa           >  2:5132915/1‑150 (MQ=255)
                                                                                                                                            aaTGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGAGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGAGATGTTCAATGGTGCCAAGCGAGGTAAGAAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAAc          <  1:2894147/150‑1 (MQ=255)
                                                                                                                                               gtgCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATg       <  2:2909477/150‑1 (MQ=255)
                                                                                                                                                   gcACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCa   >  2:4270935/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCATCCCTTTATTGGCGGAGCGTAACATGCGCAC  >  CP009273/590416‑590713

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: