Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,851,821 G→A G274D (GGC→GAC)  hypB → hydrogenase isoenzymes nickel incorporation protein HypB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,8210GA98.3% 182.1 / NA 60G274D (GGC→GAC) hypBhydrogenase isoenzymes nickel incorporation protein HypB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (27/32);  total (27/33)
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

CGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGAC  >  NC_000913/2851576‑2852061
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
cGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTTCGAAGGGGAAGACAAACCACCGAAATATCCGCATATGTTTTCCGCCCCCTCGCTGATGCTGCTCAACAAAGATGACCTTTTGCCGTATCTCCACTTT‑GACGTTGAAAAAAGCATCGCCTGC‑GCCCGCGAAGTCAATCCCGAAATTGAAATCATACTTATTTCCGCCACCAGcgacga                                                                                                                                                                                                                                               >  1:553735/1‑249 (MQ=255)
                 gCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTg                                                                                                                                                                                                                            <  1:109819/251‑1 (MQ=255)
                    cGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGTTACCGGAGGGGAAGAAAAAACACTGAAATATCCCCCTATGTTTTCCCCCGCCCCCCTGGTGCTGCTCCACCAAGATGACCTGTTGCCGTATCTCCACTTTTGA‑GTTGAAAAATTCCTCCCCTCC‑GCCCCCGAAGTCCATCCCGAAAATGAAATCATCCTTTTTTTCGCCCCCCCCGACGAAGGGATGGACCCGTGGCTGAAc                                                                                                                                                                                                                         >  1:644192/1‑251 (MQ=255)
                      gCCATCTTCGATCTCGGTGAAAAACACAAAGAGGCGGTGCTTTCCGTGACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCGCGCTGATGCTGCTCAACAAAGTTGACTTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCTTCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAActg                                                                                                                                                                                                                       <  1:449601/251‑1 (MQ=255)
                               gATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGt                                                                                                                                                                                                                                 <  1:365545/232‑1 (MQ=255)
                                 catcGGTGAACAACACAAAGGGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTCTCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTACCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCTACCAGCGACGAAGGGATGGACCAGTGGCTGAACGGGCTGGAGacac                                                                                                                                                                                                            <  2:258282/249‑1 (MQ=255)
                                    cGGTGAAAAACACAAAGTGTCGGTGCTTTTCTGTACCGAAGTTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGTCCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGaca                                                                                                                                                                                                             <  2:547994/247‑1 (MQ=255)
                                         aaaaaCACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCCCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGAtgt                                                                                                                                                                                                    <  1:679509/251‑1 (MQ=255)
                                                      ggCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTAGTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTc                                                                                                                                                                                        <  1:260714/250‑1 (MQ=255)
                                                       gCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTGccc                                                                                                                                                                                      <  1:457122/251‑1 (MQ=255)
                                                          gTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGc                                                                                                                                                                                   >  1:601573/1‑251 (MQ=255)
                                                                        cGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCAt                                                                                                                                                                     >  1:502508/1‑251 (MQ=255)
                                                                                   aCAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAAcc                                                                                                                                                          >  1:335629/1‑251 (MQ=255)
                                                                                   aCAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAAcc                                                                                                                                                          <  1:595640/251‑1 (MQ=255)
                                                                                     aacccACGGAAATAGCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAg                                                                                                                                                        <  2:33321/248‑1 (MQ=255)
                                                                                        ccACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTAGCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCTCAAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCg                                                                                                                                                     <  1:628215/251‑1 (MQ=255)
                                                                                         cACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGACATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACCGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGCCGGCCACCAgg                                                                                                                                                       >  2:352217/1‑248 (MQ=255)
                                                                                              aaaTATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGAGAAAgcc                                                                                                                                               >  2:242531/1‑249 (MQ=255)
                                                                                                   tCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                                                                          >  1:401948/1‑251 (MQ=255)
                                                                                                    ccGCATATGTTTGCCGCCGCCGCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCATGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGt                                                                                                                                          <  2:644192/250‑1 (MQ=255)
                                                                                                        atatGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCg                                                                                                                                     >  1:104840/1‑251 (MQ=255)
                                                                                                         tatGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCgg                                                                                                                                    >  1:476248/1‑251 (MQ=255)
                                                                                                          atGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGc                                                                                                                                   >  1:388066/1‑251 (MQ=255)
                                                                                                             tttGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGACGACGTCTGCGGCAtt                                                                                                                                >  1:534282/1‑251 (MQ=255)
                                                                                                                gctgccgccTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCAGTGACGGCAACCAGGCGAAAGTCGCCGTCTGCGGCATTCa                                                                                                                              <  2:403061/247‑1 (MQ=255)
                                                                                                                 ccgccgccTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAgc                                                                                                                            <  1:374114/251‑1 (MQ=255)
                                                                                                                    ccgctcatcTGACGCTTCTCAACAAAGTTAACCTGTGGCCGTAGCTCTACATT‑TACTTTGAAAAGTTCATCGCCTGC‑GCCCGAGATGTCAATCCAGAAATGGAAATCATCCTTATTTCCGCCACCAGCGATGAAAGGATGGACTAGTGGCTAAACTGGCTGGAGACACAGCGATGTGCAGAGGCGTTCCCGGCTAGCTCCGCTCCAGTGACGGCAACCAGTCGAAAGTCGACGTCTGCGGCATTCAGcgcg                                                                                                                         <  2:74810/243‑1 (MQ=255)
                                                                                                                      gcctctcggATGCTGCTCAACAAAGTTGTCCTTTCGCTGTATCTCAACTTT‑GACGTGGAAAAGTGCATTGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCTTTTAGCGCGAt                                                                                                                       <  2:581985/243‑1 (MQ=255)
                                                                                                                       ccTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATg                                                                                                                      >  1:650770/1‑251 (MQ=255)
                                                                                                                              gATGCTGCTCAACAAAGTTGCCCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTGCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGAGCTATGTCGAttt                                                                                                               <  2:551257/251‑1 (MQ=255)
                                                                                                                              gATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAttt                                                                                                               >  1:188331/1‑251 (MQ=255)
                                                                                                                                  ctgctCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCg                                                                                                                                                                           <  1:391673/187‑1 (MQ=255)
                                                                                                                                    gctCAAAAAAGGTGACCTGATGCAGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGAGAAGTCACTCCAGAAATTGAAATCATCCTTATTTCCTCCACCAGCGACGAAGGGATGGACCAGTGGCTGCACTGGCAGGAGACCAAGAGCTGGGCATAGGCGTTCCCGGACAGATCCGCACCATTGACGGCGACCAGGCGAAAGGAGACGTCTGAGGCATTCAGCGCGCTGTAGATTGAACggt                                                                                                         >  2:546289/1‑249 (MQ=255)
                                                                                                                                           aaaGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGc                                                                                                  >  1:138797/1‑251 (MQ=255)
                                                                                                                                                      gTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCCTCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGTCTGACCTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGTCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGAt                                                                                         <  1:499899/249‑1 (MQ=255)
                                                                                                                                                      gTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGAt                                                                                         <  1:352217/249‑1 (MQ=255)
                                                                                                                                                                 ccAACTTT‑GACGTTGAAAAGTTCATCGC‑TGTTGCCCGCGAAGGCAATCCAGAAAGTGAAATCATCCTTATTCCCGCCACCCGCGACGATGGGATGGACCAGTGGCTTAACTGGCTGGAGACCCAGCGGTGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGCCGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGACAACGGTCAGc                                                                            <  2:406508/250‑1 (MQ=255)
                                                                                                                                                                  cAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTc                                                                               <  1:138641/247‑1 (MQ=255)
                                                                                                                                                                      ttt‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCcgcgc                                                                       >  1:66676/1‑251 (MQ=255)
                                                                                                                                                                           actttGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTgg                                                                   <  1:313239/248‑1 (MQ=255)
                                                                                                                                                                           aCGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCCGATCCGCACCATTGACGGCAACCAGGCGAAAGACGACGTCTGCGGCCTTCAGCGCGAAGTCGATTTAACGTTAGTCGGCAGATGCGATGCAAACGGTCAGCCGCGCTTgg                                                                   >  2:151055/1‑251 (MQ=255)
                                                                                                                                                                            cGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTggg                                                                  <  1:584572/251‑1 (MQ=255)
                                                                                                                                                                              ttGAAAAGGGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCCTATTTCCGCCACCAGCGGCGCATGGATGGACCAGTGGATGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGGCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGc                                                                 <  2:661715/250‑1 (MQ=255)
                                                                                                                                                                                aaaaaGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGGGGCGGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAg                                                              <  2:502508/250‑1 (MQ=255)
                                                                                                                                                                                         aTCGCCGGC‑GCCCGCGAAGTTAATCCAGAAAGTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTg                                                     <  1:304162/251‑1 (MQ=255)
                                                                                                                                                                                            gCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCCATGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGCTCCGCACCATTGACGGCAACCAGGCGACAGTCGCCGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTCGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTAGGtcttggt                                                   >  2:635699/1‑244 (MQ=255)
                                                                                                                                                                                                     ccGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTgg                                          >  1:104262/1‑251 (MQ=255)
                                                                                                                                                                                                         gAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCttt                                      <  1:242531/251‑1 (MQ=255)
                                                                                                                                                                                                           aGTCAATCCAGAAATGGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACGGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGc                                    <  2:650770/251‑1 (MQ=255)
                                                                                                                                                                                                               aaTCCAGAAATTGAAATCATCCTTATTTACGCCACCAGCGACGAAGGGATGAACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGAAACCAGGCGACAGGCGACGTCTGCGGACTTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTAAACGTTGGCTTTGCCATg                                >  2:240307/1‑251 (MQ=255)
                                                                                                                                                                                                                    aGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTTAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGt                           >  1:129349/1‑251 (MQ=255)
                                                                                                                                                                                                                          tGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGCTGTGCATAGGCGTTCTCGTCCAGATCCGCACCATTGACGGCAACCAGGCGACAGTCGACGTCTGCGGCATTCAGCGAGATGTCGCTTTACCGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCTGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGtaatta                      >  1:656368/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ttATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGa            >  1:29236/1‑249 (MQ=255)
                                                                                                                                                                                                                                     ttATTTACGCCACCAGCGACGAAGGGATGGCACAGTGGCAGAACTGGCTGGAGACACAGCGATGTGAAGAGGCGTTCCCGGCCAAATCCGCACCATTGACGGCAACAAGGCGAACGTCGACGTCTGCGGAATTCAGCGCGATGGCGATTTAACATTAGACGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGGGGGTACAGGAACACGGTGGCTTGGCCATTGGCGTAATTAATGAAGcctaa           >  2:538625/1‑247 (MQ=255)
                                                                                                                                                                                                                                      tATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATAGTCGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTa                                                  <  1:454182/210‑1 (MQ=255)
                                                                                                                                                                                                                                       aTTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCa        <  1:549402/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ttCCGCCCCCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCAc       <  1:350559/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ttCCGCCACCAGCGACGAAGGGATGGCCCAGTGGCTGAACTGGCTGGAGACACAGCGATGTTCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCGGCTGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCAc       <  2:388066/250‑1 (MQ=255)
                                                                                                                                                                                                                                          tCCGCCACCCGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCCGGCGAACGTCGACGTCGGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCTGCAGCTGCGATGACAACTGTCAGCCGCGCGTGGGCCAGTGGGTAATGGTACACGTTGGCTTTTCCATGAGCGTAATTACGTACGCCGAAGCAcgc     >  2:279403/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGac  >  1:494639/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ccaccaGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGac  >  1:152223/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
CGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTT‑GACGTTGAAAAGTGCATCGCCTGC‑GCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGAC  >  NC_000913/2851576‑2852061

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑