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

Only 100 of 594 total aligned reads displayed.

TATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCG  >  NC_008463/5791790‑5792086
                                                                                                                                                    |                                                                                                                                                    
tATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGAc                                                                                                                                                    <  2:528811/151‑1 (MQ=255)
  tGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCCCAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACgg                                                                                                                                                  <  1:7406925/151‑1 (MQ=255)
  tGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACgg                                                                                                                                                  >  1:2065489/1‑151 (MQ=255)
   gCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACggg                                                                                                                                                 >  1:10847203/1‑151 (MQ=255)
   gCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACggg                                                                                                                                                 <  1:11051442/151‑1 (MQ=255)
     cGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGcc                                                                                                                                               >  1:485978/1‑151 (MQ=255)
      gagaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCt                                                                                                                                              >  1:8304428/1‑151 (MQ=255)
      gagaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCt                                                                                                                                              >  2:4714769/1‑151 (MQ=255)
      gagaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCt                                                                                                                                              >  1:6523227/1‑151 (MQ=255)
       agaTCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCtt                                                                                                                                             <  1:2538102/151‑1 (MQ=255)
            gACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcg                                                                                                                                        <  2:1351276/151‑1 (MQ=255)
             aCGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTcgc                                                                                                                                       <  2:4365737/151‑1 (MQ=255)
               gCCTTCCAGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAACACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCACGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCgcgg                                                                                                                                     <  1:5545153/151‑1 (MQ=255)
                ccTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCgcggg                                                                                                                                    <  1:11409167/151‑1 (MQ=255)
                   tCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACCCGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCgcgggcgg                                                                                                                                 <  1:358176/151‑1 (MQ=255)
                   tCCTGCACGGCCAGCCGCTGCGCGAGGAAGCAGCGCGGGACATCTGAGACACATACCACAAGACCCAGCACATGAGCGTTCTGCCGAAGGCGCCCTGAGCCTCCCGGCGCGGGCGGCAGGCCCGTCGTCGACGGGGGTTTGTCgcgggcgg                                                                                                                                 >  1:9615008/1‑151 (MQ=255)
                    ccTGCACGGCCAGCCGCTGCGCGAGGAAGACGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGcgggcggg                                                                                                                                >  2:8046777/1‑151 (MQ=255)
                     cTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGa                                                                                                                               <  2:1469778/151‑1 (MQ=255)
                     cTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGa                                                                                                                               >  2:6789196/1‑151 (MQ=255)
                      tGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGaa                                                                                                                              >  1:3263307/1‑151 (MQ=255)
                         aCGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACGGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgt                                                                                                                           <  2:9523250/151‑1 (MQ=255)
                          cGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAtgtg                                                                                                                          <  2:3079388/151‑1 (MQ=255)
                             ccAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCt                                                                                                                       <  2:10118617/151‑1 (MQ=255)
                                  cGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTATGTCGCGGGCGGGAATGTGCCTGTcgc                                                                                                                  >  2:7803111/1‑149 (MQ=255)
                                       cgcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCa                                                                                                            <  1:1812722/151‑1 (MQ=255)
                                         cgAGGACGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAgg                                                                                                          >  2:9123733/1‑151 (MQ=255)
                                          gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGt                                                                                                         >  1:3038297/1‑151 (MQ=255)
                                          gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                         <  2:2198157/151‑1 (MQ=255)
                                          gAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGt                                                                                                         >  1:1670935/1‑151 (MQ=255)
                                            ggACGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCGTTTGTCGCGGGCGGGAATGTGCCTGTCGCACCCAGGTCt                                                                                                       >  1:6941386/1‑151 (MQ=255)
                                            ggAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCt                                                                                                       <  1:3365245/151‑1 (MQ=255)
                                            ggAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCt                                                                                                       <  2:8459596/151‑1 (MQ=255)
                                                    cgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                               <  2:4940647/151‑1 (MQ=255)
                                                     gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcg                                                                                              >  2:7247996/1‑151 (MQ=255)
                                                      cgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgc                                                                                             <  2:538603/151‑1 (MQ=255)
                                                       gggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCACCCAGGTCTTCCAGTCcgcg                                                                                            >  2:6523310/1‑151 (MQ=255)
                                                       gggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                            >  1:7917477/1‑151 (MQ=255)
                                                       gggTGATCGGCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                            <  1:3974766/151‑1 (MQ=255)
                                                         gTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                          <  2:9950462/151‑1 (MQ=255)
                                                          tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcgg                                                                                         <  2:3251055/151‑1 (MQ=255)
                                                          tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcgg                                                                                         <  1:10258536/151‑1 (MQ=255)
                                                           gATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                        <  2:8824646/151‑1 (MQ=255)
                                                             tcgtcgACACCTACCACAAGACCCAGCACACGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                      <  2:4255139/151‑1 (MQ=255)
                                                             tcgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                      <  1:6145431/151‑1 (MQ=255)
                                                               gtcgACCCCTACCACACGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                    <  2:5994627/151‑1 (MQ=255)
                                                                tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                   <  2:2419887/151‑1 (MQ=255)
                                                                tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                   >  2:7439475/1‑151 (MQ=255)
                                                                 cgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcg                                                                                  <  2:7928532/151‑1 (MQ=255)
                                                                  gACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                 >  1:9194864/1‑151 (MQ=255)
                                                                   acacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCt                                                                                <  1:9028190/151‑1 (MQ=255)
                                                                   acacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCt                                                                                <  1:5885023/151‑1 (MQ=255)
                                                                    cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                               >  1:5339937/1‑151 (MQ=255)
                                                                     cccTACCACACGACCCCGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTag                                                                              <  2:708920/150‑2 (MQ=255)
                                                                     acCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTag                                                                              >  2:3908518/1‑150 (MQ=255)
                                                                       cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                            <  2:545932/151‑1 (MQ=255)
                                                                       cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                            <  1:8583630/151‑1 (MQ=255)
                                                                           cccaAGACCCCGCACTAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGACTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                        <  1:10128765/149‑1 (MQ=255)
                                                                             cacgACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                      <  2:2624135/148‑1 (MQ=255)
                                                                             caAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                      <  2:175193/151‑1 (MQ=255)
                                                                              aaGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCt                                                                     <  1:8772134/151‑1 (MQ=255)
                                                                                  cccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAAt                                                                 <  1:6791826/151‑1 (MQ=255)
                                                                                     aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                              <  2:3347992/151‑1 (MQ=255)
                                                                                     aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                              >  1:7478731/1‑151 (MQ=255)
                                                                                     aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAAAccc                                                              <  1:4955957/151‑1 (MQ=255)
                                                                                      gCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCTACCAGGTCTTCAAGTCCGCGCGGCGGGGCGCTAGAAGGAAGCTGACTCCCt                                                             >  2:11356661/1‑151 (MQ=255)
                                                                                      gCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCt                                                             >  2:4013310/1‑151 (MQ=255)
                                                                                       cacaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                            <  1:4206119/151‑1 (MQ=255)
                                                                                        acaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGc                                                           <  1:1350605/151‑1 (MQ=255)
                                                                                           cgcgcGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCACGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGaa                                                        <  2:10179195/150‑1 (MQ=255)
                                                                                           aGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGaa                                                        <  2:10785575/151‑1 (MQ=255)
                                                                                             cgcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgg                                                      <  1:5703890/151‑1 (MQ=255)
                                                                                              gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgc                                                     >  1:6392841/1‑151 (MQ=255)
                                                                                               cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcc                                                    >  1:7763212/1‑151 (MQ=255)
                                                                                               cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCAGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcc                                                    <  2:5782116/151‑1 (MQ=255)
                                                                                                gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                   >  2:4062096/1‑151 (MQ=255)
                                                                                                gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAGGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                   >  1:733981/1‑151 (MQ=255)
                                                                                                  aCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagc                                                 <  1:9315158/151‑1 (MQ=255)
                                                                                                   cTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcc                                                >  2:10636658/1‑151 (MQ=255)
                                                                                                      ccGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATg                                             <  1:9890575/151‑1 (MQ=255)
                                                                                                       cGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgc                                            >  2:3715493/1‑151 (MQ=255)
                                                                                                       cGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgc                                            >  1:6710986/1‑151 (MQ=255)
                                                                                                         aaGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCACCCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGTCCAGCCATgcgc                                          >  1:7167425/1‑151 (MQ=255)
                                                                                                           ggCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAc                                        >  1:4443400/1‑151 (MQ=255)
                                                                                                            gcgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAcc                                       <  2:2963619/151‑1 (MQ=255)
                                                                                                             cgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                      <  2:11588473/151‑1 (MQ=255)
                                                                                                              gcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCg                                     >  1:2295157/1‑151 (MQ=255)
                                                                                                                ccTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtc                                   >  2:7642797/1‑151 (MQ=255)
                                                                                                                 cTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtct                                  <  2:11194694/151‑1 (MQ=255)
                                                                                                                 cTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtct                                  >  1:8098030/1‑151 (MQ=255)
                                                                                                                  tGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtctc                                 <  2:8250597/151‑1 (MQ=255)
                                                                                                                    aGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTccc                               >  1:1133882/1‑151 (MQ=255)
                                                                                                                       cTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGc                            >  2:10949573/1‑151 (MQ=255)
                                                                                                                            ggCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                       >  2:5850543/1‑151 (MQ=255)
                                                                                                                            ggCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                       >  1:10656367/1‑151 (MQ=255)
                                                                                                                               cactGCGGCAGGCCCGGCGACGACGGGCCATTGACGCCGGCGGGAATGTGCCTCTCACAACCAGCTCTTCAAGTACGCTCGCCGGCGCGATAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCACCGCCCAGctcc                    >  2:9678624/4‑151 (MQ=255)
                                                                                                                                  tGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCg                 <  2:1836877/151‑1 (MQ=255)
                                                                                                                                      gCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCAGGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtgt             >  1:378342/1‑151 (MQ=255)
                                                                                                                                         ggCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCACCCCGGTCTTCACGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGCAGGCCAGCCATGCGCACCCGTCGCCCCGCCCTGCTCCTCGGTGTCct          <  1:4635513/151‑1 (MQ=255)
                                                                                                                                            ccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgct       <  1:11559409/151‑1 (MQ=255)
                                                                                                                                             cgtcgtCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgctg      >  2:1897575/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATGCCGAGATCGACGCCTTCCTGCACGGCCAGCCGCTGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCG  >  NC_008463/5791790‑5792086

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: