Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 518,820 C→A 88.7% A35E (GCG→GAG)  rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273518,8200CA88.7% 637.0 / 76.5 246A35E (GCG→GAG) rhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base A (122/96);  total (134/112)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGAC  >  CP009273/518674‑518964
                                                                                                                                                  |                                                                                                                                                
tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggagt                                                                                                                                                <  1:4383557/149‑1 (MQ=255)
tgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggagt                                                                                                                                                <  1:4588915/149‑1 (MQ=255)
  cTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGc                                                                                                                                              >  1:1338797/1‑149 (MQ=255)
   tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCt                                                                                                                                             <  1:4164266/149‑1 (MQ=255)
   tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCt                                                                                                                                             >  1:602902/1‑149 (MQ=255)
         aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtgt                                                                                                                                       <  2:3095595/149‑1 (MQ=255)
         aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgt                                                                                                                                       >  1:426344/1‑149 (MQ=255)
         aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgt                                                                                                                                       >  2:2026877/1‑149 (MQ=255)
          gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtgtg                                                                                                                                      <  2:1296777/149‑1 (MQ=255)
          gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtg                                                                                                                                      >  2:3066583/1‑149 (MQ=255)
           ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtgt                                                                                                                                     >  1:1325572/1‑149 (MQ=255)
           ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtgt                                                                                                                                     >  1:683979/1‑149 (MQ=255)
           ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtgt                                                                                                                                     >  2:4140719/1‑149 (MQ=255)
           ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtgt                                                                                                                                     >  2:1416413/1‑149 (MQ=255)
            gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTc                                                                                                                                    <  2:3832866/149‑1 (MQ=255)
               atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgg                                                                                                                                 >  2:680020/1‑149 (MQ=255)
               atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgg                                                                                                                                 >  2:268427/1‑149 (MQ=255)
                tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcggg                                                                                                                                >  2:4728742/1‑149 (MQ=255)
                 gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggc                                                                                                                               <  2:3825096/149‑1 (MQ=255)
                   aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgggcgg                                                                                                                             >  1:146326/1‑149 (MQ=255)
                     aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGa                                                                                                                           <  1:2329828/149‑1 (MQ=255)
                        aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGa                                                                                                                        >  2:2420575/1‑149 (MQ=255)
                          cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACt                                                                                                                      <  1:502969/149‑1 (MQ=255)
                           tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACtt                                                                                                                     >  1:1668428/1‑149 (MQ=255)
                            aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTc                                                                                                                    >  2:1594349/1‑149 (MQ=255)
                            aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTc                                                                                                                    >  2:4218928/1‑149 (MQ=255)
                              aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                                                  >  1:802182/1‑149 (MQ=255)
                                 ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                                               >  2:2520988/1‑149 (MQ=255)
                                        tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGt                                                                                                        <  2:2778817/149‑1 (MQ=255)
                                          tatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                      >  1:4299399/1‑149 (MQ=255)
                                          tatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                      >  1:3136629/1‑149 (MQ=255)
                                          tatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                      >  2:1844214/1‑149 (MQ=255)
                                           atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaaa                                                                                                     >  2:2698391/1‑149 (MQ=255)
                                            tGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAAt                                                                                                    <  1:1494786/149‑1 (MQ=255)
                                             gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATc                                                                                                   >  1:570398/1‑149 (MQ=255)
                                             gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATc                                                                                                   <  1:1167756/149‑1 (MQ=255)
                                             gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATc                                                                                                   >  1:2165306/1‑149 (MQ=255)
                                               gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCg                                                                                                 >  2:2057931/1‑149 (MQ=255)
                                                cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTTCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgc                                                                                                <  2:2849046/149‑1 (MQ=255)
                                                cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgc                                                                                                <  2:4490581/149‑1 (MQ=255)
                                                 ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgct                                                                                               <  1:3614899/149‑1 (MQ=255)
                                                  gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctg                                                                                              >  2:552340/1‑149 (MQ=255)
                                                    aaaCCCGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAGTGTCCCGGGTTCGGCAGGTGTCAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGGACATCCgctgct                                                                                            >  2:3287490/1‑149 (MQ=255)
                                                    aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                                                                                            <  1:4513449/149‑1 (MQ=255)
                                                      aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgg                                                                                          >  2:1609833/1‑149 (MQ=255)
                                                             gcgcGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGaa                                                                                   <  1:1942861/149‑1 (MQ=255)
                                                              cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                                                  <  1:1910554/149‑1 (MQ=255)
                                                                cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                                                <  2:614514/149‑1 (MQ=255)
                                                                  tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgc                                                                              <  1:3513396/149‑1 (MQ=255)
                                                                  tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgc                                                                              <  1:3290935/149‑1 (MQ=255)
                                                                   cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                                             >  2:506549/1‑149 (MQ=255)
                                                                    aGGGAGATATGACTCAGTATGGCGGTCCCAGTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                                                            <  1:866425/149‑1 (MQ=255)
                                                                          atatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGc                                                                      >  2:4301506/1‑149 (MQ=255)
                                                                           tatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTCAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg                                                                     >  2:2301066/1‑149 (MQ=255)
                                                                           tatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg                                                                     >  2:3637775/1‑149 (MQ=255)
                                                                            atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                                                                    <  2:4505193/149‑1 (MQ=255)
                                                                            atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga                                                                    <  2:3415533/149‑1 (MQ=255)
                                                                              gACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                                                                  >  1:2075703/1‑149 (MQ=255)
                                                                              gACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga                                                                  <  1:2256718/149‑1 (MQ=255)
                                                                              gACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGCGCTGCCCGGCgaga                                                                  <  2:3533150/149‑1 (MQ=255)
                                                                                 tCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                                               <  2:754187/149‑1 (MQ=255)
                                                                                 tCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACgg                                                               >  1:1823960/1‑149 (MQ=255)
                                                                                        ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                                                        >  1:269521/1‑149 (MQ=255)
                                                                                          cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                                                      <  2:1036008/149‑1 (MQ=255)
                                                                                           ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCt                                                     <  2:1199036/149‑1 (MQ=255)
                                                                                            gTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                                                    >  2:3326357/1‑149 (MQ=255)
                                                                                              cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcc                                                  <  1:1704796/149‑1 (MQ=255)
                                                                                              cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcc                                                  >  1:2091690/1‑149 (MQ=255)
                                                                                                  ttGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggc                                              <  1:2283648/149‑1 (MQ=255)
                                                                                                  ttGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggc                                              <  2:2165306/149‑1 (MQ=255)
                                                                                                   tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                                             <  1:4204698/149‑1 (MQ=255)
                                                                                                   tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcc                                             <  1:235505/149‑1 (MQ=255)
                                                                                                    gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                                            <  1:3720096/149‑1 (MQ=255)
                                                                                                    gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                                            <  1:667623/149‑1 (MQ=255)
                                                                                                     tCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                           >  1:4696468/1‑149 (MQ=255)
                                                                                                       cAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                                         >  1:2761142/1‑149 (MQ=255)
                                                                                                        aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTg                                        >  2:227428/1‑149 (MQ=255)
                                                                                                         ggggTCGGCAGGTGTACGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                       <  2:365284/149‑1 (MQ=255)
                                                                                                         gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                       <  1:298059/149‑1 (MQ=255)
                                                                                                         gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                       >  2:3600504/1‑149 (MQ=255)
                                                                                                         gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                       >  1:4221054/1‑149 (MQ=255)
                                                                                                          ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc                                      >  1:3632253/1‑149 (MQ=255)
                                                                                                          ggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc                                      <  1:4451000/149‑1 (MQ=255)
                                                                                                           gTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg                                     >  2:3029845/1‑149 (MQ=255)
                                                                                                            ttCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGt                                    >  1:3748029/1‑149 (MQ=255)
                                                                                                                gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt                                <  2:3669342/149‑1 (MQ=255)
                                                                                                                gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt                                >  1:871129/1‑149 (MQ=255)
                                                                                                                  aGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc                              >  1:1607175/1‑149 (MQ=255)
                                                                                                                   ggTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGttctttct                             >  2:1999750/1‑149 (MQ=255)
                                                                                                                       tAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTcc                         >  1:3919927/1‑149 (MQ=255)
                                                                                                                        aaGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTccc                        <  1:2157822/149‑1 (MQ=255)
                                                                                                                          gAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGc                      >  2:3022247/1‑149 (MQ=255)
                                                                                                                          gAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATACACTCCCGc                      >  1:3031877/1‑149 (MQ=255)
                                                                                                                              tGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                  >  2:967127/1‑149 (MQ=255)
                                                                                                                              tGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                  >  2:155068/1‑149 (MQ=255)
                                                                                                                                  gcgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcag              >  2:3986353/1‑149 (MQ=255)
                                                                                                                                    gcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca            <  1:104281/149‑1 (MQ=255)
                                                                                                                                     cccACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagcag           >  1:692610/1‑149 (MQ=255)
                                                                                                                                       cACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCt         >  1:3391992/1‑149 (MQ=255)
                                                                                                                                        aCCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTa        >  1:4721715/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGAC  >  CP009273/518674‑518964

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: