Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,791,938 A→G intergenic (+32/‑103) nadE → / → choE ammonia‑dependent NAD(+) synthetase/cholinesterase ChoE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,791,9380AG100.0% 2268.7 / NA 580intergenic (+32/‑103)nadE/choEammonia‑dependent NAD(+) synthetase/cholinesterase ChoE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (298/282);  total (298/282)

Only 100 of 582 total aligned reads displayed.

TATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCG  >  NC_008463/5791790‑5792086
                                                                                                                                                    |                                                                                                                                                    
tATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGAc                                                                                                                                                    >  1:7763654/1‑151 (MQ=255)
 aTGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACg                                                                                                                                                   >  2:1982736/1‑151 (MQ=255)
  tGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACgg                                                                                                                                                  >  1:6220078/1‑151 (MQ=255)
  tGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACgg                                                                                                                                                  >  2:10922018/1‑151 (MQ=255)
        gaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCttt                                                                                                                                            <  2:10310748/151‑1 (MQ=255)
        gaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCttt                                                                                                                                            <  1:6749317/151‑1 (MQ=255)
         aTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTg                                                                                                                                           >  2:9924353/1‑151 (MQ=255)
         aTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTg                                                                                                                                           <  1:10677676/151‑1 (MQ=255)
          tCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGt                                                                                                                                          <  2:3278919/151‑1 (MQ=255)
            gACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcg                                                                                                                                        >  2:4398285/1‑151 (MQ=255)
              cGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGCAGCCGCGCGGGTGAGCGTCGACCCCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGGGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcgcg                                                                                                                                      >  1:719148/1‑151 (MQ=255)
                       gCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAt                                                                                                                             <  1:4502462/151‑1 (MQ=255)
                       gCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAt                                                                                                                             <  1:9376280/151‑1 (MQ=255)
                        cACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCCCAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACtg                                                                                                                            >  2:4535606/1‑151 (MQ=255)
                        cACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCCCAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtg                                                                                                                            <  1:436265/151‑1 (MQ=255)
                        cACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACtg                                                                                                                            >  2:100127/1‑151 (MQ=255)
                         aCGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgt                                                                                                                           >  2:4980499/1‑151 (MQ=255)
                          cGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgtg                                                                                                                          >  1:7597210/1‑151 (MQ=255)
                           ggCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGc                                                                                                                         >  1:9548311/1‑151 (MQ=255)
                            gCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGcc                                                                                                                        <  1:9778117/151‑1 (MQ=255)
                             ccAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCt                                                                                                                       <  2:9668929/151‑1 (MQ=255)
                                   gCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcgca                                                                                                                 >  1:4732231/1‑148 (MQ=255)
                                   gCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcgca                                                                                                                 >  2:9186436/1‑148 (MQ=255)
                                   gCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcgca                                                                                                                 <  1:299722/151‑4 (MQ=255)
                                    cTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcgcaa                                                                                                                >  2:10509067/1‑147 (MQ=255)
                                    cTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcgcaa                                                                                                                <  1:2292428/151‑5 (MQ=255)
                                     tGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAAc                                                                                                              <  2:4085702/151‑1 (MQ=255)
                                        gcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAg                                                                                                           >  2:7100249/1‑151 (MQ=255)
                                          gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGt                                                                                                         >  2:1910183/1‑151 (MQ=255)
                                              aaGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTc                                                                                                     >  2:5237346/1‑151 (MQ=255)
                                              aaGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTc                                                                                                     >  1:3450248/1‑151 (MQ=255)
                                              aaGCCGCGCGGGTGATCGTCGACACCGACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGCGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTc                                                                                                     <  2:8015599/151‑1 (MQ=255)
                                               aGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCa                                                                                                    <  1:9221129/151‑1 (MQ=255)
                                                gCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACCAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTTTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCaa                                                                                                   >  2:5027855/1‑151 (MQ=255)
                                                 ccGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAg                                                                                                  <  1:9632650/151‑1 (MQ=255)
                                                  cgcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCTGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGt                                                                                                 >  2:4302966/1‑151 (MQ=255)
                                                   gcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCTAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTc                                                                                                >  1:573640/1‑151 (MQ=255)
                                                   gcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTc                                                                                                >  1:8899339/1‑151 (MQ=255)
                                                    cgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTcc                                                                                               >  1:11028942/1‑151 (MQ=255)
                                                     gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCcg                                                                                              <  2:3539395/151‑1 (MQ=255)
                                                      cgGGTGATCGTCGACACCTACCACCAGACCCCGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCcgc                                                                                             <  1:10835655/151‑1 (MQ=255)
                                                      cgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCcgc                                                                                             <  2:940298/151‑1 (MQ=255)
                                                          tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcgg                                                                                         <  1:2909790/151‑1 (MQ=255)
                                                           gATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggc                                                                                        <  1:10645456/151‑1 (MQ=255)
                                                           gATCGTCGACACCTAACACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCAGCgcggc                                                                                        <  2:1767839/151‑1 (MQ=255)
                                                           gATCGTCGAAACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGAgcggc                                                                                        >  2:6812959/1‑151 (MQ=255)
                                                            ctcgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTTCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcg                                                                                       <  2:7303307/150‑1 (MQ=255)
                                                             tcgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                      >  2:6641283/1‑151 (MQ=255)
                                                              cgtcgACACCTACCACAAGACCCAGCCCAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                     <  1:789344/151‑1 (MQ=255)
                                                              cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                     >  1:9698230/1‑151 (MQ=255)
                                                              cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                     >  1:400261/1‑151 (MQ=255)
                                                               gtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                    >  1:4403519/1‑151 (MQ=255)
                                                                 cgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcg                                                                                  <  1:603928/151‑1 (MQ=255)
                                                                  gACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGGCGACGGGCCGTTGTCGCGGGCGGGAATGAGCCAGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                 >  2:6364893/1‑151 (MQ=255)
                                                                   acacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCt                                                                                >  1:387496/1‑151 (MQ=255)
                                                                    cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGTATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                               <  1:2051696/151‑1 (MQ=255)
                                                                    cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgcga                                                                               >  2:2531563/1‑149 (MQ=255)
                                                                    cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                               >  1:1303849/1‑151 (MQ=255)
                                                                     acCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTag                                                                              <  2:4721780/151‑2 (MQ=255)
                                                                     acCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGCATGTGCCTGTC‑GAAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTag                                                                              >  1:846957/1‑150 (MQ=255)
                                                                      ccTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTaga                                                                             <  1:5620640/151‑3 (MQ=255)
                                                                       cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCtcgaa                                                                            >  2:7619211/1‑147 (MQ=255)
                                                                        tACCACAAGTCCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCTGCCCGTGCGGCTGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAg                                                                           <  2:11214053/151‑1 (MQ=255)
                                                                           cacaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                        <  2:151311/151‑1 (MQ=255)
                                                                            acaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAgaaggaag                                                                       <  2:5876426/151‑1 (MQ=255)
                                                                             caAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                      <  1:8571186/151‑1 (MQ=255)
                                                                              aaGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCt                                                                     <  1:6325763/151‑1 (MQ=255)
                                                                                gACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGa                                                                   <  1:6617656/151‑1 (MQ=255)
                                                                                gACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGa                                                                   <  1:1027818/151‑1 (MQ=255)
                                                                                   ccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATc                                                                >  1:5000011/1‑151 (MQ=255)
                                                                                    cAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATcc                                                               <  2:7177888/151‑1 (MQ=255)
                                                                                        acaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGc                                                           >  2:7644961/1‑151 (MQ=255)
                                                                                         caAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCg                                                          >  1:4181960/1‑151 (MQ=255)
                                                                                          aaGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGa                                                         >  1:1688812/1‑151 (MQ=255)
                                                                                                gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                   <  2:10762222/151‑1 (MQ=255)
                                                                                                gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                   >  1:2528430/1‑151 (MQ=255)
                                                                                                 aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                  <  1:7404489/151‑1 (MQ=255)
                                                                                                   cTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcc                                                <  1:830175/151‑1 (MQ=255)
                                                                                                   cTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcc                                                <  2:2136947/151‑1 (MQ=255)
                                                                                                    tGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcca                                               >  2:3120931/1‑151 (MQ=255)
                                                                                                      ccGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                             >  2:8190329/1‑151 (MQ=255)
                                                                                                      ccGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                             <  1:8128713/151‑1 (MQ=255)
                                                                                                      ccGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                             <  1:308745/151‑1 (MQ=255)
                                                                                                             cgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                      <  1:10066213/151‑1 (MQ=255)
                                                                                                             cgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                      <  1:8729273/151‑1 (MQ=255)
                                                                                                               cccTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGt                                    <  1:5392633/151‑1 (MQ=255)
                                                                                                                        tCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGcc                           <  2:1293577/151‑1 (MQ=255)
                                                                                                                        tCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGcc                           <  2:7763654/151‑1 (MQ=255)
                                                                                                                         cccGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGccc                          >  1:373339/1‑151 (MQ=255)
                                                                                                                             gCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGct                      >  2:8489355/1‑151 (MQ=255)
                                                                                                                               ccGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcc                    <  2:8861189/151‑1 (MQ=255)
                                                                                                                                 gTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcctc                  <  2:2321791/151‑1 (MQ=255)
                                                                                                                                 gTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcctc                  >  2:2970174/1‑151 (MQ=255)
                                                                                                                                     ggCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtg              >  2:5443627/1‑151 (MQ=255)
                                                                                                                                     ggCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtg              <  1:2723218/151‑1 (MQ=255)
                                                                                                                                      gCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtgt             >  1:10401381/1‑151 (MQ=255)
                                                                                                                                      gCACGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTCGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtgt             <  2:6148460/151‑1 (MQ=255)
                                                                                                                                           cccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgc        >  2:8913118/1‑151 (MQ=255)
                                                                                                                                              gtcgtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTgg     >  1:1471447/1‑151 (MQ=255)
                                                                                                                                               tcgtcgACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTC‑GCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggc    >  2:3834409/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCG  >  NC_008463/5791790‑5792086

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: ‑