Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,858,963 A→T 48.7% V224E (GTA→GAA)  yidE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,858,9630AT48.7% 14.5 / inf 255V224E (GTA→GAA) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base A (55/75);  new base T (64/60);  total (119/136)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 258 total aligned reads displayed.

CACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858815‑3859112
                                                                                                                                                    |                                                                                                                                                     
cACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTAc                                                                                                                                                      >  1:2681954/1‑150 (MQ=255)
 aCCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCa                                                                                                                                                     >  1:3756966/1‑150 (MQ=255)
 aCCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCa                                                                                                                                                     <  2:4221558/150‑1 (MQ=255)
 aCCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                                                                                                                     <  1:1081302/150‑1 (MQ=255)
   cAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATc                                                                                                                                                   >  1:2369169/1‑150 (MQ=255)
      ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCttt                                                                                                                                                <  1:727379/150‑1 (MQ=255)
      ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCttt                                                                                                                                                >  2:4015169/1‑150 (MQ=255)
      ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCttt                                                                                                                                                <  1:1936877/150‑1 (MQ=255)
      ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                                                                <  1:907509/150‑1 (MQ=255)
        gcgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTaa                                                                                                                                              <  1:4122533/150‑1 (MQ=255)
        gcgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTaa                                                                                                                                              <  1:1852669/150‑1 (MQ=255)
             ttATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCt                                                                                                                                         <  2:3294173/150‑1 (MQ=255)
              tATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTa                                                                                                                                        <  1:2547271/150‑1 (MQ=255)
              tATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTa                                                                                                                                        <  1:2547189/150‑1 (MQ=255)
               aTGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaa                                                                                                                                       >  1:1765615/1‑150 (MQ=255)
                 gTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAAt                                                                                                                                     >  2:4432405/1‑150 (MQ=255)
                 gTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAAt                                                                                                                                     <  1:3523310/150‑1 (MQ=255)
                    aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                                                                                                                  >  2:2239741/1‑150 (MQ=255)
                       ccGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCa                                                                                                                               <  1:1003646/150‑1 (MQ=255)
                              ctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAgg                                                                                                                        >  1:459765/1‑150 (MQ=255)
                                gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      >  1:2415292/1‑150 (MQ=255)
                                gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      >  1:1988041/1‑150 (MQ=255)
                                gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      >  1:375066/1‑150 (MQ=255)
                                  cccACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtc                                                                                                                    <  2:2598718/150‑1 (MQ=255)
                                   ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                   >  2:3586482/1‑150 (MQ=255)
                                   ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                   >  2:3251213/1‑150 (MQ=255)
                                    caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                  <  2:551695/150‑1 (MQ=255)
                                     accaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCa                                                                                                                 >  2:4285630/1‑150 (MQ=255)
                                      ccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCaa                                                                                                                <  2:2369033/150‑1 (MQ=255)
                                      ccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCGTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCaa                                                                                                                <  2:2369624/150‑1 (MQ=255)
                                       caGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                               >  2:1076148/1‑150 (MQ=255)
                                       caGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                               <  2:1131348/150‑1 (MQ=255)
                                       caGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                               >  1:2269606/1‑150 (MQ=255)
                                           tGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGa                                                                                                           <  1:1381019/150‑1 (MQ=255)
                                             cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAat                                                                                                         >  1:237664/1‑150 (MQ=255)
                                             cagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAat                                                                                                         >  1:237645/1‑150 (MQ=255)
                                              agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                        <  1:3644126/150‑1 (MQ=255)
                                              agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                        >  1:4488260/1‑150 (MQ=255)
                                              agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                        >  1:2566688/1‑150 (MQ=255)
                                              agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                        >  1:1234500/1‑150 (MQ=255)
                                                  aaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa                                                                                                    >  2:3788738/1‑150 (MQ=255)
                                                    tCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATa                                                                                                  <  2:675938/150‑1 (MQ=255)
                                                       cccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                                                               >  2:1242026/1‑150 (MQ=255)
                                                       cccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                                                               <  1:3433445/150‑1 (MQ=255)
                                                         cAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt                                                                                             <  1:1895896/150‑1 (MQ=255)
                                                          aGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTg                                                                                            <  1:2463385/150‑1 (MQ=255)
                                                          aGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTg                                                                                            <  1:2463351/150‑1 (MQ=255)
                                                          aGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTg                                                                                            <  1:2463115/150‑1 (MQ=255)
                                                            ttGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGAt                                                                                          >  2:81090/1‑150 (MQ=255)
                                                            ttGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGAt                                                                                          <  1:3117195/150‑1 (MQ=255)
                                                               gATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAg                                                                                       >  1:3271728/1‑150 (MQ=255)
                                                               gATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAg                                                                                       >  1:3271410/1‑150 (MQ=255)
                                                                aTAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgc                                                                                      <  1:3983970/150‑1 (MQ=255)
                                                                     ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                 >  2:3772963/1‑150 (MQ=255)
                                                                     ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                 <  1:3789503/150‑1 (MQ=255)
                                                                     ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                 <  2:2292166/150‑1 (MQ=255)
                                                                          cTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCTCCGCCag                                                                            >  2:3516640/1‑149 (MQ=255)
                                                                           tGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCAtt                                                                           <  2:3822473/150‑1 (MQ=255)
                                                                            ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTg                                                                          >  1:860286/1‑150 (MQ=255)
                                                                            ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTg                                                                          >  2:4488106/1‑150 (MQ=255)
                                                                            ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTg                                                                          <  2:1407071/150‑1 (MQ=255)
                                                                             gCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTgg                                                                         >  1:2965975/1‑150 (MQ=255)
                                                                                aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGc                                                                      >  1:2292095/1‑150 (MQ=255)
                                                                                aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGc                                                                      >  1:720456/1‑150 (MQ=255)
                                                                                    aaCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGa                                                                  <  2:1211075/150‑1 (MQ=255)
                                                                                       ttttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATCGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                                                               >  1:2972141/1‑150 (MQ=255)
                                                                                       ttttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                                                               >  1:557278/1‑150 (MQ=255)
                                                                                           aGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTCCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCg                                                           <  2:2419973/150‑1 (MQ=255)
                                                                                            gggTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGt                                                          >  2:4461288/1‑150 (MQ=255)
                                                                                             ggTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtg                                                         <  2:113644/150‑1 (MQ=255)
                                                                                              gTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgc                                                        >  1:2742024/1‑150 (MQ=255)
                                                                                               tttCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgct                                                       >  1:1137039/1‑150 (MQ=255)
                                                                                                   ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   >  2:2962476/1‑150 (MQ=255)
                                                                                                   ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   <  2:2760375/150‑1 (MQ=255)
                                                                                                        gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt                                              <  1:597832/150‑1 (MQ=255)
                                                                                                        gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt                                              >  2:1442810/1‑150 (MQ=255)
                                                                                                        gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt                                              >  2:2783335/1‑150 (MQ=255)
                                                                                                          ttCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCt                                            <  2:1617716/150‑1 (MQ=255)
                                                                                                          ttCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCt                                            <  2:2928757/150‑1 (MQ=255)
                                                                                                            cAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGt                                          <  2:1151520/150‑1 (MQ=255)
                                                                                                             aGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtc                                         <  2:3644168/150‑1 (MQ=255)
                                                                                                              gACGCGAGCAGATAATTTTGTCGCCGTTGAGACTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtct                                        <  2:708767/150‑1 (MQ=255)
                                                                                                              gACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtct                                        <  1:1573990/150‑1 (MQ=255)
                                                                                                              gACGCGAGCAGATAATTGTGTCGCCGTTGAGAATCGGTACATCTTGAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTGGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtct                                        <  2:1997563/150‑1 (MQ=255)
                                                                                                               aCGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtctc                                       >  2:3925983/1‑150 (MQ=255)
                                                                                                               aCGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtctc                                       <  2:4351408/150‑1 (MQ=255)
                                                                                                                cgcgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCg                                      >  1:2492385/1‑150 (MQ=255)
                                                                                                                  cgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGAc                                    <  1:1468158/150‑1 (MQ=255)
                                                                                                                   gAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACa                                   >  1:1054037/1‑150 (MQ=255)
                                                                                                                    aGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACAt                                  <  2:1429795/150‑1 (MQ=255)
                                                                                                                        gATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAc                              <  1:3199576/150‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGaa                        <  1:1723907/150‑1 (MQ=255)
                                                                                                                                        ttGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGc              <  1:4335976/150‑1 (MQ=255)
                                                                                                                                        ttGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGc              <  1:4336374/150‑1 (MQ=255)
                                                                                                                                           agaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc           <  1:2334135/150‑1 (MQ=255)
                                                                                                                                             aaTCGGTACATCTTTAATGGCTAAATCATGCAGGTGAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACAt         <  2:2369192/150‑1 (MQ=255)
                                                                                                                                               tCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATcc       <  1:4488080/150‑1 (MQ=255)
                                                                                                                                                cGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCca      >  2:312150/1‑150 (MQ=255)
                                                                                                                                                   tACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAt   >  1:1962090/1‑150 (MQ=255)
                                                                                                                                                    aCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATg  >  1:452710/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858815‑3859112

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: