Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,882,377 C→A 78.4% T96T (ACC→ACA tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,882,3770CA78.4% 528.1 / 198.5 301T96T (ACC→ACAtnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (24/41);  new base A (130/106);  total (154/147)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 303 total aligned reads displayed.

ATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATC  >  CP009273/3882231‑3882524
                                                                                                                                                  |                                                                                                                                                   
aTCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAAtt                                                                                                                                                  <  2:3798390/150‑1 (MQ=255)
 tCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTc                                                                                                                                                 <  1:4484948/150‑1 (MQ=255)
  cGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTcc                                                                                                                                                >  2:780801/1‑150 (MQ=255)
   gATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCg                                                                                                                                               >  2:2377092/1‑150 (MQ=255)
   gATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCg                                                                                                                                               <  1:4544597/150‑1 (MQ=255)
    aTTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGa                                                                                                                                              >  2:2196310/1‑150 (MQ=255)
        aCTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCa                                                                                                                                          <  1:2965394/150‑1 (MQ=255)
        aCTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCa                                                                                                                                          <  2:5158607/150‑1 (MQ=255)
           gACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCa                                                                                                                                       >  1:3416924/1‑150 (MQ=255)
           gACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCa                                                                                                                                       <  2:1777858/150‑1 (MQ=255)
            aCCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAg                                                                                                                                      <  2:2500599/150‑1 (MQ=255)
             ccGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAgg                                                                                                                                     <  2:3015807/150‑1 (MQ=255)
             ccGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAgg                                                                                                                                     >  1:3305550/1‑150 (MQ=255)
               gACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGc                                                                                                                                   <  1:4970797/150‑1 (MQ=255)
               gACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGc                                                                                                                                   >  2:4971462/1‑150 (MQ=255)
                aCAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGcc                                                                                                                                  <  1:4029541/150‑1 (MQ=255)
                aCAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGcc                                                                                                                                  <  1:5055441/150‑1 (MQ=255)
                  aGCGGCACCGGGGCGGTGACGCCGAGCCTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGGTATCAATACACAATTCCGACTCACCAGGGCCGt                                                                                                                                >  2:2258085/1‑150 (MQ=255)
                  aGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGt                                                                                                                                >  1:1279899/1‑150 (MQ=255)
                  aGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGt                                                                                                                                <  2:4052189/150‑1 (MQ=255)
                     ggCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGgc                                                                                                                             >  1:4388567/1‑150 (MQ=255)
                        aCCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCa                                                                                                                          >  2:2019057/1‑150 (MQ=255)
                        aCCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCa                                                                                                                          >  2:4590153/1‑150 (MQ=255)
                         ccGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCag                                                                                                                         >  1:3990237/1‑150 (MQ=255)
                          cGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCaga                                                                                                                        <  1:4751486/150‑1 (MQ=255)
                           ggggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                                       >  1:5154944/1‑150 (MQ=255)
                           ggggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCagag                                                                                                                       >  1:2755238/1‑150 (MQ=255)
                             ggCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCa                                                                                                                     >  1:2921710/1‑150 (MQ=255)
                              gCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCaa                                                                                                                    <  1:24223/150‑1 (MQ=255)
                               cGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCaaa                                                                                                                   >  1:1796438/1‑150 (MQ=255)
                                ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAt                                                                                                                  <  1:3604215/150‑1 (MQ=255)
                                ggTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAAt                                                                                                                  >  2:2091514/1‑150 (MQ=255)
                                 gTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATc                                                                                                                 <  1:1089413/150‑1 (MQ=255)
                                     cGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCtata                                                                                                             >  2:1910724/1‑150 (MQ=255)
                                        agagCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTc                                                                                                          >  1:684223/1‑150 (MQ=255)
                                            cATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGt                                                                                                      <  2:1883213/150‑1 (MQ=255)
                                             aTGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTa                                                                                                     <  1:1617458/150‑1 (MQ=255)
                                                cAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTg                                                                                                  >  1:1449251/1‑150 (MQ=255)
                                                 aGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCACATCTATATTCCGGTACTGa                                                                                                 >  2:1799471/1‑150 (MQ=255)
                                                            atgCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAAcgc                                                                                      <  2:2709646/150‑1 (MQ=255)
                                                                gcggcgACGCAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGc                                                                                  <  2:5040314/150‑1 (MQ=255)
                                                                 cggcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCa                                                                                 <  1:2290139/150‑1 (MQ=255)
                                                                  ggcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAg                                                                                >  2:1165005/1‑150 (MQ=255)
                                                                  ggcgACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAg                                                                                <  1:3240137/150‑1 (MQ=255)
                                                                     gacgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGac                                                                             >  2:4777476/1‑149 (MQ=255)
                                                                     gacgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaa                                                                             <  1:3202925/150‑1 (MQ=255)
                                                                      acgaAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGaaa                                                                            >  1:3777584/1‑150 (MQ=255)
                                                                            ccTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGcc                                                                      >  1:4117048/1‑150 (MQ=255)
                                                                             cTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCt                                                                     <  1:3389352/150‑1 (MQ=255)
                                                                              tACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTg                                                                    <  1:1900717/150‑1 (MQ=255)
                                                                               aCAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTgg                                                                   >  1:1980073/1‑150 (MQ=255)
                                                                                   cGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCg                                                               >  2:4133429/1‑150 (MQ=255)
                                                                                   cGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCg                                                               >  2:4195231/1‑150 (MQ=255)
                                                                                    ggCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCgc                                                              >  2:1385751/1‑150 (MQ=255)
                                                                                      cAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCgcag                                                            >  1:1008920/1‑150 (MQ=255)
                                                                                          cGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaa                                                        >  1:2642097/1‑150 (MQ=255)
                                                                                          cGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCaaa                                                        <  2:3343618/150‑1 (MQ=255)
                                                                                             aGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtg                                                     >  1:2320885/1‑150 (MQ=255)
                                                                                               ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCt                                                                     >  2:353995/1‑132 (MQ=255)
                                                                                               ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCt                                                                     <  1:353995/132‑1 (MQ=255)
                                                                                               ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGTATCGCAGCAAAAtggt                                                   >  2:2969150/1‑150 (MQ=255)
                                                                                               ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggt                                                   <  2:2844233/150‑1 (MQ=255)
                                                                                               ctactaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggt                                                   >  2:4726975/1‑150 (MQ=255)
                                                                                                  ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGc                                                <  2:2509397/150‑1 (MQ=255)
                                                                                                  ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGc                                                <  1:2918556/150‑1 (MQ=255)
                                                                                                  ctaTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGc                                                >  2:1064270/1‑150 (MQ=255)
                                                                                                   taTGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCg                                               <  1:3841392/150‑1 (MQ=255)
                                                                                                      gCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtc                                            >  1:2093274/1‑150 (MQ=255)
                                                                                                      gCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtc                                            >  1:756434/1‑150 (MQ=255)
                                                                                                      gCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtc                                            <  1:4232486/150‑1 (MQ=255)
                                                                                                       cGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                                           >  2:4462760/1‑150 (MQ=255)
                                                                                                       cGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtct                                           <  1:4640207/150‑1 (MQ=255)
                                                                                                        gTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                          <  2:4310509/150‑1 (MQ=255)
                                                                                                        gTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtctt                                          <  2:1617417/150‑2 (MQ=255)
                                                                                                        gTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                          >  2:2580885/1‑150 (MQ=255)
                                                                                                        gTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                          >  2:3659808/1‑150 (MQ=255)
                                                                                                          tAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                        <  2:4618651/150‑1 (MQ=255)
                                                                                                          tAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                        >  2:3922625/1‑150 (MQ=255)
                                                                                                             ccGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACt                                     >  1:883343/1‑150 (MQ=255)
                                                                                                              cGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTa                                    >  2:1716436/1‑150 (MQ=255)
                                                                                                              cGAGTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTa                                    >  1:4862694/1‑150 (MQ=255)
                                                                                                                 gTCAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAttt                                 >  1:2195153/1‑150 (MQ=255)
                                                                                                                   cAGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                               >  2:1519153/1‑150 (MQ=255)
                                                                                                                    aGTGAAAAATATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                              <  1:4175239/150‑1 (MQ=255)
                                                                                                                      tGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTg                            >  1:365386/1‑150 (MQ=255)
                                                                                                                          aaaTATCTTTGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                        >  1:1115334/1‑150 (MQ=255)
                                                                                                                                cTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                  <  2:3181413/150‑1 (MQ=255)
                                                                                                                                   tGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg               <  2:330854/150‑1 (MQ=255)
                                                                                                                                   tGGTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg               >  1:1725676/1‑150 (MQ=255)
                                                                                                                                     gTTATCAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGcc             <  1:4787049/150‑1 (MQ=255)
                                                                                                                                          cAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc        >  2:4862110/1‑150 (MQ=255)
                                                                                                                                          cAATACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc        >  2:889800/1‑150 (MQ=255)
                                                                                                                                           aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc       <  2:3692303/150‑1 (MQ=255)
                                                                                                                                           aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc       >  1:1240728/1‑150 (MQ=255)
                                                                                                                                           aaTACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc       >  2:2646962/1‑150 (MQ=255)
                                                                                                                                           aaTACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc       >  2:882821/1‑150 (MQ=255)
                                                                                                                                            aTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa      <  1:1334784/150‑1 (MQ=255)
                                                                                                                                            aTACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa      <  1:928555/150‑1 (MQ=255)
                                                                                                                                             tACACAATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAg     >  2:2393619/1‑150 (MQ=255)
                                                                                                                                                acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc  <  2:3232749/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGCGTTAGCCGAGTCAGTGAAAAATATCTTTGGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATC  >  CP009273/3882231‑3882524

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: