Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,791,973 Δ1 bp intergenic (+67/‑68) 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,9730A.100.0% 2731.3 / NA 549intergenic (+67/‑68)nadE/choEammonia‑dependent NAD(+) synthetase/cholinesterase ChoE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (292/257);  total (292/257)

Only 100 of 549 total aligned reads displayed.

GCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCT  >  NC_008463/5791828‑5792117
                                                                                                                                                 |                                                                                                                                                
 cgcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAAcc                                                                                                                                            <  1:1990460/150‑1 (MQ=255)
   cgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAg                                                                                                                                          >  2:2751939/1‑150 (MQ=255)
   cgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAg                                                                                                                                          >  2:2463433/1‑150 (MQ=255)
    gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGCCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAgg                                                                                                                                         >  1:894491/1‑150 (MQ=255)
     aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                                                        <  2:6932442/150‑1 (MQ=255)
         aGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTc                                                                                                                                    >  2:223152/1‑150 (MQ=255)
           ccGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCaa                                                                                                                                  <  1:2604645/150‑1 (MQ=255)
            cgcgcgGGTGATCGTCGACACCTACCCCAAGACCCAGCCCAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAg                                                                                                                                 <  2:7607823/150‑1 (MQ=255)
              cgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTc                                                                                                                               <  2:2615436/150‑1 (MQ=255)
               gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                                                              >  2:10228557/1‑150 (MQ=255)
                  ggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                                                           >  1:1271864/1‑150 (MQ=255)
                   gTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgc                                                                                                                          <  2:9360104/150‑1 (MQ=255)
                    tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                                                         <  1:10336197/150‑1 (MQ=255)
                    tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                                                         >  1:8640039/1‑150 (MQ=255)
                     gATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcgg                                                                                                                        <  1:8021653/150‑1 (MQ=255)
                      aTCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                                                       >  1:3373594/1‑150 (MQ=255)
                        cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                                                     <  1:276755/150‑1 (MQ=255)
                        cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                                                     <  1:5656617/150‑1 (MQ=255)
                         gtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                                                    <  2:87786/150‑1 (MQ=255)
                          tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                                                   >  1:5771727/1‑150 (MQ=255)
                           cgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                                                  >  2:1670271/1‑150 (MQ=255)
                            gACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcg                                                                                                                 >  2:4270142/1‑150 (MQ=255)
                             acacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                                                <  1:9484334/150‑1 (MQ=255)
                             acacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                                                <  2:1194682/150‑1 (MQ=255)
                               acCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                              <  2:3189074/150‑1 (MQ=255)
                               acCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                              >  2:5160370/1‑150 (MQ=255)
                               acCTACCACAAGACCCAGCACAAGCGCGAACAGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                              <  1:5476698/150‑1 (MQ=255)
                                  tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAGTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                                                           >  1:514520/1‑150 (MQ=255)
                                  tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                                                           >  2:5777132/1‑150 (MQ=255)
                                      acaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                                                       <  2:1809182/150‑1 (MQ=255)
                                         aGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCGTCAAGTCCGCGCGGCGGCGCGCGAGACGGAAGCt                                                                                                    >  2:4017651/1‑150 (MQ=255)
                                         aGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGTATGTGCCTGTCGCAACCAGGTCTTCATGTCCGCGCGGCGGCGCGCGAGAAggtcgct                                                                                                    >  1:181532/1‑145 (MQ=255)
                                          gACCCAGCTCATGCGCGACCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAGGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTg                                                                                                   <  2:6309354/150‑1 (MQ=255)
                                           aCCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGa                                                                                                  >  2:4462096/1‑150 (MQ=255)
                                            cccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGaa                                                                                                 <  1:7583284/150‑1 (MQ=255)
                                                gCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                                                             <  1:8813507/150‑1 (MQ=255)
                                                  acaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                                                           <  1:6126812/150‑1 (MQ=255)
                                                  acaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                                                           <  1:5556594/150‑1 (MQ=255)
                                                  acaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                                                           <  1:134840/150‑1 (MQ=255)
                                                    aaGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCg                                                                                         >  1:4082418/1‑150 (MQ=255)
                                                        gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgg                                                                                     <  2:5657826/150‑1 (MQ=255)
                                                        gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgg                                                                                     >  2:3051756/1‑150 (MQ=255)
                                                         cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGACGCTGAATCCCTGCGAAGgc                                                                                    >  2:5240739/1‑150 (MQ=255)
                                                           aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCTGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCGAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                  >  1:4991700/1‑150 (MQ=255)
                                                           aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCTTCGTCGACGGGCCTTTGTCGCGGGCGGGTATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                  <  2:7331228/150‑1 (MQ=255)
                                                           aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                  >  2:1998838/1‑150 (MQ=255)
                                                            aCTGCCGAAGGCGCCCAGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                                                 <  1:1651708/150‑1 (MQ=255)
                                                               gCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcca                                                                              >  1:6837564/1‑150 (MQ=255)
                                                                ccGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCAt                                                                             >  2:6144783/1‑150 (MQ=255)
                                                                 cGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                                                            >  1:9794235/1‑150 (MQ=255)
                                                                 cGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTGTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                                                            >  2:6593162/1‑150 (MQ=255)
                                                                   aaGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                                                          <  1:4487850/150‑1 (MQ=255)
                                                                    aGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcgc                                                                         <  2:8282130/150‑1 (MQ=255)
                                                                    aGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcgc                                                                         <  2:3297431/150‑1 (MQ=255)
                                                                        gcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                                                     <  2:8623052/150‑1 (MQ=255)
                                                                        gcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                                                     <  1:10000877/150‑1 (MQ=255)
                                                                         cccTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCg                                                                    <  1:10131298/150‑1 (MQ=255)
                                                                         cccTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCg                                                                    >  1:1188107/1‑150 (MQ=255)
                                                                           cTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtc                                                                  >  1:2329746/1‑150 (MQ=255)
                                                                             gAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtctc                                                                <  1:2277938/150‑1 (MQ=255)
                                                                             gAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtctc                                                                <  1:148520/150‑1 (MQ=255)
                                                                                    ccGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGccc                                                         <  2:4105505/150‑1 (MQ=255)
                                                                                     cGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCt                                                        >  1:7094411/1‑150 (MQ=255)
                                                                                       gCCCGTGCGGCAGGCCCGGCGTCGACGGGCCTTTGTCGCGGGCGGGAAGGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                                                      >  2:3638189/1‑150 (MQ=255)
                                                                                         ccGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctc                                                    >  1:5257385/1‑150 (MQ=255)
                                                                                             gCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGATTCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCg                                                <  2:5823065/150‑1 (MQ=255)
                                                                                               ggCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGACCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgt                                              <  2:4973525/150‑1 (MQ=255)
                                                                                                  aGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCGAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTc                                           >  1:6001396/1‑150 (MQ=255)
                                                                                                    gCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCct                                         >  1:7187873/1‑150 (MQ=255)
                                                                                                    gCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCct                                         >  2:3902991/1‑150 (MQ=255)
                                                                                                     cccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctg                                        >  1:6617490/1‑150 (MQ=255)
                                                                                                        gtcgtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgctg                                     >  2:8501880/1‑150 (MQ=255)
                                                                                                         tcgtcgACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTgg                                    >  1:6374075/1‑150 (MQ=255)
                                                                                                         tcgtcgACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTgg                                    >  1:927007/1‑150 (MQ=255)
                                                                                                          cgtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggc                                   <  2:1620014/150‑1 (MQ=255)
                                                                                                           gtcGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggcc                                  >  2:3453578/1‑150 (MQ=255)
                                                                                                           gtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggcc                                  >  1:3250372/1‑150 (MQ=255)
                                                                                                            tcgaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccg                                 >  2:4204961/4‑150 (MQ=255)
                                                                                                            tcgaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccg                                 >  1:3024545/4‑150 (MQ=255)
                                                                                                              gaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccggc                               >  1:4746411/2‑150 (MQ=255)
                                                                                                              gaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccggc                               <  2:10649175/149‑1 (MQ=255)
                                                                                                                cGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                             <  2:8025459/150‑1 (MQ=255)
                                                                                                                cGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                             >  1:8365393/1‑150 (MQ=255)
                                                                                                                cGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                             >  2:6427964/1‑150 (MQ=255)
                                                                                                                 gggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAg                            >  2:9614413/1‑150 (MQ=255)
                                                                                                                  ggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAgg                           >  2:4626146/1‑150 (MQ=255)
                                                                                                                  ggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAgg                           <  1:8181754/150‑1 (MQ=255)
                                                                                                                        tGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCg                     <  1:10496821/150‑1 (MQ=255)
                                                                                                                        tGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCg                     >  2:9421355/1‑150 (MQ=255)
                                                                                                                           cgcgGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGac                  >  1:2460350/1‑150 (MQ=255)
                                                                                                                            gcgggcggGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGaca                 >  1:4281136/1‑150 (MQ=255)
                                                                                                                                    gAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTccc         >  1:4018980/1‑150 (MQ=255)
                                                                                                                                    gAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTccc         >  1:7718253/1‑150 (MQ=255)
                                                                                                                                    gAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTccc         >  2:8875044/1‑150 (MQ=255)
                                                                                                                                    gAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTccc         >  2:9083593/1‑150 (MQ=255)
                                                                                                                                      aTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCg       >  2:6737264/1‑150 (MQ=255)
                                                                                                                                       tgtgCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCAGGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgc      <  1:2091380/150‑1 (MQ=255)
                                                                                                                                       ggtgCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTCGACGGAAGCTGACTCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCGGGCCGGCCAGGCCTGCGGACACACGTCCCCgc      <  1:10633313/149‑1 (MQ=255)
                                                                                                                                           ccTGTCGCACCCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgctgct  <  2:713921/150‑1 (MQ=255)
                                                                                                                                           ccTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgctgct  >  2:7750909/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCT  >  NC_008463/5791828‑5792117

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: