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,8210GA97.8% 135.5 / ‑4.2 47G274D (GGC→GAC) hypBhydrogenase isoenzymes nickel incorporation protein HypB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/17);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

AACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTC  >  NC_000913/2851583‑2852067
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
aaCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTGGACGTTGAAAAGTGCAGCGCCTGCGCCCGCGAAGTCGATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGa                                                                                                                                                                                                                                            <  1:226114/251‑1 (MQ=255)
 aCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGAc                                                                                                                                                                                                                                           <  1:37538/251‑1 (MQ=255)
              ggCCATCTTCGCTCTCGGTGACAAACACAAAGTTGCGGTGTTTTCCGTTTCCCAAGGGTAATACAAACCACTGAAATATCCGCATATTTTTGCCGCCGCCTCTCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTTACGTTGAAAAGTGCAGCGCCTGCGCCCGCGAAGTCAATCCAGACATTGAAATCATCCTGATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTTAAct                                                                                                                                                                                                                              <  2:468295/251‑1 (MQ=255)
               gCCAGCTTCGAGCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGa                                                                                                                                                                                                                                            <  1:212063/236‑1 (MQ=255)
                  aGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGCGCTTTCCGTTTACGAAGGTGGAGAAAAAACACTTAAATATCCACATATGTTTTCCGCCGCCTCGCTGATGCTGTTCAAAAAAGTTGAACTGTTTCCGTATCTCAACTTTTAAGTTGAAAAAGGCATCGCCTGCGCCCGCGGAGTCAAACCAGAAATTGAAATCATACTTTTTTCCGCCCCCAGCGAGGAAGGGATGGACCAGTGGCTGGActggct                                                                                                                                                                                                                          >  1:359898/1‑251 (MQ=255)
                             cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCGGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGacac                                                                                                                                                                                                                  <  1:326818/248‑1 (MQ=255)
                                ggAAAAACACAAAGTTGCGGTGCTTCCCGTTACCGAAGGTGAAGACAAACCACTGAAATTACCGCATATGTTTGCTGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCGCAACGTTGACGTTGAAAATTGCATCGCCGGCGCCCGTGAAGTCAATCCAGAAATTGAAATCTTCCTTAGTTCCGCCACCAGCGCCGATGGGATGGAGCAGTGGCTGAACTGGCTGGGGACACAGCGAt                                                                                                                                                                                                            <  2:248728/250‑1 (MQ=255)
                                                               aCCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACCCAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCAcc                                                                                                                                                                             >  2:289748/1‑251 (MQ=255)
                                                                               aaCCACTTAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAgg                                                                                                                                                             <  2:347503/251‑1 (MQ=255)
                                                                                      aaaaTATCCCCATATGTTTGCCGCCGCCTCGCTGATGCTGCCCAACAAAGTTGACCTGTTGCCGTAGCCCAACTTGGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGt                                                                                                                                                      <  2:123734/250‑1 (MQ=255)
                                                                                         atatCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCg                                                                                                                                                    >  1:140705/1‑250 (MQ=255)
                                                                                             ccGCAGATGTTTGCCGCCGCCTCGCTGATGCTTCGCAACAAAGTTGACCTGTTGCCGTATCTCAACTGTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTc                                                                                                                                               <  1:473646/251‑1 (MQ=255)
                                                                                                              cgccTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAACAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCCCCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGcgcg                                                                                                                               >  1:214945/1‑250 (MQ=255)
                                                                                                                     cTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAg                                                                                                                                                                                                    <  2:164965/174‑1 (MQ=255)
                                                                                                                     cTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAg                                                                                                                                                                                                    >  1:164965/1‑174 (MQ=255)
                                                                                                                     cTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAt                                                                                                                       >  1:355340/1‑251 (MQ=255)
                                                                                                                             gctCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAGTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGTTTTAACggt                                                                                                               >  2:80692/1‑249 (MQ=255)
                                                                                                                                  acaaAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCg                                                                                                          >  1:349510/1‑251 (MQ=255)
                                                                                                                                      aGTTGACCTGTTTCCGTATCGCAACGTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTAGTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAg                                                                                                      <  2:419496/251‑1 (MQ=255)
                                                                                                                                      aGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATCCAGCGCGATGTCGATTTAACGTTAGTCGGc                                                                                                        >  1:23899/1‑249 (MQ=255)
                                                                                                                                              tGTTTCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGAAGAAGGGATGGAACAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCACGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGAAACTGCGATg                                                                                              >  2:422755/1‑251 (MQ=255)
                                                                                                                                                    cGTATCTCAAATTTGACGTTGAAAAGTGCATCCCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCCCCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGGGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACg                                                                                        <  2:417676/251‑1 (MQ=255)
                                                                                                                                                     gTATCTCCACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGCCGTCTGCGGCCTTCAGCGCGCTGTCGATTTCACGGTCGTCGGCAGCTGCGATGAAAAc                                                                                         >  2:236172/1‑249 (MQ=255)
                                                                                                                                                        tctcAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTACCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGCTTTAACGTTAGTCGGCAGCTGAGATGAAAACGGTc                                                                                     >  1:431717/1‑250 (MQ=255)
                                                                                                                                                            aaCTTTGACGTTGAAACGTGCATCGCCTGCGAACGCGACGTCAATCCCGAAATTGAAATCATACTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGAGATGTGCATAGGCGTTCCCGGCCAGATCCGCAACATTGTCGGCAACCAGGCGAAAGTCGACGCCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGACGGGAGCGGAGATGaaaaaggcaagc                                                                                  >  1:176139/1‑241 (MQ=255)
                                                                                                                                                              cTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGCTGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCcgcg                                                                              >  1:415463/1‑251 (MQ=255)
                                                                                                                                                               tttGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCcgcgc                                                                             >  2:30319/1‑251 (MQ=255)
                                                                                                                                                                ttGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTGCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCcgcg                                                                              <  2:214288/249‑1 (MQ=255)
                                                                                                                                                                   aCGTTGCAAAGTGAATCGCCTGCACCCGCGAAGTCAATCCAGAAATTGAAATAGTCCTTAATTCCGCCACCAGCGACGAAGGGATGGACAAGTGGCTGAACTGGCTGGAGACACAGCGATGTGAATAAGCGTTCCCGGCCAGATCCGCCCCCTTGACGGCACCCAGGTGAAAGTCGACGTCTGCGGCATTCAGCGCGATGACGATTTACCGTTCATCGGCAGTGGCGATGACAACCGTCAgacgggcgggg                                                                         >  2:281945/1‑241 (MQ=255)
                                                                                                                                                                   aCGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTgg                                                                         <  1:369906/251‑1 (MQ=255)
                                                                                                                                                                         aaaaGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAgg                                                                   >  1:241929/1‑250 (MQ=255)
                                                                                                                                                                         aaaaGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATCTCCGCCACCAGCGACGAAGGGATGGAACAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAGACGGTCAGCCGCGCGTGGGcaagg                                                                   >  1:285883/1‑247 (MQ=255)
                                                                                                                                                                           aaGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTgg                                                                 <  1:289748/251‑1 (MQ=255)
                                                                                                                                                                                 aTCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTg                                                           >  1:152672/1‑251 (MQ=255)
                                                                                                                                                                                                gAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTCATCTGGCTGGAAACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGAGGCATTCCGCGCGATGTCGATTTAACGTTGGTAGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCACGTGGGTACTTGTACACATTGGCttt                                            >  2:52780/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGTCAATCCAGCAATTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGCTCCGCACCCTTGACGCCAACAAGGCGAAAGTCGTGGTCTGCGGCATTCAGCGCGCTGGCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGAGTGGGCCAGTGGGTACTGGGACACGGTAGCtttt                                           >  2:300204/1‑250 (MQ=255)
                                                                                                                                                                                                             aaaTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGCGACACAGCGATGTGCATAGGCGTTCCAGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGCTGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCGTGAGCGta                                >  2:72802/1‑250 (MQ=255)
                                                                                                                                                                                                             aaaTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGt                                                                                      >  2:168124/1‑196 (MQ=255)
                                                                                                                                                                                                             aaaTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGt                                                                                      <  1:168124/196‑1 (MQ=255)
                                                                                                                                                                                                               aTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAttt                                                                                                                     >  1:97758/1‑163 (MQ=255)
                                                                                                                                                                                                               aTTGAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGAttt                                                                                                                     <  2:97758/163‑1 (MQ=255)
                                                                                                                                                                                                                  gAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGtaattaat                          <  1:30319/251‑1 (MQ=255)
                                                                                                                                                                                                                  gAAATCATCCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGAGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGtaattaat                          >  1:223799/1‑251 (MQ=255)
                                                                                                                                                                                                                        acacTTATTTCCGCCCCCAGCGACGAAGGTATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCTGCAGCTGCGATGAACACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGcc                    <  1:406192/248‑1 (MQ=255)
                                                                                                                                                                                                                         tcCTTATTTCCGCCACCAGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCg                   >  1:50673/1‑251 (MQ=255)
                                                                                                                                                                                                                                      caccaGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGCGCGTAATTAATGAAGCCGAAGCAAGCGacac      >  2:463617/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aGCGACGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCAGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGGGGGTACTGGGCCACGTTGGCTTTGCCATGAGCGTAATTAATGGAGCCGAAGCACGCGACActcac  >  2:379772/1‑249 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
AACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTGACGGCAACCAGGCGAAAGTCGACGTCTGCGGCATTCAGCGCGATGTCGATTTAACGTTAGTCGGCAGCTGCGATGAAAACGGTCAGCCGCGCGTGGGCCAGTGGGTACTGGTACACGTTGGCTTTGCCATGAGCGTAATTAATGAAGCCGAAGCACGCGACACTCTC  >  NC_000913/2851583‑2852067

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: