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.5% 2862.4 / 2.5 581intergenic (+67/‑68)nadE/choEammonia‑dependent NAD(+) synthetase/cholinesterase ChoE
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (298/280);  total (300/281)
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 581 total aligned reads displayed.

TGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCTC  >  NC_008463/5791827‑5792118
                                                                                                                                                  |                                                                                                                                                 
  cgcgAGGACGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCa                                                                                                                                            >  2:10035813/1‑151 (MQ=255)
  cgcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCa                                                                                                                                            >  2:3670045/1‑151 (MQ=255)
  cgcgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCa                                                                                                                                            >  1:2366386/1‑151 (MQ=255)
    cgAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAgg                                                                                                                                          <  2:5879928/151‑1 (MQ=255)
      aGGAAGCCGCGCGGGTGATCGTCGACACCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTc                                                                                                                                        <  2:6216072/151‑1 (MQ=255)
       ggAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCt                                                                                                                                       >  1:3467677/1‑151 (MQ=255)
        gAAGCCGCGCGGGTGATCGTCGACACCTACCACCAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCtt                                                                                                                                      <  1:3817481/151‑1 (MQ=255)
             cgcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGt                                                                                                                                 <  1:6335094/151‑1 (MQ=255)
              gcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTc                                                                                                                                <  1:9732914/151‑1 (MQ=255)
              gcgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTc                                                                                                                                <  2:7648937/151‑1 (MQ=255)
               cgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                                                               >  1:8453128/1‑151 (MQ=255)
               cgcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTcc                                                                                                                               <  1:7763365/151‑1 (MQ=255)
                gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcg                                                                                                                              >  2:3921670/1‑151 (MQ=255)
                gcgGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcg                                                                                                                              >  1:3208831/1‑151 (MQ=255)
                  gggTGGTCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                                                            >  1:4025438/1‑151 (MQ=255)
                  gggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTAGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcg                                                                                                                            >  2:974795/1‑151 (MQ=255)
                   ggTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgc                                                                                                                           <  1:267132/151‑1 (MQ=255)
                    gTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                                                          <  1:8663638/151‑1 (MQ=255)
                    gTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCcgcgcg                                                                                                                          >  2:5181178/1‑151 (MQ=255)
                     tGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcgg                                                                                                                         <  1:5920654/151‑1 (MQ=255)
                      gATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggc                                                                                                                        >  1:6381875/1‑151 (MQ=255)
                       aTCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcg                                                                                                                       >  1:9424687/1‑151 (MQ=255)
                         cgtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggc                                                                                                                     >  2:4358677/1‑151 (MQ=255)
                          gtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGACTGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                                                    >  2:8884372/1‑151 (MQ=255)
                          gtcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCgcggcggcg                                                                                                                    >  1:1611215/1‑151 (MQ=255)
                           tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                                                   <  1:8351237/151‑1 (MQ=255)
                           tcgACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgc                                                                                                                   <  1:2099649/151‑1 (MQ=255)
                            cgACCCCTACCACAAGACCCAGCACAAGCGCGAACTGCCGACGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCGTTGTCGCGGGCGGGAATGTGCCTGGCGCACCCAGGGCTTCAAGTCCGCGCGGCGgcgcg                                                                                                                  >  2:10512939/1‑151 (MQ=255)
                             gACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGgcgcgc                                                                                                                 >  2:4893894/1‑151 (MQ=255)
                               cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                               <  1:9067638/151‑1 (MQ=255)
                               cacCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTa                                                                                                               <  1:6793957/151‑1 (MQ=255)
                                 ccTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTaga                                                                                                             >  2:7937205/1‑149 (MQ=255)
                                  cTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTagaa                                                                                                            >  1:4343881/1‑151 (MQ=255)
                                     ccACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGa                                                                                                         <  1:1164802/151‑1 (MQ=255)
                                      cacaAGACCCAGCACAAGCGCGAACTGCCGCAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                                                        <  2:3183596/151‑1 (MQ=255)
                                      cacaAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGaa                                                                                                        <  2:2179093/151‑1 (MQ=255)
                                        caAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGc                                                                                                      >  2:8159846/1‑151 (MQ=255)
                                           gACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGa                                                                                                   >  1:9882824/1‑151 (MQ=255)
                                                aGCACAAGCGCGAACTGCCGACGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCAGTCGCAACCCGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATccc                                                                                              <  1:4487174/151‑1 (MQ=255)
                                                 gCACAAGCGCGAACTGCCGAAGGCGACCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTGGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCt                                                                                             <  2:9132451/151‑1 (MQ=255)
                                                  cacaAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTg                                                                                            >  1:1522915/1‑151 (MQ=255)
                                                   acaAGCGCGAACAGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGc                                                                                           >  2:8424389/1‑151 (MQ=255)
                                                      aGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGaa                                                                                        <  1:328219/151‑1 (MQ=255)
                                                        cgcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAgg                                                                                      <  1:919479/151‑1 (MQ=255)
                                                         gcgAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgc                                                                                     <  2:192546/151‑1 (MQ=255)
                                                           gAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgcca                                                                                   >  2:4754168/1‑151 (MQ=255)
                                                            aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                                                  <  2:9188712/151‑1 (MQ=255)
                                                            aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                                                  <  1:10329310/151‑1 (MQ=255)
                                                            aaCTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGgccag                                                                                  <  1:10664736/151‑1 (MQ=255)
                                                              cTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGCCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCAACAATTACGCGCGGCGGCGCGCTAGAAGGACGCTGAATCCCAGCGACGgccagcc                                                                                >  1:5599971/1‑151 (MQ=255)
                                                                   gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                                                           <  1:7450666/151‑1 (MQ=255)
                                                                   gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                                                           >  2:6963928/1‑151 (MQ=255)
                                                                   gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATgcg                                                                           >  2:7363030/1‑151 (MQ=255)
                                                                   gAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCGAGAAGGAAGCTGACTCCCTGCGAAGGCCAGCCATgcg                                                                           >  1:3036238/1‑151 (MQ=255)
                                                                        cgcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCATTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTACAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCAccc                                                                      <  1:10865521/151‑1 (MQ=255)
                                                                         gcCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCg                                                                     <  1:6887270/151‑1 (MQ=255)
                                                                           ccTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtc                                                                   >  2:508936/1‑151 (MQ=255)
                                                                             tGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGtctc                                                                 <  2:7196536/151‑1 (MQ=255)
                                                                                gCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTcccc                                                              <  2:10203824/151‑1 (MQ=255)
                                                                                  cTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGc                                                            >  1:3350441/1‑151 (MQ=255)
                                                                                  cTCCCGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGc                                                            >  1:11132235/1‑151 (MQ=255)
                                                                                    cccGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGccc                                                          >  1:1618757/1‑151 (MQ=255)
                                                                                    cccGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGccc                                                          >  1:2815945/1‑151 (MQ=255)
                                                                                      cGGCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTg                                                        >  1:1320822/1‑151 (MQ=255)
                                                                                      cGGCCCGTGCGGCAGGCCCGTCGACGACGGGCCTTTGTCGCGGGCGGGAAAGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTg                                                        >  2:6705877/1‑151 (MQ=255)
                                                                                        gCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGct                                                      <  2:1764358/151‑1 (MQ=255)
                                                                                        gCCCGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGct                                                      <  1:2567131/151‑1 (MQ=255)
                                                                                           cGTGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTAGCAACCAGGTCTTCAATTCAGCGCGGCTGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGACATGCGCACCCGTCTCCCCGCCCAGctcct                                                   >  1:5335254/1‑151 (MQ=255)
                                                                                             tGCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCg                                                 >  1:7148810/1‑151 (MQ=255)
                                                                                              gCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCgg                                                <  2:6522344/151‑1 (MQ=255)
                                                                                              gCGGCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCgg                                                <  2:2932958/151‑1 (MQ=255)
                                                                                                ggCAGGCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGgtg                                              <  1:8714923/151‑1 (MQ=255)
                                                                                                    ggCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCct                                          >  1:9602973/1‑151 (MQ=255)
                                                                                                     gCCCGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGACGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctg                                         >  2:3852284/1‑151 (MQ=255)
                                                                                                       ccGTCGTCGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCctgct                                       >  2:5651357/1‑151 (MQ=255)
                                                                                                           cgtcGACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggcc                                   >  1:10771170/1‑151 (MQ=255)
                                                                                                             tcgaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCCTGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccgg                                 >  2:2686852/4‑151 (MQ=255)
                                                                                                               gaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccggcc                               >  1:8916981/2‑151 (MQ=255)
                                                                                                               gaCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTggccggcc                               <  2:11036371/150‑1 (MQ=255)
                                                                                                                aCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                              >  2:5341541/1‑151 (MQ=255)
                                                                                                                aCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                              >  2:8760187/1‑151 (MQ=255)
                                                                                                                aCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCa                              >  1:11057546/1‑151 (MQ=255)
                                                                                                                aCGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGATGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCATCGGTGTCCTGCTGGCCGGCCa                              >  1:10005344/1‑151 (MQ=255)
                                                                                                                    gCCTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGcc                          <  1:6149911/151‑1 (MQ=255)
                                                                                                                     ccTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCt                         >  2:2089254/1‑151 (MQ=255)
                                                                                                                      cTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTg                        >  2:9475188/1‑151 (MQ=255)
                                                                                                                      cTTTGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTg                        >  1:3738978/1‑151 (MQ=255)
                                                                                                                        ttGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCg                      <  2:11097494/151‑1 (MQ=255)
                                                                                                                         tGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCgg                     >  2:2329235/1‑151 (MQ=255)
                                                                                                                         tGTCGCGGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCgg                     >  2:1290957/1‑151 (MQ=255)
                                                                                                                            cgcgGGCGGGAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGaca                  >  2:4235498/1‑151 (MQ=255)
                                                                                                                              cgggcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGacaca                <  2:6616763/151‑1 (MQ=255)
                                                                                                                                 gcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGt             >  2:10781523/1‑151 (MQ=255)
                                                                                                                                 gcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGt             <  2:3526839/151‑1 (MQ=255)
                                                                                                                                 gcgggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGt             <  1:5106728/151‑1 (MQ=255)
                                                                                                                                   gggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTcc           <  2:6363312/151‑1 (MQ=255)
                                                                                                                                   gggAATGTGCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTcc           >  1:4758770/1‑151 (MQ=255)
                                                                                                                                     gAATGTGCCTGTCGCACCCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGACGGAAGCTTAATCCCTGCGATGCCCAGCCTTGCGCACCCGTCACCCCGCCCGGCTCCTGGGTGTCCTGGTGGCCGGCCTGGCCTGCGGACACACGTcccc         <  1:1956833/151‑1 (MQ=255)
                                                                                                                                        tgtgCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgct      <  1:8217688/151‑1 (MQ=255)
                                                                                                                                        tgtgCCTGTCGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAGAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCgct      >  2:6026159/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCGCGAGGAAGCCGCGCGGGTGATCGTCGACACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGAAGGCGCCCTGAGCCTCCCGGCCCGTGCGGCAGGCCCGTCGTCAACGGGCCTTTGTCGCGGGCGGGAATGTGCCTGTCAGCAACCAGGTCTTCAAGTCCGCGCGGCGGCGCGCTAAAAGGAAGCTGAATCCCTGCGAAGGCCAGCCATGCGCACCCGTCTCCCCGCCCTGCTCCTCGGTGTCCTGCTGGCCGGCCAGGCCTGCGGACACACGTCCCCGCTGCTC  >  NC_008463/5791827‑5792118

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