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,9380AG99.8% 2215.5 / NA 566intergenic (+32/‑103)nadE/choEammonia‑dependent NAD(+) synthetase/cholinesterase ChoE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (274/291);  minor base T (1/0);  total (275/291)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 570 total aligned reads displayed.

ATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCC  >  NC_008463/5791791‑5792085
                                                                                                                                                   |                                                                                                                                                   
 tGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACg                                                                                                                                                  <  1:3802291/150‑1 (MQ=255)
   ccGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACggg                                                                                                                                                >  1:1386451/1‑150 (MQ=255)
         tCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTg                                                                                                                                          >  2:506382/1‑150 (MQ=255)
          cGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGt                                                                                                                                         <  1:10389765/150‑1 (MQ=255)
            aCGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcg                                                                                                                                       <  1:8410507/150‑1 (MQ=255)
            aCGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcg                                                                                                                                       >  1:9853773/1‑150 (MQ=255)
            aCGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcg                                                                                                                                       >  2:8456863/1‑150 (MQ=255)
               ccTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCgcgg                                                                                                                                    >  2:9410190/1‑150 (MQ=255)
                    cTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGcgggcggg                                                                                                                               >  1:10021910/1‑150 (MQ=255)
                     agCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGa                                                                                                                              <  2:6995431/149‑1 (MQ=255)
                      gCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGac                                                                                                                             >  1:9996599/1‑149 (MQ=255)
                       cACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAt                                                                                                                            >  2:2599316/1‑150 (MQ=255)
                        aCGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtg                                                                                                                           <  1:512869/150‑1 (MQ=255)
                         cGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACACGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgt                                                                                                                          <  2:5740778/150‑1 (MQ=255)
                         cGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgt                                                                                                                          >  1:1813567/1‑150 (MQ=255)
                           gCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGc                                                                                                                        >  1:4180705/1‑150 (MQ=255)
                            ccAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGcc                                                                                                                       >  1:6280868/1‑150 (MQ=255)
                              aGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTg                                                                                                                     <  1:1610680/150‑1 (MQ=255)
                               gCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGt                                                                                                                    <  2:4453292/150‑1 (MQ=255)
                               gCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGt                                                                                                                    <  2:3414952/150‑1 (MQ=255)
                                 cGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTcg                                                                                                                  >  1:10556985/1‑149 (MQ=255)
                                     gcgcGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAAc                                                                                                             <  1:1917749/150‑1 (MQ=255)
                                        cgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAg                                                                                                          <  2:2849237/150‑1 (MQ=255)
                                         gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAgg                                                                                                         >  2:3071902/1‑150 (MQ=255)
                                          aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACATGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                        <  2:633563/150‑1 (MQ=255)
                                          aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                        <  1:8163212/150‑1 (MQ=255)
                                          aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                        >  2:11038386/1‑150 (MQ=255)
                                           ggAAGCCGCGCGGGTGATCGTCGACACCTACCCCAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTc                                                                                                       <  2:3845628/150‑1 (MQ=255)
                                            gAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCt                                                                                                      <  1:8636743/150‑1 (MQ=255)
                                             aaGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCtt                                                                                                     >  2:7941849/1‑150 (MQ=255)
                                                   cgcgGGTGATCGTCGACACCTACCCCAAAACCCAGCACAACCGCCAACATCCGTAGGCGCCCTGTTCCTCCCCGCCCGTTCGGCAGGCCCGTCGTCGACGGGCCGTTTTCGCGGGCGGGAAATTGCCAGTCGCAAACAAGTCATAAAGTc                                                                                               >  2:1810960/1‑150 (MQ=255)
                                                    gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                              <  2:156255/150‑1 (MQ=255)
                                                    gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                              >  1:978413/1‑150 (MQ=255)
                                                      gggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgc                                                                                            >  2:6158436/1‑150 (MQ=255)
                                                       ggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                           <  1:643692/150‑1 (MQ=255)
                                                         tGATCGTCGACACCTACCACAAGACCCAGCACACGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                         >  2:9224806/1‑150 (MQ=255)
                                                          gATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcgg                                                                                        <  1:7510804/150‑1 (MQ=255)
                                                           aTCGTCGACCCCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGCGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                       <  2:2134785/150‑1 (MQ=255)
                                                           aTCGTCGACACCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                       <  2:1649575/150‑1 (MQ=255)
                                                           aTCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                       >  1:7738878/1‑150 (MQ=255)
                                                           aTCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                       <  1:3540594/150‑1 (MQ=255)
                                                             cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTNGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                     <  1:5876184/150‑1 (MQ=255)
                                                              gtcgACACCTACCACACGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                    <  1:5613197/150‑1 (MQ=255)
                                                               tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                   >  1:1018406/1‑150 (MQ=255)
                                                                cgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                  >  1:3320017/1‑150 (MQ=255)
                                                                    acCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                              >  1:9667620/1‑150 (MQ=255)
                                                                     ccTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTag                                                                             >  1:6656395/1‑149 (MQ=255)
                                                                      cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTaga                                                                            <  1:4022188/150‑3 (MQ=255)
                                                                       tACCACACGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGCCTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgcgagaa                                                                           >  2:528974/1‑145 (MQ=255)
                                                                       tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                           >  1:457764/1‑150 (MQ=255)
                                                                       tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                           <  1:2762327/150‑1 (MQ=255)
                                                                       tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                           >  1:6176759/1‑150 (MQ=255)
                                                                         ccACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAgg                                                                         >  2:9609767/1‑150 (MQ=255)
                                                                         ccACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAgg                                                                         <  1:2068408/150‑1 (MQ=255)
                                                                           acaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                       >  1:4367980/1‑150 (MQ=255)
                                                                               gACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTg                                                                   <  1:4628157/150‑1 (MQ=255)
                                                                                 cccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGaa                                                                 >  2:1963775/1‑150 (MQ=255)
                                                                                 cccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGaa                                                                 <  2:7698517/150‑1 (MQ=255)
                                                                                  ccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAAt                                                                >  1:5983171/1‑150 (MQ=255)
                                                                                   cAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATc                                                               >  2:8570415/1‑150 (MQ=255)
                                                                                   cAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATc                                                               >  2:5528566/1‑150 (MQ=255)
                                                                                    aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATcc                                                              <  2:6553649/150‑1 (MQ=255)
                                                                                    aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATcc                                                              >  1:9340319/1‑150 (MQ=255)
                                                                                     gctcaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                             <  1:4932671/147‑1 (MQ=255)
                                                                                      cacaCGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCt                                                            <  2:4017553/150‑1 (MQ=255)
                                                                                        caAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGc                                                          <  1:3167917/150‑1 (MQ=255)
                                                                                          aGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGa                                                        >  1:6507255/1‑150 (MQ=255)
                                                                                             gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCACGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgt                                                     <  2:4845383/150‑2 (MQ=255)
                                                                                             gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgg                                                     <  2:10556589/150‑1 (MQ=255)
                                                                                              cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgc                                                    >  1:7299063/1‑150 (MQ=255)
                                                                                              cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCGGAATCCCTGCGAAGgc                                                    >  2:7425960/1‑150 (MQ=255)
                                                                                               gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcc                                                   <  2:7589389/150‑1 (MQ=255)
                                                                                                aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                  <  2:820986/150‑1 (MQ=255)
                                                                                                  cTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagc                                                >  1:5468370/1‑150 (MQ=255)
                                                                                                   tGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcc                                               <  2:4293318/150‑1 (MQ=255)
                                                                                                   tGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcc                                               >  1:2323134/1‑150 (MQ=255)
                                                                                                    gCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcca                                              <  2:10315548/150‑1 (MQ=255)
                                                                                                    gCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcca                                              >  2:1186627/1‑150 (MQ=255)
                                                                                                      cGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                            <  1:887179/150‑1 (MQ=255)
                                                                                                       gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgc                                           <  2:9480165/150‑1 (MQ=255)
                                                                                                       gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgc                                           <  1:7828952/150‑1 (MQ=255)
                                                                                                        aaGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGACGCTGAATCCCTGCGAAGGCCAGCCATgcg                                          >  1:9524589/1‑150 (MQ=255)
                                                                                                        aaGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                          <  1:6768748/150‑1 (MQ=255)
                                                                                                          ggCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCa                                        >  1:2485681/1‑150 (MQ=255)
                                                                                                            cgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAcc                                      <  2:2749131/150‑1 (MQ=255)
                                                                                                             gcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                     <  1:7016004/150‑1 (MQ=255)
                                                                                                                    gCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTccc                              <  1:9206397/150‑1 (MQ=255)
                                                                                                                         ccGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGccc                         >  1:4214242/1‑150 (MQ=255)
                                                                                                                           ggCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTg                       >  1:10283262/1‑150 (MQ=255)
                                                                                                                            gCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                      <  2:903462/150‑1 (MQ=255)
                                                                                                                               cGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcc                   >  2:3643115/1‑150 (MQ=255)
                                                                                                                               cGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcc                   >  1:1175259/1‑150 (MQ=255)
                                                                                                                                gTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcct                  <  1:9024847/150‑1 (MQ=255)
                                                                                                                                  gCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCg                <  2:8329762/150‑1 (MQ=255)
                                                                                                                                     gCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCGAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtg             >  2:8788539/1‑150 (MQ=255)
                                                                                                                                           ccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgc       <  2:2556291/150‑1 (MQ=255)
                                                                                                                                            cgtcgtCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgct      <  1:10856952/150‑1 (MQ=255)
                                                                                                                                              tcgtcgACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTgg    >  2:4593846/1‑150 (MQ=255)
                                                                                                                                               cgtcGACGGGCCAGTGTCGCGGGCGGGAATGTGCCTGTCGCACCCAGGTCTTCAAGGCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAAGCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggc   <  2:6786013/150‑1 (MQ=255)
                                                                                                                                                gtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggcc  >  2:10960492/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCC  >  NC_008463/5791791‑5792085

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: