Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,469,418 C→A 17.2% R23R (CGA→AGA)  yfdT → CPS‑53 (KpLE1) prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,469,4180CA17.2% 800.3 / 113.5 264R23R (CGA→AGA) yfdTCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (118/100);  new base A (34/12);  total (152/112)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

GACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATAC  >  CP009273/2469271‑2469567
                                                                                                                                                   |                                                                                                                                                     
  cATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCgaga                                                                                                                                                   <  1:3883121/150‑1 (MQ=255)
    tttCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACg                                                                                                                                                 >  1:2213952/1‑150 (MQ=255)
     ttCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGa                                                                                                                                                >  2:3223052/1‑150 (MQ=255)
     ttCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGa                                                                                                                                                <  1:3458812/150‑1 (MQ=255)
     ttCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGa                                                                                                                                                >  2:4234169/1‑150 (MQ=255)
      tCAGTTTGGGGACCTTGAAATCTTCAATATCACTCACTTTGCAGTTATGCCAAAAGCACTTGAGGCTGTATCATAATGACCACATTTACAAATAAAGATCTGATTAAAGAAATAAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGat                                                                                                                                               <  1:3782890/150‑1 (MQ=255)
         ggTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATAtt                                                                                                                                            <  1:2258720/150‑1 (MQ=255)
         ggTTGGGGAAATTGAAAAAGTAAATATCACTCACTGGTAAGTAATGCCAAAAGCACCGGAGGATTGAGCATAATGACCAAATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATAtt                                                                                                                                            <  1:3159158/150‑1 (MQ=255)
          gTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTg                                                                                                                                           >  2:537425/1‑150 (MQ=255)
                       aaaaCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTa                                                                                                                              >  1:3096345/1‑150 (MQ=255)
                       aaaaCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTa                                                                                                                              >  2:4411600/1‑150 (MQ=255)
                          aCGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCTCCGTGCTTATGa                                                                                                                           <  2:3197301/150‑1 (MQ=255)
                           cGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGaa                                                                                                                          <  2:1010478/150‑1 (MQ=255)
                            gTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGaaa                                                                                                                         <  2:4092817/150‑1 (MQ=255)
                            gTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGaaa                                                                                                                         >  2:1254048/1‑150 (MQ=255)
                             tCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAAt                                                                                                                        <  2:1234545/150‑1 (MQ=255)
                               aaTATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCg                                                                                                                      >  1:4239120/1‑150 (MQ=255)
                                 tatCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCa                                                                                                                    >  1:4612581/1‑150 (MQ=255)
                                   tcactcacTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACt                                                                                                                  <  2:3043624/150‑1 (MQ=255)
                                    cactcactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTc                                                                                                                 <  1:3393106/150‑1 (MQ=255)
                                     actcactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCg                                                                                                                >  1:2992817/1‑150 (MQ=255)
                                        cactGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCAAATTTACCAATAAATAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTAt                                                                                                             <  2:2847793/150‑1 (MQ=255)
                                          ctGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTAtct                                                                                                           <  1:2058613/150‑1 (MQ=255)
                                                aGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGaa                                                                                                     >  1:3027045/1‑150 (MQ=255)
                                                   aaTGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTa                                                                                                  <  2:4206084/150‑1 (MQ=255)
                                                         aaaaGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCa                                                                                            >  2:313238/1‑150 (MQ=255)
                                                         aaaaGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCa                                                                                            >  1:2608099/1‑150 (MQ=255)
                                                           aaGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGa                                                                                          >  1:2377391/1‑150 (MQ=255)
                                                             gCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATg                                                                                        <  1:2511399/150‑1 (MQ=255)
                                                             gCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATg                                                                                        >  2:3201409/1‑150 (MQ=255)
                                                              cACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGa                                                                                       <  1:4639498/150‑1 (MQ=255)
                                                              cACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGa                                                                                       >  1:1787547/1‑150 (MQ=255)
                                                               aCCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGaa                                                                                      <  1:2575105/150‑1 (MQ=255)
                                                                   gAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAAcgcg                                                                                  >  1:3586736/1‑150 (MQ=255)
                                                                   gAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAAcgcg                                                                                  <  1:1650533/150‑1 (MQ=255)
                                                                   gAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAAcgcg                                                                                  <  1:1644514/150‑1 (MQ=255)
                                                                   gAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAAcgcg                                                                                  >  1:3281143/1‑150 (MQ=255)
                                                                    aGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGa                                                                                 >  2:3489857/1‑150 (MQ=255)
                                                                        tGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGcc                                                                             >  1:2081669/1‑150 (MQ=255)
                                                                        tGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGcc                                                                             <  1:1566172/150‑1 (MQ=255)
                                                                           aGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTAt                                                                          <  1:980466/150‑1 (MQ=255)
                                                                           aGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTAt                                                                          <  1:322241/150‑1 (MQ=255)
                                                                           aGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTAt                                                                          <  1:338539/150‑1 (MQ=255)
                                                                                aaTGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAatt                                                                     <  1:443488/150‑1 (MQ=255)
                                                                                 aTGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAatta                                                                    >  2:2910846/1‑150 (MQ=255)
                                                                                 aTGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAatta                                                                    >  2:4632334/1‑150 (MQ=255)
                                                                                  tGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAattat                                                                   <  1:4325905/150‑1 (MQ=255)
                                                                                  tGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAattat                                                                   <  2:4541869/150‑1 (MQ=255)
                                                                                  tGACCACATTTACCAATAAAGAACAGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCTTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAattat                                                                   >  1:1950006/1‑150 (MQ=255)
                                                                                   gACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAattatt                                                                  <  1:3155372/150‑1 (MQ=255)
                                                                                     ccACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCa                                                                >  2:2435967/1‑150 (MQ=255)
                                                                                        caTTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATgg                                                             <  1:4309124/150‑1 (MQ=255)
                                                                                            tACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGaaaa                                                         <  1:1960233/150‑1 (MQ=255)
                                                                                             aCCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAg                                                        <  1:4808455/150‑1 (MQ=255)
                                                                                                aaTAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGt                                                     >  1:3819099/1‑150 (MQ=255)
                                                                                                aaTAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGt                                                     >  1:2153247/1‑150 (MQ=255)
                                                                                                 aTAAAGAACTGATTAAAGAAATCAAAGACCGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGtt                                                    <  2:4727716/150‑1 (MQ=255)
                                                                                                 aTAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGACTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCAAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGtt                                                    <  1:3328265/150‑1 (MQ=255)
                                                                                                  tAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGttt                                                   >  1:469093/1‑150 (MQ=255)
                                                                                                  tAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGttt                                                   >  1:629832/1‑150 (MQ=255)
                                                                                                    aaGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCg                                                 >  2:4500917/1‑150 (MQ=255)
                                                                                                      gAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGt                                               >  1:3956656/1‑150 (MQ=255)
                                                                                                       aaCTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTg                                              >  1:514382/1‑150 (MQ=255)
                                                                                                          tGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACt                                           >  2:1796470/1‑150 (MQ=255)
                                                                                                          tGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACt                                           <  2:1582998/150‑1 (MQ=255)
                                                                                                          tGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACt                                           <  2:1310181/150‑1 (MQ=255)
                                                                                                            aTTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTta                                         <  1:1359574/150‑1 (MQ=255)
                                                                                                             ttAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTtag                                        >  1:759034/1‑150 (MQ=255)
                                                                                                                 aGAACTCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTAGGAACTCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGa                                    <  2:1760480/150‑1 (MQ=255)
                                                                                                                  gAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGAt                                   >  1:3130287/1‑150 (MQ=255)
                                                                                                                   aaaTCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGAtc                                  >  2:2194324/1‑150 (MQ=255)
                                                                                                                   aaaTCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGAtc                                  >  2:170709/1‑150 (MQ=255)
                                                                                                                   aaaTCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGAtc                                  >  2:2602588/1‑150 (MQ=255)
                                                                                                                        aaaGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCgg                             >  2:4145492/1‑150 (MQ=255)
                                                                                                                         aaGAACGACTCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGga                            <  2:832493/150‑1 (MQ=255)
                                                                                                                         aaGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGga                            <  2:832532/150‑1 (MQ=255)
                                                                                                                          aGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTTATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGgag                           >  2:3168684/1‑150 (MQ=255)
                                                                                                                           gAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGgaga                          <  1:254372/150‑1 (MQ=255)
                                                                                                                           gAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGgaga                          <  1:4004925/150‑1 (MQ=255)
                                                                                                                              cGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCAAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGAGACTTAGTAGATCGTCGGAGAGCa                       >  2:88784/1‑150 (MQ=255)
                                                                                                                              cGAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCa                       >  1:290220/1‑150 (MQ=255)
                                                                                                                               gAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaa                      >  1:2817881/1‑150 (MQ=255)
                                                                                                                               gAATCAGCAGCCTAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaa                      >  2:3546364/1‑150 (MQ=255)
                                                                                                                                aaTCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaaa                     >  1:1168147/1‑150 (MQ=255)
                                                                                                                                  tCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaaaaa                   >  2:4409997/1‑150 (MQ=255)
                                                                                                                                   cagcagCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaaaaaa                  >  1:4508231/1‑150 (MQ=255)
                                                                                                                                   cagcagCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCaaaaaa                  >  2:2169831/1‑150 (MQ=255)
                                                                                                                                    agcagcCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAAc                 >  1:2303892/1‑150 (MQ=255)
                                                                                                                                      cagcCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACgg               >  2:715917/1‑150 (MQ=255)
                                                                                                                                       agcCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGc              <  2:661663/150‑1 (MQ=255)
                                                                                                                                           tAGAGGTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATa          >  2:3844809/1‑150 (MQ=255)
                                                                                                                                             gtggTTCGTGTCTATATTGAGCGCGGTGATTATTAAATCGCACTCTTATCTCTGGATGTAGAGCCAGATGAACGCGAAGCCTATTAATTTTTCTTGTAAATGGTTTTCGGTTACTTTGTAGATCGTCGGAGAGAAAAAAACGGCGATAAc        <  2:3314653/148‑1 (MQ=255)
                                                                                                                                               ggTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGa      >  2:970540/1‑150 (MQ=255)
                                                                                                                                               ggTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGa      >  2:2057475/1‑150 (MQ=255)
                                                                                                                                               ggTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGa      >  1:4556903/1‑150 (MQ=255)
                                                                                                                                               ggTTAGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGa      >  1:27581/1‑150 (MQ=255)
                                                                                                                                                gTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGaa     >  2:3665855/1‑150 (MQ=255)
                                                                                                                                                 ttCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAAt    >  2:2334439/1‑150 (MQ=255)
                                                                                                                                                   cGAGACGTTATTGAGCGCCTTGCTTATGAAATCGCACTCTTATCTCTGGAAGTCGAGCCAGTTGAACGCGAAGCCTATTAATTATTCATGTAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATAc  <  2:2964203/150‑1 (MQ=255)
                                                                                                                                                   cGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATAc  >  1:4840811/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACATTTCAGGTTGGGGACCTTGAAAACGTCAATATCACTCACTGGGCAGTAATGCCAAAAGCACCGGAGGCTGGAGCATAATGACCACATTTACCAATAAAGAACTGATTAAAGAAATCAAAGAACGAATCAGCAGCCTAGAGGTTCGAGACGATATTGAGCGCCGTGCTTATGAAATCGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATAC  >  CP009273/2469271‑2469567

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: