Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,350,891 G→T 84.2% S159R (AGC→AGA) ‡ sapA ← antimicrobial peptide transport ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,350,8910GT84.2% 621.4 / 132.2 278S159R (AGC→AGA) ‡sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base T (123/111);  total (144/134)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 281 total aligned reads displayed.

AGTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACC  >  CP009273/1350745‑1351036
                                                                                                                                                  |                                                                                                                                                 
aGTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGActct                                                                                                                                                 <  2:339459/149‑1 (MQ=255)
 gTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCtt                                                                                                                                                >  1:4324513/1‑149 (MQ=255)
 gTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCtt                                                                                                                                                <  1:903129/149‑1 (MQ=255)
  ttGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCttt                                                                                                                                               <  2:291889/149‑1 (MQ=255)
  ttGCTCCTGGCGATCTTCTTTCTCGACCTTCCGGGCATATTCGGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAACAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCttt                                                                                                                                               <  2:4082836/149‑1 (MQ=255)
      tCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCAc                                                                                                                                           <  2:616189/149‑1 (MQ=255)
        cTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGt                                                                                                                                         <  2:2851690/149‑1 (MQ=255)
        cTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGt                                                                                                                                         >  1:3048211/1‑149 (MQ=255)
           gCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTAt                                                                                                                                      >  1:109635/1‑149 (MQ=255)
            cGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATc                                                                                                                                     <  1:2107890/149‑1 (MQ=255)
            cGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATc                                                                                                                                     <  2:4709845/149‑1 (MQ=255)
             gATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCg                                                                                                                                    <  1:1784536/149‑1 (MQ=255)
             gATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCg                                                                                                                                    >  2:4251618/1‑149 (MQ=255)
               tcttctTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGc                                                                                                                                  >  2:3045373/1‑149 (MQ=255)
               tcttctTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGc                                                                                                                                  >  1:1102674/1‑149 (MQ=255)
                 ttctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaa                                                                                                                                >  2:2123390/1‑149 (MQ=255)
                 ttctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaa                                                                                                                                >  1:3220931/1‑149 (MQ=255)
                 ttctttctCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCaa                                                                                                                                >  1:785514/1‑149 (MQ=255)
                        tctAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAg                                                                                                                         >  2:671318/1‑149 (MQ=255)
                        tctAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAg                                                                                                                         <  1:3487599/149‑1 (MQ=255)
                          tAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGc                                                                                                                       >  1:190109/1‑149 (MQ=255)
                           aaCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCt                                                                                                                      <  2:2515761/149‑1 (MQ=255)
                           aaCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCt                                                                                                                      <  1:531825/149‑1 (MQ=255)
                           aaCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCt                                                                                                                      >  1:260418/1‑149 (MQ=255)
                           aaCTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTACACGTTATCGGCAAATTGCAGGCt                                                                                                                      >  1:315264/1‑149 (MQ=255)
                                ccGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGa                                                                                                                 >  1:3389432/1‑149 (MQ=255)
                                 cGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGaa                                                                                                                <  1:21200/149‑1 (MQ=255)
                                    gCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTa                                                                                                             <  2:1983591/149‑1 (MQ=255)
                                    gCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTa                                                                                                             <  1:1201436/149‑1 (MQ=255)
                                    gCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTa                                                                                                             <  2:1289876/149‑1 (MQ=255)
                                       tatTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAggg                                                                                                          >  2:4575026/1‑149 (MQ=255)
                                       tatTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAggg                                                                                                          >  1:2813832/1‑149 (MQ=255)
                                        atTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAgggg                                                                                                         >  1:183990/1‑149 (MQ=255)
                                        atTCTGCCGACATGAACGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGGTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAgggg                                                                                                         >  2:3320231/1‑149 (MQ=255)
                                            tGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGt                                                                                                     <  2:4686666/149‑1 (MQ=255)
                                               cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAGGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgc                                                                                                  <  2:3907847/149‑1 (MQ=255)
                                               cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgc                                                                                                  <  2:3837945/149‑1 (MQ=255)
                                               cGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgc                                                                                                  <  2:1091687/149‑1 (MQ=255)
                                                 aCATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgctg                                                                                                >  1:3678592/1‑149 (MQ=255)
                                                  cATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGACGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTtgctgc                                                                                               >  2:2927978/1‑149 (MQ=255)
                                                     gACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGt                                                                                            >  2:2248506/1‑149 (MQ=255)
                                                      aCCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGtt                                                                                           >  2:1788143/1‑149 (MQ=255)
                                                      aCCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGACTGTCGAAGTAGGGGAAGTTGCTGCCGtt                                                                                           <  2:1649691/149‑1 (MQ=255)
                                                        cGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGa                                                                                         <  1:3701263/149‑1 (MQ=255)
                                                        cGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGa                                                                                         >  1:1008395/1‑149 (MQ=255)
                                                          aaGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACg                                                                                       >  1:2996947/1‑149 (MQ=255)
                                                             cATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTa                                                                                    >  1:998016/1‑149 (MQ=255)
                                                              aTAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCAGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTAt                                                                                   >  2:1389583/1‑149 (MQ=255)
                                                                 aTGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGcc                                                                                >  2:1921716/1‑149 (MQ=255)
                                                                 aTGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGcc                                                                                >  2:1037471/1‑149 (MQ=255)
                                                                 aTGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGcc                                                                                >  2:3810757/1‑149 (MQ=255)
                                                                  tGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCa                                                                               >  1:1266523/1‑149 (MQ=255)
                                                                     gTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACgg                                                                            <  2:957264/149‑1 (MQ=255)
                                                                     gTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACgg                                                                            >  2:211690/1‑149 (MQ=255)
                                                                      ttGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACggg                                                                           <  2:3597429/149‑1 (MQ=255)
                                                                       tGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTGATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGgt                                                                          <  1:2187751/149‑1 (MQ=255)
                                                                           aGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGgttgtt                                                                      >  1:267365/1‑149 (MQ=255)
                                                                             gTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGACGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGc                                                                    >  2:4157274/1‑149 (MQ=255)
                                                                                ccACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGt                                                                 >  1:4686800/1‑149 (MQ=255)
                                                                                         agaagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTc                                                        <  1:2735914/149‑1 (MQ=255)
                                                                                         agaagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTc                                                        <  2:4487596/149‑1 (MQ=255)
                                                                                           aagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGc                                                      >  1:3450824/1‑149 (MQ=255)
                                                                                           aagCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGc                                                      <  1:4729286/149‑1 (MQ=255)
                                                                                             gCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTg                                                    >  2:3986262/1‑149 (MQ=255)
                                                                                             gCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTg                                                    >  1:3241582/1‑149 (MQ=255)
                                                                                                  cGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAgg                                               <  1:1611509/149‑1 (MQ=255)
                                                                                                  cGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAgg                                               <  1:3189765/149‑1 (MQ=255)
                                                                                                    gCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTa                                             <  1:2299029/149‑1 (MQ=255)
                                                                                                    gCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTa                                             <  2:3736643/149‑1 (MQ=255)
                                                                                                       gggCGAGTCGGAACTCAACGGTATGATTATCCAGGTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGGAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAac                                          >  2:3157653/1‑149 (MQ=255)
                                                                                                         gCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAacac                                        <  2:3221550/149‑1 (MQ=255)
                                                                                                           gAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAAcacca                                      <  2:260409/149‑1 (MQ=255)
                                                                                                             gTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACcaca                                    >  1:714259/1‑149 (MQ=255)
                                                                                                                ggAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcg                                 >  2:79105/1‑149 (MQ=255)
                                                                                                                 gAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgt                                >  1:3624839/1‑149 (MQ=255)
                                                                                                                 gAACTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgt                                <  2:2641934/149‑1 (MQ=255)
                                                                                                                 gAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgt                                <  1:4211731/149‑1 (MQ=255)
                                                                                                                 gAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgt                                <  1:1822390/149‑1 (MQ=255)
                                                                                                                  aaCTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgtc                               <  1:2660498/149‑1 (MQ=255)
                                                                                                                  aaCTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGTGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgtc                               <  2:1554312/149‑1 (MQ=255)
                                                                                                                   aCTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACAtcgtcg                              >  2:2961661/1‑149 (MQ=255)
                                                                                                                   aCTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAATTCGCTGGAAGGTAAACACCACATCGTCgg                             >  1:606820/1‑149 (MQ=255)
                                                                                                                    cTCAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCgg                             >  2:2407061/1‑149 (MQ=255)
                                                                                                                      cAACGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGca                           >  2:2776873/1‑149 (MQ=255)
                                                                                                                       aaCGGTATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcat                          >  1:324870/1‑149 (MQ=255)
                                                                                                                       aaCGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATAGGCAAATTGCAGGCTGTAGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATTCCAAGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAAAACCACATAGTCGGcat                          >  1:2445207/1‑149 (MQ=255)
                                                                                                                            tATGATTATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcatt                     <  1:2019010/149‑1 (MQ=255)
                                                                                                                            tATGATTATCCAGTTTGCGGACTATTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcatt                     <  2:694385/149‑1 (MQ=255)
                                                                                                                            tATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGcattcatt                     <  2:1584489/149‑1 (MQ=255)
                                                                                                                                 ttATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTAct                <  1:4023622/149‑2 (MQ=255)
                                                                                                                                 ttATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTAct                >  2:3927677/1‑148 (MQ=255)
                                                                                                                                  tATCCAGTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTActa               <  1:3952852/149‑3 (MQ=255)
                                                                                                                                        gTTTGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGa         <  2:869794/149‑1 (MQ=255)
                                                                                                                                          ttGCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGt       <  1:3295082/149‑1 (MQ=255)
                                                                                                                                            gCGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTaa     >  1:3796983/1‑149 (MQ=255)
                                                                                                                                             cGGACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTaaa    <  1:1563635/149‑1 (MQ=255)
                                                                                                                                              ggACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAc   <  1:2290587/149‑1 (MQ=255)
                                                                                                                                              ggACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAc   <  1:4575026/149‑1 (MQ=255)
                                                                                                                                               gACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc  >  1:4299232/1‑149 (MQ=255)
                                                                                                                                               gACTCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACTAGTGGGAGTAAAcc  >  1:2949979/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGTTGCTCCTGGCGATCTTCTTTCTCTAACTTCCGGGCATATTCTGCCGACATGACCGAAGCATAATGGGTTGCGAGGTGCCACAAAAAAGAAGCATCCGGCTGGGCGAGTCGGAACTCAACGGTATGATTATCCAGTTTGCGGACGCTTTTCACGTTATCGGCAAATTGCAGGCTGTCGAAGTAGGGGAAGTTGCTGCCGTTGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACC  >  CP009273/1350745‑1351036

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: