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.99.8% 2975.3 / ‑3.8 598intergenic (+67/‑68)nadE/choEammonia‑dependent NAD(+) synthetase/cholinesterase ChoE
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (284/313);  total (285/313)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 598 total aligned reads displayed.

TGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCTC  >  NC_008463/5791827‑5792118
                                                                                                                                                  |                                                                                                                                                 
 gcgcGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACCACACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAAcc                                                                                                                                             <  1:3051720/151‑1 (MQ=255)
   gcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAg                                                                                                                                           >  2:6997281/1‑151 (MQ=255)
      aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTc                                                                                                                                        <  1:10745877/151‑1 (MQ=255)
      aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTc                                                                                                                                        <  1:1317923/151‑1 (MQ=255)
       ggAAGCCGCGCGGGTGATCGTCGACACCTACCACACGACCCAGCACACGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCt                                                                                                                                       <  2:7869823/151‑1 (MQ=255)
         aaGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTc                                                                                                                                     <  2:1375722/151‑1 (MQ=255)
           gCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCaa                                                                                                                                   <  2:533053/151‑1 (MQ=255)
                gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGGCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcg                                                                                                                              >  1:3526136/1‑151 (MQ=255)
                 cgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCcgc                                                                                                                             >  1:10943920/1‑151 (MQ=255)
                 cgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgc                                                                                                                             <  1:9651558/151‑1 (MQ=255)
                 cgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgc                                                                                                                             <  2:4621292/151‑1 (MQ=255)
                    gTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                                                          <  2:5551933/151‑1 (MQ=255)
                       aTCGTCGACACCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcg                                                                                                                       <  2:4038389/151‑1 (MQ=255)
                        tcgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcgg                                                                                                                      >  2:1888615/1‑151 (MQ=255)
                             gACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                                                 >  1:9194864/1‑151 (MQ=255)
                               cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                               <  1:8155252/151‑1 (MQ=255)
                               cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                               >  1:5339937/1‑151 (MQ=255)
                                  cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                                                            <  1:8583630/151‑1 (MQ=255)
                                   tACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAg                                                                                                           >  1:1031826/1‑151 (MQ=255)
                                      cacaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                                                        <  1:4623440/151‑1 (MQ=255)
                                        cacgACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                                                      <  2:2624135/148‑1 (MQ=255)
                                        caAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                                                      >  2:996760/1‑151 (MQ=255)
                                        caAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                                                      <  2:175193/151‑1 (MQ=255)
                                          aGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGACGCTg                                                                                                    >  2:10154405/1‑151 (MQ=255)
                                          aGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTg                                                                                                    <  2:297181/151‑1 (MQ=255)
                                            aCCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGaa                                                                                                  <  1:7790190/151‑1 (MQ=255)
                                             cccAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAAt                                                                                                 >  2:7149898/1‑151 (MQ=255)
                                               cAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATcc                                                                                               <  2:2133109/151‑1 (MQ=255)
                                                aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                                                              >  2:515538/1‑151 (MQ=255)
                                                aGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                                                              >  1:7452207/1‑151 (MQ=255)
                                                 gCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCt                                                                                             <  2:5901206/151‑1 (MQ=255)
                                                  cacaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                                                            <  1:4206119/151‑1 (MQ=255)
                                                    caAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCg                                                                                          <  1:6135183/151‑1 (MQ=255)
                                                     aaGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGa                                                                                         >  1:178480/1‑151 (MQ=255)
                                                     aaGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGa                                                                                         <  1:7139784/151‑1 (MQ=255)
                                                          cgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCAGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcc                                                                                    <  2:5782116/151‑1 (MQ=255)
                                                           gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                   >  2:8241716/1‑151 (MQ=255)
                                                           gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAAGGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                   >  1:733981/1‑151 (MQ=255)
                                                            aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                                                  >  2:10710650/1‑151 (MQ=255)
                                                             aCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagc                                                                                 <  1:9312922/151‑1 (MQ=255)
                                                               tGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccagcca                                                                               <  2:9609806/151‑1 (MQ=255)
                                                                   gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGACAGTCGTCGACGGGCCTTTGTCGCGGGAGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                                                           <  1:870940/151‑1 (MQ=255)
                                                                     aGGCGCCCTGAGCCTCCCGGCCCGTGCGGCTGGCCCGTCGTCGCCGGGCCTTTGTCGCGGGCGGGAATGTGCCAGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCa                                                                         <  2:542937/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:3407765/1‑151 (MQ=255)
                                                                            cTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtct                                                                  >  2:8023528/1‑151 (MQ=255)
                                                                            cTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtct                                                                  <  2:11194694/151‑1 (MQ=255)
                                                                               aGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTccc                                                               >  1:1133882/1‑151 (MQ=255)
                                                                                gCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTcccc                                                              <  1:4864766/151‑1 (MQ=255)
                                                                                gCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTcccc                                                              >  1:8071359/1‑151 (MQ=255)
                                                                                 ccTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCg                                                             <  2:11384608/151‑1 (MQ=255)
                                                                                  cTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGc                                                            >  2:10949573/1‑151 (MQ=255)
                                                                                   tCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCTAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGcc                                                           >  1:1494022/1‑151 (MQ=255)
                                                                                   tCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGcc                                                           >  1:4008689/1‑151 (MQ=255)
                                                                                       ggCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                                                       <  2:4870151/151‑1 (MQ=255)
                                                                                       ggCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGc                                                       >  1:3293088/1‑151 (MQ=255)
                                                                                         cccGTGCGGCAGGCCCGTCGTCGACGGCCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctc                                                     <  2:8430148/151‑1 (MQ=255)
                                                                                          ccGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcc                                                    >  1:9915270/1‑151 (MQ=255)
                                                                                           cGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcct                                                   <  2:7436634/151‑1 (MQ=255)
                                                                                            gTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGctcctc                                                  >  2:8596715/1‑151 (MQ=255)
                                                                                             tGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCg                                                 <  2:1836681/151‑1 (MQ=255)
                                                                                                 gCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtgt                                             >  1:378377/1‑151 (MQ=255)
                                                                                                   aGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTcc                                           <  1:7593703/151‑1 (MQ=255)
                                                                                                    ggCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCct                                          >  1:10221565/1‑151 (MQ=255)
                                                                                                     gCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctg                                         <  1:6864030/151‑1 (MQ=255)
                                                                                                       ccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgct                                       <  1:11559409/151‑1 (MQ=255)
                                                                                                       ccGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgct                                       <  1:11442566/151‑1 (MQ=255)
                                                                                                         gtcgtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTgg                                     >  2:6009936/1‑151 (MQ=255)
                                                                                                          tcgtcgACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggc                                    >  2:886759/1‑151 (MQ=255)
                                                                                                             tcgaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccgg                                 <  2:9798973/148‑1 (MQ=255)
                                                                                                              cgaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccggc                                >  2:10798791/3‑151 (MQ=255)
                                                                                                                aCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                              <  2:6221372/151‑1 (MQ=255)
                                                                                                                 cGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAg                             >  2:3303918/1‑151 (MQ=255)
                                                                                                                 cGGGCCAGTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGACGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAg                             <  1:7906285/151‑1 (MQ=255)
                                                                                                                  gggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAgg                            >  2:6476687/1‑151 (MQ=255)
                                                                                                                   ggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGc                           <  1:6728026/151‑1 (MQ=255)
                                                                                                                   ggCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGc                           >  1:4985683/1‑151 (MQ=255)
                                                                                                                    gCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGcc                          >  1:10096236/1‑151 (MQ=255)
                                                                                                                      cTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTg                        <  1:6331214/151‑1 (MQ=255)
                                                                                                                      cTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTg                        >  2:6718627/1‑151 (MQ=255)
                                                                                                                       tttGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGc                       <  1:11118892/151‑1 (MQ=255)
                                                                                                                        ttGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCACGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCg                      <  2:3865252/151‑1 (MQ=255)
                                                                                                                         tGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCgg                     >  1:4393929/1‑151 (MQ=255)
                                                                                                                           tCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGac                   >  2:3769316/1‑151 (MQ=255)
                                                                                                                             gcgggcggGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacac                 <  1:2085459/151‑1 (MQ=255)
                                                                                                                              cgggcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacaca                <  1:6636913/151‑1 (MQ=255)
                                                                                                                              cgggcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacaca                <  1:257881/151‑1 (MQ=255)
                                                                                                                               gggcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGTAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacacac               <  2:9204103/151‑1 (MQ=255)
                                                                                                                               gggcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacacac               <  2:2325613/151‑1 (MQ=255)
                                                                                                                                  cgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTc            >  1:5596769/1‑151 (MQ=255)
                                                                                                                                   gggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTcc           >  2:10525513/1‑151 (MQ=255)
                                                                                                                                    ggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTccc          >  1:9022946/1‑151 (MQ=255)
                                                                                                                                    ggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACACGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCAGGCCAGGCCTGCGGACACACGTccc          >  2:2512468/1‑151 (MQ=255)
                                                                                                                                     gattgtgCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCGGAATCCCTGCGAAGGCCAGCCATGCGCACCCGACTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTcccc         <  1:10494168/148‑1 (MQ=255)
                                                                                                                                     gAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTcccc         >  2:4064467/1‑151 (MQ=255)
                                                                                                                                       aTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgc       >  1:5719887/1‑151 (MQ=255)
                                                                                                                                         gtgCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgctg     >  2:4322528/1‑151 (MQ=255)
                                                                                                                                           gCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgctgct   <  2:312739/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCTC  >  NC_008463/5791827‑5792118

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: