Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,771,366 G→A 17.3% intergenic (+225/+139) ypjF → / ← ypjA CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,771,3660GA17.3% 726.7 / 121.2 273intergenic (+225/+139)ypjF/ypjACP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base G (124/100);  new base A (28/19);  total (153/120)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 274 total aligned reads displayed.

CTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAC  >  CP009273/2771219‑2771514
                                                                                                                                                   |                                                                                                                                                    
ctccccCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACACGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGac                                                                                                                                                    <  2:2071553/150‑1 (MQ=255)
 tccccCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCGGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGaca                                                                                                                                                   >  2:3760754/1‑150 (MQ=255)
    cccGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGacaaca                                                                                                                                                <  2:2901333/150‑1 (MQ=255)
     ccGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAacac                                                                                                                                               >  2:4622352/1‑150 (MQ=255)
      cGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACAcc                                                                                                                                              <  2:1277843/150‑1 (MQ=255)
       gAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACca                                                                                                                                             >  2:1948424/1‑150 (MQ=255)
        aaaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTACGGACAACACcaa                                                                                                                                            >  2:2475113/1‑150 (MQ=255)
         aaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACcaac                                                                                                                                           >  2:1779099/1‑150 (MQ=255)
         aaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACcaac                                                                                                                                           <  2:3128021/150‑1 (MQ=255)
         aaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACcaac                                                                                                                                           <  2:1982136/150‑1 (MQ=255)
            aaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACaaa                                                                                                                                        >  1:1510342/1‑150 (MQ=255)
            aaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAAATTAACAACAAGTTACACAAACAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACaaa                                                                                                                                        <  2:766536/150‑1 (MQ=255)
              ccTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGCCAACACCAACaaaaa                                                                                                                                      >  2:1407263/1‑150 (MQ=255)
                 gACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAg                                                                                                                                   <  1:2859461/150‑1 (MQ=255)
                    aTTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGaa                                                                                                                                <  2:4224144/150‑1 (MQ=255)
                    aTTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACTAAAACAGGaa                                                                                                                                <  2:750655/150‑1 (MQ=255)
                      tttCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGt                                                                                                                              >  1:4249392/1‑150 (MQ=255)
                         cTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAg                                                                                                                           <  1:3543550/150‑1 (MQ=255)
                             tAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGt                                                                                                                       <  1:1053084/150‑1 (MQ=255)
                              aGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGtc                                                                                                                      <  2:3009647/150‑1 (MQ=255)
                               ggACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGtct                                                                                                                     <  1:3175336/150‑1 (MQ=255)
                                   caacaaTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGc                                                                                                                 <  1:3811593/150‑1 (MQ=255)
                                        aTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAac                                                                                                            >  2:1632158/1‑150 (MQ=255)
                                           ggACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACAcc                                                                                                         <  1:1204120/150‑1 (MQ=255)
                                            gACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCg                                                                                                        >  2:233368/1‑150 (MQ=255)
                                             aCCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGa                                                                                                       >  1:963383/1‑150 (MQ=255)
                                               cAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGAcc                                                                                                     >  1:2531755/1‑150 (MQ=255)
                                                aaaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCa                                                                                                    >  2:4090773/1‑150 (MQ=255)
                                                aaaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCa                                                                                                    >  1:3745529/1‑150 (MQ=255)
                                                aaaaTGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCa                                                                                                    >  2:3572571/1‑150 (MQ=255)
                                                   aTGAAAAAAGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATACAGACTCCGCCAGCCCACCAATCATGAATGGACAGTGAAACTACAACACAAACAAAAACAGGAAGTTATAAGTCTCAGCAAAACACCGAACAGAc                                                                                                 >  2:4072961/1‑150 (MQ=255)
                                                    tGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACg                                                                                                <  1:3870156/150‑1 (MQ=255)
                                                    tGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACg                                                                                                >  2:1196195/1‑150 (MQ=255)
                                                    tGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACAAAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACg                                                                                                >  2:4281140/1‑150 (MQ=255)
                                                      aaaaTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACCGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGt                                                                                              <  2:1930877/150‑1 (MQ=255)
                                                         aTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgag                                                                                           >  1:1896300/1‑150 (MQ=255)
                                                         aTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgag                                                                                           >  1:1896366/1‑150 (MQ=255)
                                                         aTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgag                                                                                           <  2:683657/150‑1 (MQ=255)
                                                         aTTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgag                                                                                           >  1:4428103/1‑150 (MQ=255)
                                                           tGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgagga                                                                                         <  1:630724/150‑1 (MQ=255)
                                                            gAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTgaggag                                                                                        <  1:855709/150‑1 (MQ=255)
                                                             aaCTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGgaga                                                                                       >  2:1175084/1‑150 (MQ=255)
                                                              aCTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGAc                                                                                      <  2:2827816/150‑1 (MQ=255)
                                                              aCTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGAc                                                                                      <  2:2215720/150‑1 (MQ=255)
                                                               cTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACa                                                                                     >  2:2955478/1‑150 (MQ=255)
                                                                 gAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCCACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATa                                                                                   >  2:4346014/1‑150 (MQ=255)
                                                                   aGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATaaa                                                                                 >  1:4038138/1‑150 (MQ=255)
                                                                   aGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATaaa                                                                                 >  1:3902056/1‑150 (MQ=255)
                                                                       tACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGa                                                                             >  1:2132830/1‑150 (MQ=255)
                                                                           atcaatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACCTAAAAGGATACg                                                                         >  2:3500263/1‑150 (MQ=255)
                                                                              aatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCaa                                                                      <  2:61115/150‑1 (MQ=255)
                                                                              aatTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACACCACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCaa                                                                      >  2:1751935/1‑150 (MQ=255)
                                                                               atTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCaaa                                                                     >  1:3046679/1‑150 (MQ=255)
                                                                                 ttAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAgg                                                                   >  1:1489008/1‑150 (MQ=255)
                                                                                  tAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGa                                                                  >  1:2096763/1‑150 (MQ=255)
                                                                                   aacaacAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAg                                                                 >  1:3423890/1‑150 (MQ=255)
                                                                                    acaacaAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACATACACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGc                                                                >  1:2717974/1‑150 (MQ=255)
                                                                                    acaacaAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGc                                                                <  1:2517448/150‑1 (MQ=255)
                                                                                    acaacaAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGc                                                                <  1:1384138/150‑1 (MQ=255)
                                                                                        caaGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCcgcg                                                            <  2:987956/150‑1 (MQ=255)
                                                                                          aGTTACACACCCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAACAGGATACGCAAAGGAGCCGCGGc                                                          <  2:1510323/150‑1 (MQ=255)
                                                                                          aGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGc                                                          >  1:2075504/1‑150 (MQ=255)
                                                                                           gTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCt                                                         >  1:2704865/1‑150 (MQ=255)
                                                                                           gTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCt                                                         >  2:1391924/1‑150 (MQ=255)
                                                                                               cacaACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTg                                                     >  1:2771377/1‑150 (MQ=255)
                                                                                                acaACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTgg                                                    >  1:4216061/1‑150 (MQ=255)
                                                                                                  aaCCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTa                                                  <  2:1664981/150‑1 (MQ=255)
                                                                                                  aaCCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTa                                                  >  2:1867941/1‑150 (MQ=255)
                                                                                                      aaTTCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCACCAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCCGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACAt                                              <  2:3943922/150‑1 (MQ=255)
                                                                                                         tCAGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGaa                                           >  1:2366630/1‑150 (MQ=255)
                                                                                                         tCAGACTCCGCCAGCCCACCAATCATGATGGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGaa                                           <  2:3032707/150‑1 (MQ=255)
                                                                                                           aGACTCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAg                                         >  1:2214689/1‑150 (MQ=255)
                                                                                                               tCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCca                                     >  2:4161968/1‑150 (MQ=255)
                                                                                                               tCCGCCAGCCCACCAATCATGATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCca                                     >  2:3290428/1‑150 (MQ=255)
                                                                                                                ccGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCcac                                    <  2:1530950/150‑1 (MQ=255)
                                                                                                                 cGCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCcaca                                   <  1:2480289/150‑1 (MQ=255)
                                                                                                                  gCCAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAg                                  <  2:3273603/150‑1 (MQ=255)
                                                                                                                   ccAGACCACCAATCATGATTGGACAGTGTAAGTACAACACCAACAAAAACAGGAATTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCAGCGGCTCCTGGTAACATGAAAGCCCACAGa                                 <  2:2155147/150‑1 (MQ=255)
                                                                                                                    cAGCCCACCAATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAt                                <  1:4178963/150‑1 (MQ=255)
                                                                                                                        ccaccaATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                            >  1:4515965/1‑150 (MQ=255)
                                                                                                                        ccaccaATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                            >  1:1982434/1‑150 (MQ=255)
                                                                                                                        ccaccaATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                            <  1:2653459/150‑1 (MQ=255)
                                                                                                                            caATCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCtt                        >  1:1882409/1‑150 (MQ=255)
                                                                                                                              aTCATGATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCtttt                      >  2:3583062/1‑150 (MQ=255)
                                                                                                                               tCATGAGTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttttt                     <  1:4621939/150‑1 (MQ=255)
                                                                                                                                   gATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGtt                 >  2:5063/1‑150 (MQ=255)
                                                                                                                                   gATTGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGtt                 <  2:767099/150‑1 (MQ=255)
                                                                                                                                   gATTGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGtt                 >  2:3646751/1‑150 (MQ=255)
                                                                                                                                     ttGGACAGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGa               <  2:2338807/150‑1 (MQ=255)
                                                                                                                                      tGGACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGAt              <  2:1349676/150‑1 (MQ=255)
                                                                                                                                       ggACAGTGTAAGGACAACACCAACAAAAACATGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCAGGTAACAAGAAAGCCCACAGATGTGGTCTTTTTCGTTTATg             >  1:629313/1‑150 (MQ=255)
                                                                                                                                       ggACAGTGTAAGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg             <  1:3464354/150‑1 (MQ=255)
                                                                                                                                            gtgtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAg        >  2:4186228/1‑150 (MQ=255)
                                                                                                                                            gtgtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAg        <  2:2326413/150‑1 (MQ=255)
                                                                                                                                              gtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGaa      >  1:3341726/1‑150 (MQ=255)
                                                                                                                                              gtAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGaa      >  2:2826278/1‑150 (MQ=255)
                                                                                                                                               tAAGGACAACACCACCAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAAc     <  1:574487/150‑1 (MQ=255)
                                                                                                                                                aaGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACg    >  1:3375235/1‑150 (MQ=255)
                                                                                                                                                aaGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACg    >  2:3280519/1‑150 (MQ=255)
                                                                                                                                                aaGAACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCTTTGATGGTCAGAACg    >  1:1772142/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCCCCCGAAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAC  >  CP009273/2771219‑2771514

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: