Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,963,075 G→T 77.1% R183R (CGG→AGG)  cheR ← chemotaxis regulator, protein‑glutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,963,0750GT77.1% 510.5 / 226.0 319R183R (CGG→AGG) cheRchemotaxis regulator, protein‑glutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/30);  new base T (124/122);  total (167/152)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

Only 100 of 321 total aligned reads displayed.

CGTTACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCG  >  CP009273/1962928‑1963221
                                                                                                                                                   |                                                                                                                                                  
cGTTACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGtt                                                                                                                                                  <  1:1212881/150‑1 (MQ=255)
   tACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTGTGCa                                                                                                                                               <  2:2692214/150‑1 (MQ=255)
    aCGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAg                                                                                                                                              <  1:1194258/150‑1 (MQ=255)
     cGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGt                                                                                                                                             >  2:2029825/1‑150 (MQ=255)
      gACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGtt                                                                                                                                            >  1:4903454/1‑150 (MQ=255)
      gACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGtt                                                                                                                                            >  2:555518/1‑150 (MQ=255)
        cAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTtgc                                                                                                                                          >  2:2848443/1‑150 (MQ=255)
         agaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTtgct                                                                                                                                         <  2:1784951/150‑1 (MQ=255)
         agaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTtgct                                                                                                                                         <  2:2080084/150‑1 (MQ=255)
         agaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTtgct                                                                                                                                         <  1:1040593/150‑1 (MQ=255)
         agaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTtgct                                                                                                                                         >  1:3196243/1‑150 (MQ=255)
          gaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTtgctg                                                                                                                                        <  1:2226165/150‑1 (MQ=255)
            agaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTgcg                                                                                                                                      <  1:1910694/150‑1 (MQ=255)
                 gCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTc                                                                                                                                 >  1:3902339/1‑150 (MQ=255)
                  cATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCa                                                                                                                                <  1:4018969/150‑1 (MQ=255)
                   aTCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAg                                                                                                                               >  1:4556643/1‑150 (MQ=255)
                    tCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAgg                                                                                                                              >  1:2976471/1‑150 (MQ=255)
                    tCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAgg                                                                                                                              >  2:4000905/1‑150 (MQ=255)
                    tCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAgg                                                                                                                              <  1:3685396/150‑1 (MQ=255)
                     cAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGt                                                                                                                             <  1:4079502/150‑1 (MQ=255)
                     cAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGt                                                                                                                             >  1:3866712/1‑150 (MQ=255)
                      aaaCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGtt                                                                                                                            >  1:4354772/1‑150 (MQ=255)
                       aaCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGttt                                                                                                                           >  1:1097894/1‑150 (MQ=255)
                          ggCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTc                                                                                                                        >  2:3340044/1‑150 (MQ=255)
                           gCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCa                                                                                                                       <  1:672123/150‑1 (MQ=255)
                               cGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAActc                                                                                                                   <  1:998509/150‑1 (MQ=255)
                                 gCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCtt                                                                                                                 <  2:2596244/150‑1 (MQ=255)
                                   aCGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCa                                                                                                               <  1:478458/150‑1 (MQ=255)
                                   aCGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCa                                                                                                               <  1:3539844/150‑1 (MQ=255)
                                      gtgtACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATgg                                                                                                            <  2:5147322/150‑1 (MQ=255)
                                        gtACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCg                                                                                                          >  2:3037489/1‑150 (MQ=255)
                                         tACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAACTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCga                                                                                                         >  2:3900167/1‑150 (MQ=255)
                                            tGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatag                                                                                                      >  2:3757592/1‑150 (MQ=255)
                                             gTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgataga                                                                                                     >  1:874725/1‑150 (MQ=255)
                                             gTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgataga                                                                                                     <  1:1355405/150‑1 (MQ=255)
                                              tttCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatagat                                                                                                    <  2:4876570/150‑1 (MQ=255)
                                              tttCGCCAGTAGAGTCAGCGGGGCAAACTCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGGTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatagat                                                                                                    <  2:1538261/150‑1 (MQ=255)
                                               ttCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatagata                                                                                                   >  1:1708047/1‑150 (MQ=255)
                                                tCGCCAGTAGATTCAGCGGGGCAAAATCAACATATTTCGCCAGCTCCTGTAGCACGCGTACCAGCCCTTCATGCTGCCCCGTCCCTCGCATGAAATTCCTTTGAAGTTTCTGCGGAGTCAGGTTTTTAAACTCTTCATGTCGATAGATAc                                                                                                  <  2:1031694/150‑1 (MQ=255)
                                                tCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATAc                                                                                                  <  2:2446960/150‑1 (MQ=255)
                                                tCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATAc                                                                                                  <  1:150998/150‑1 (MQ=255)
                                                 agCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCACTCGCATGAAATACATTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATAcc                                                                                                 >  2:5132243/2‑150 (MQ=255)
                                                  gCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCg                                                                                                >  2:4272119/1‑150 (MQ=255)
                                                   ccAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGc                                                                                               >  1:4669770/1‑150 (MQ=255)
                                                   ccAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGc                                                                                               >  1:2851899/1‑150 (MQ=255)
                                                    cAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCt                                                                                              <  1:555518/150‑1 (MQ=255)
                                                       tAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCt                                                                                           <  2:4450598/150‑1 (MQ=255)
                                                       tAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCt                                                                                           <  1:1956039/150‑1 (MQ=255)
                                                             cAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCttt                                                                                     >  2:2535225/1‑150 (MQ=255)
                                                              aGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCtttt                                                                                    <  2:3255101/150‑1 (MQ=255)
                                                               gCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCttttt                                                                                   <  1:2448940/150‑1 (MQ=255)
                                                               gCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCttttt                                                                                   <  1:5008199/150‑1 (MQ=255)
                                                                  gggCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCCTGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCa                                                                                <  2:3106525/150‑1 (MQ=255)
                                                                   ggCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAg                                                                               <  2:1179047/150‑1 (MQ=255)
                                                                      aaaaTCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCAc                                                                            >  2:631959/1‑150 (MQ=255)
                                                                       aaaTCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACt                                                                           <  2:1402400/150‑1 (MQ=255)
                                                                        aaTCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACtt                                                                          >  2:2322191/1‑150 (MQ=255)
                                                                          tCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCg                                                                        <  2:715480/150‑1 (MQ=255)
                                                                           cAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCgg                                                                       >  2:845246/1‑150 (MQ=255)
                                                                             aCATAGTTCGCCAGCTCCTTACTCACGCGGACCATCCCTTCATGCGGACCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGAGTCATGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGgtg                                                                     <  1:2527141/150‑1 (MQ=255)
                                                                               atcgTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTGTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTc                                                                   <  2:3015806/147‑1 (MQ=255)
                                                                                tAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCg                                                                  <  1:2939507/150‑1 (MQ=255)
                                                                                 aGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGa                                                                 >  1:1310851/1‑150 (MQ=255)
                                                                                    tCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGatat                                                              >  1:2832121/1‑150 (MQ=255)
                                                                                      gCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCa                                                            >  1:471480/1‑150 (MQ=255)
                                                                                        cAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACt                                                          <  1:2040326/150‑1 (MQ=255)
                                                                                        cAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACt                                                          >  2:359609/1‑150 (MQ=255)
                                                                                         aGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTg                                                         <  1:260474/150‑1 (MQ=255)
                                                                                             ccTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCaa                                                     >  1:440456/1‑150 (MQ=255)
                                                                                               tGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAac                                                   <  1:523471/150‑1 (MQ=255)
                                                                                                  cgcacgcGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACt                                                <  1:3070175/150‑1 (MQ=255)
                                                                                                  cgcacgcGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGAGAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACt                                                >  2:4762761/1‑150 (MQ=255)
                                                                                                   gaacgcGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACtt                                               >  2:4632516/3‑150 (MQ=255)
                                                                                                    cacgcGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACttt                                              <  1:4174990/150‑1 (MQ=255)
                                                                                                    cacgcGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACttt                                              <  1:3958009/150‑1 (MQ=255)
                                                                                                         gTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGc                                         <  1:3313409/150‑1 (MQ=255)
                                                                                                          tACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCTCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCg                                        >  1:2812070/1‑150 (MQ=255)
                                                                                                          tACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCg                                        >  1:2002236/1‑150 (MQ=255)
                                                                                                           aCCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGt                                       >  1:4632819/1‑150 (MQ=255)
                                                                                                           aCCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGt                                       <  2:4986989/150‑1 (MQ=255)
                                                                                                              aGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCg                                    <  2:4069801/150‑1 (MQ=255)
                                                                                                              aGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCg                                    <  1:2419261/150‑1 (MQ=255)
                                                                                                                   ttCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGgcgc                               >  2:2210059/1‑150 (MQ=255)
                                                                                                                   ttCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGgcgc                               >  1:424586/1‑150 (MQ=255)
                                                                                                                      gtGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGt                            >  1:165368/2‑150 (MQ=255)
                                                                                                                       tGCGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTg                           <  2:1651149/150‑1 (MQ=255)
                                                                                                                         cGGCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGcc                         <  2:1944452/150‑1 (MQ=255)
                                                                                                                         cGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGACTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGcc                         >  1:2174960/1‑150 (MQ=255)
                                                                                                                          ggCCCCGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGccc                        >  2:2006786/1‑150 (MQ=255)
                                                                                                                             cccGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAt                     <  1:3757592/150‑1 (MQ=255)
                                                                                                                             cccGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCACACACTTTCCAGCGTCCGGGCGCGGTGCCCAAt                     >  2:4379272/1‑150 (MQ=255)
                                                                                                                               tgTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAtgt                   >  1:160280/2‑150 (MQ=255)
                                                                                                                               cGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAtgt                   <  2:1146064/150‑1 (MQ=255)
                                                                                                                               cGTCCCTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAtgt                   >  1:4491169/1‑150 (MQ=255)
                                                                                                                               cGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAAtgt                   >  2:4330309/1‑150 (MQ=255)
                                                                                                                                   ccTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTca               >  2:3963285/1‑150 (MQ=255)
                                                                                                                                    cTCGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcag              >  1:1773853/1‑150 (MQ=255)
                                                                                                                                      cGCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGACAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcagcc            >  2:3526506/1‑150 (MQ=255)
                                                                                                                                       gCATGAAATACCTTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcagcca           <  1:3187990/150‑1 (MQ=255)
                                                                                                                                       gCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTcagcca           >  2:1303755/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGATACCGCTTCTGGCTTTTTCCAGCACTTCGGTGTCGATATCACTGGCAAACACTTTCCAGCGTCCGGGCGCGGTGCCCAATGTGTCAGCCAGCGTCATCG  >  CP009273/1962928‑1963221

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: