Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,184,128 C→T 6.8% S46S (AGC→AGT yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,184,1280CT 6.8% 936.2 / 37.6 266S46S (AGC→AGTyqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (123/125);  new base T (10/8);  total (133/133)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 272 total aligned reads displayed.

AGATCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGACT  >  CP009273/3183980‑3184276
                                                                                                                                                    |                                                                                                                                                    
agaTCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                                                                                                                                                      >  1:3416785/1‑149 (MQ=255)
agaTCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                                                                                                                                                      <  2:3937277/149‑1 (MQ=255)
agaTCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTCTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGc                                                                                                                                                      >  1:803621/1‑149 (MQ=255)
 gaTCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTAAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCg                                                                                                                                                     >  2:3424845/1‑149 (MQ=255)
   tCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGAt                                                                                                                                                   <  2:1239228/149‑1 (MQ=255)
    cAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATg                                                                                                                                                  <  1:328907/149‑1 (MQ=255)
     aataatATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATgg                                                                                                                                                 >  1:641208/1‑149 (MQ=255)
      ataataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGt                                                                                                                                                >  1:789232/1‑149 (MQ=255)
      ataataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGt                                                                                                                                                >  2:557221/1‑149 (MQ=255)
        aataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaa                                                                                                                                              <  2:330221/149‑1 (MQ=255)
        aataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaa                                                                                                                                              >  1:2424963/1‑149 (MQ=255)
             gcgcTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAg                                                                                                                                         <  1:3769124/149‑1 (MQ=255)
               gcTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGt                                                                                                                                       >  1:3503307/1‑149 (MQ=255)
                  aCTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGAt                                                                                                                                    <  1:1994251/149‑1 (MQ=255)
                   cTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATa                                                                                                                                   >  1:791687/1‑149 (MQ=255)
                    ttGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATaa                                                                                                                                  >  1:550524/1‑149 (MQ=255)
                        tAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATc                                                                                                                              <  2:3652006/149‑1 (MQ=255)
                        tAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATc                                                                                                                              <  1:1801870/149‑1 (MQ=255)
                         aaTTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCt                                                                                                                             <  1:3387208/149‑1 (MQ=255)
                           ttGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTgg                                                                                                                           <  1:2228223/149‑1 (MQ=255)
                             gTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCa                                                                                                                         >  1:2583764/1‑149 (MQ=255)
                                 ttaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGc                                                                                                                     >  1:2506864/1‑149 (MQ=255)
                                 ttaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGc                                                                                                                     >  1:4243396/1‑149 (MQ=255)
                                  taatttttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCa                                                                                                                    <  1:3600997/146‑1 (MQ=255)
                                  taCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCa                                                                                                                    <  1:1206607/149‑1 (MQ=255)
                                  taCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCa                                                                                                                    <  1:1654515/149‑1 (MQ=255)
                                   aCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAg                                                                                                                   >  1:269137/1‑149 (MQ=255)
                                     tttttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACGTGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTg                                                                                                                 >  1:4673965/1‑149 (MQ=255)
                                     tttttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTg                                                                                                                 <  1:1856168/149‑1 (MQ=255)
                                        ttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAat                                                                                                              <  1:3888710/149‑1 (MQ=255)
                                         tttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAata                                                                                                             <  1:4370620/149‑1 (MQ=255)
                                           tATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtt                                                                                                           <  1:932491/149‑1 (MQ=255)
                                            aTGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCACAGCAGTGAATAttt                                                                                                          >  1:2269737/1‑149 (MQ=255)
                                            aTGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAttt                                                                                                          >  2:3347732/1‑149 (MQ=255)
                                            aTGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAttt                                                                                                          <  1:2326345/149‑1 (MQ=255)
                                             tGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtttt                                                                                                         >  1:2929211/1‑149 (MQ=255)
                                             tGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtttt                                                                                                         >  1:29974/1‑149 (MQ=255)
                                               ggTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTAt                                                                                                       <  2:355627/149‑1 (MQ=255)
                                                 ttttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGa                                                                                                     >  1:3670618/1‑149 (MQ=255)
                                                 ttttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGa                                                                                                     >  1:2251045/1‑149 (MQ=255)
                                                  tttAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGaa                                                                                                    >  2:2166213/1‑149 (MQ=255)
                                                     aGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGt                                                                                                 >  2:2767480/1‑149 (MQ=255)
                                                      gTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATCTTTTATGAATGtt                                                                                                >  2:4680776/1‑149 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                               <  1:1539355/149‑1 (MQ=255)
                                                         cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTc                                                                                             >  1:4466440/1‑149 (MQ=255)
                                                         cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTc                                                                                             >  2:3368074/1‑149 (MQ=255)
                                                          aTTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCt                                                                                            >  1:467192/1‑149 (MQ=255)
                                                           ttACCTGCATGGTCAATGGACTGCTATTCTGAAAATGAAGGTGTAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCtt                                                                                           <  2:2571909/149‑2 (MQ=255)
                                                           ttACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGTAAATACCGTTGTAATTGGTTACGTACCAAGAATCTCTATCCCAAGCGATGGTAAAAAAGTTGATAAAATATGGAAAAGAAGTGAATATTTTATTAATGttttatt                                                                                           >  2:3521313/1‑146 (MQ=255)
                                                           ttACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCCGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTg                                                                                           <  2:2313012/149‑1 (MQ=255)
                                                           ttACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTg                                                                                           >  2:2222662/1‑149 (MQ=255)
                                                           ttACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTg                                                                                           <  2:1349274/149‑1 (MQ=255)
                                                              ccTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGtaa                                                                                        <  1:3275242/149‑1 (MQ=255)
                                                               cTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGtaat                                                                                       >  2:4070387/1‑149 (MQ=255)
                                                               cTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGtaat                                                                                       >  2:4403344/1‑149 (MQ=255)
                                                                 gCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGtaataa                                                                                     <  2:80116/149‑1 (MQ=255)
                                                                      ggCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCAc                                                                                <  2:1671697/149‑1 (MQ=255)
                                                                         aaTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTAc                                                                             <  2:862326/149‑1 (MQ=255)
                                                                                  cTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAt                                                                    >  2:1032723/1‑149 (MQ=255)
                                                                                   tATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAtc                                                                   <  1:3309946/149‑1 (MQ=255)
                                                                                   tATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCAATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAtc                                                                   <  1:4150215/149‑1 (MQ=255)
                                                                                     tGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAtctc                                                                 >  2:4623820/1‑149 (MQ=255)
                                                                                          aaCATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGga                                                            >  2:1371684/1‑149 (MQ=255)
                                                                                              tGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGagaag                                                        >  1:3097658/1‑149 (MQ=255)
                                                                                               gAAGGTTGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGagaaga                                                       <  2:2882973/149‑1 (MQ=255)
                                                                                                 aGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTTATAATGCACTACCAGGCCCATCTCCAGGAGAAGAAt                                                     <  2:1738722/149‑1 (MQ=255)
                                                                                                 aGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAAt                                                     <  2:3816841/149‑1 (MQ=255)
                                                                                                      gAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCa                                                >  2:1330920/1‑149 (MQ=255)
                                                                                                       aaaTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCAt                                               <  2:2946654/149‑1 (MQ=255)
                                                                                                        aaTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATc                                              <  1:984833/149‑1 (MQ=255)
                                                                                                         aTACCGTTGTCATTGGTTACGTACAAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCt                                             <  2:1841466/149‑1 (MQ=255)
                                                                                                          tACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTg                                            >  1:3137041/1‑149 (MQ=255)
                                                                                                           aCCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGc                                           >  1:1068647/1‑149 (MQ=255)
                                                                                                               ttGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATgg                                       >  1:3544256/1‑149 (MQ=255)
                                                                                                                      tGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAata                                <  2:1394769/149‑1 (MQ=255)
                                                                                                                            cGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatga                          <  1:4291756/149‑1 (MQ=255)
                                                                                                                              tACCAAGAATCTCTATCCCCAGCGATGGTCAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatgatg                        >  2:3935930/1‑149 (MQ=255)
                                                                                                                              tACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatgatg                        <  2:269097/149‑1 (MQ=255)
                                                                                                                               aCCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGAtgt                       >  1:1115685/1‑149 (MQ=255)
                                                                                                                               aCCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGAtgt                       <  1:1371684/149‑1 (MQ=255)
                                                                                                                                ccAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGAtgtt                      >  1:2102517/1‑149 (MQ=255)
                                                                                                                                   aGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATgttgtt                   >  2:4293252/1‑149 (MQ=255)
                                                                                                                                   aGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATgttgtt                   >  2:1556327/1‑149 (MQ=255)
                                                                                                                                   aGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATgttgtt                   >  2:3071358/1‑149 (MQ=255)
                                                                                                                                     aaTCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAg                 >  1:1403015/1‑149 (MQ=255)
                                                                                                                                      aTCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGc                >  2:2356814/1‑149 (MQ=255)
                                                                                                                                       tctcTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCa               >  1:3477627/1‑149 (MQ=255)
                                                                                                                                        ctctATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCCGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCAt              <  1:499128/149‑1 (MQ=255)
                                                                                                                                        ctctATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCAt              <  1:672619/149‑1 (MQ=255)
                                                                                                                                         tctATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATc             <  2:1941558/149‑1 (MQ=255)
                                                                                                                                           tATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAg           >  1:4036434/1‑149 (MQ=255)
                                                                                                                                           tATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAg           <  1:1623934/149‑1 (MQ=255)
                                                                                                                                            aTCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAgg          >  2:239057/1‑149 (MQ=255)
                                                                                                                                              ccccAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTc        <  2:789226/149‑1 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      <  1:148022/149‑1 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      <  1:2827047/149‑1 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      <  1:2827126/149‑1 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      <  2:2072556/149‑1 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      >  2:3177189/1‑149 (MQ=255)
                                                                                                                                                ccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGAAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCaa      >  1:1838177/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGATCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGACT  >  CP009273/3183980‑3184276

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: