Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 816,907 C→T intergenic (+267/‑222) unknown → / → unknown NADP‑specific glutamate dehydrogenase (EC 1.4.1.4)/Catalase KatE (EC 1.11.1.6)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1816,9070CT100.0% 159.3 / NA 46intergenic (+267/‑222)unknown/unknownNADP‑specific glutamate dehydrogenase (EC 1.4.1.4)/Catalase KatE (EC 1.11.1.6)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/18);  total (28/18)

CGCACTTTGATCTGCCCTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCCGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCAAAAATTAAGGAAACGATTATGT  >  1/816671‑817132
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
cGCACTTTGATCTGCCCTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     <  1:320887/251‑1 (MQ=255)
         aTCTGCCCTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACTTTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATc                                                                                                                                                                                                             <  1:725100/250‑1 (MQ=255)
              cccTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTg                                                                                                                                                                                                       <  1:12829/251‑1 (MQ=255)
              cccTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTg                                                                                                                                                                                                       <  1:463302/251‑1 (MQ=255)
                           gACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCAttttc                                                                                                                                                                                            <  1:149112/249‑1 (MQ=255)
                                 gccgccAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     >  1:346619/1‑218 (MQ=255)
                                  ccgccAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCAt                                                                                                                                                                                   >  1:276657/1‑251 (MQ=255)
                                                       ggTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCa                                                                                                                                                                                                 <  1:318319/216‑1 (MQ=255)
                                                       ggTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTaa                                                                                                                                                                 >  1:194541/1‑248 (MQ=255)
                                                       ggTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAAt                                                                                                                                                               >  1:406099/1‑250 (MQ=255)
                                                                 ggAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     >  1:19591/1‑186 (MQ=255)
                                                                      aCATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACt                                                                                                                                                    >  1:315226/1‑246 (MQ=255)
                                                                          tattTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGt                                                                                                                                           <  1:736452/251‑1 (MQ=255)
                                                                                    gAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     >  1:6114/1‑167 (MQ=255)
                                                                                         tttttATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGt                                                                                                                             >  1:726224/1‑250 (MQ=255)
                                                                                               ttCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAatat                                                                                                                      <  1:713256/251‑1 (MQ=255)
                                                                                                       ggTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAg                                                                                                                                            <  1:669141/221‑1 (MQ=255)
                                                                                                        gTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     <  1:238359/147‑1 (MQ=255)
                                                                                                                   tgatgaACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTgag                                                                                                  >  1:540770/1‑251 (MQ=255)
                                                                                                                                cGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTg                                                                                      <  1:245637/250‑1 (MQ=255)
                                                                                                                                     gATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATc                                                                                                                                                              >  1:742745/1‑173 (MQ=255)
                                                                                                                                        gAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     >  1:300286/1‑115 (MQ=255)
                                                                                                                                        gAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAggg                                                                                                                                                                                                                     <  1:248940/115‑1 (MQ=255)
                                                                                                                                          atatatTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAg                                                                                                                                            >  1:94732/1‑186 (MQ=255)
                                                                                                                                                   ttattaTTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGa                                                                                                                                                                                                    >  1:220994/1‑121 (MQ=255)
                                                                                                                                                   ttattaTTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATAATTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCt                                                                                                                                                                        <  1:52140/149‑1 (MQ=255)
                                                                                                                                                      ttattGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTc                                                                                                                                                                         <  1:464573/145‑1 (MQ=255)
                                                                                                                                                      ttattGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATaa                                                                            >  1:609570/1‑238 (MQ=255)
                                                                                                                                                      ttattGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTc                                                                        >  1:243857/1‑242 (MQ=255)
                                                                                                                                                                gCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAAc                                                                                                                                                     >  1:784817/1‑155 (MQ=255)
                                                                                                                                                                              aaGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTc                                                                                                                                                         >  1:214429/1‑137 (MQ=255)
                                                                                                                                                                                  tAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTAc                                   <  1:575162/251‑1 (MQ=255)
                                                                                                                                                                                              gCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCaa                       >  1:110976/1‑251 (MQ=255)
                                                                                                                                                                                                   gtgcagtgcaGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATc                                                                                                                                                              >  1:42245/1‑111 (MQ=255)
                                                                                                                                                                                                   gtgcagtgcaGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACt                                                                                                                                                    >  1:58773/1‑121 (MQ=255)
                                                                                                                                                                                                   gtgcagtgcaGAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCAAAAAtt                  >  1:287428/1‑251 (MQ=255)
                                                                                                                                                                                                     gcagtgcagAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAAt                                                                 <  1:421768/202‑1 (MQ=255)
                                                                                                                                                                                                        gtgcagAGCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTc                         >  1:111448/1‑239 (MQ=255)
                                                                                                                                                                                                             gagCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACt                                                                                                                                                    >  1:629756/1‑111 (MQ=255)
                                                                                                                                                                                                               gCTCAATATTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTa                                                                                                                                         >  1:723504/1‑120 (MQ=255)
                                                                                                                                                                                                                    atatTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATattatt                                                                                                                  >  1:871978/1‑138 (MQ=255)
                                                                                                                                                                                                                    atatTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGACAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCAAAAATTAAGGAAACGATTATGt  >  1:194101/1‑250 (MQ=255)
                                                                                                                                                                                                                      atTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAAc                                                                                               >  1:793150/1‑155 (MQ=255)
                                                                                                                                                                                                                      atTTAACTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATAtttt                                       <  1:342551/211‑1 (MQ=255)
                                                                                                                                                                                                                          aaCTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAAc                                                                                               <  1:525501/151‑1 (MQ=255)
                                                                                                                                                                                                                           aCTTAGGTTTAGCTATCTGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCAAAAATTAAgg              >  1:529649/1‑231 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
CGCACTTTGATCTGCCCTCCCAAAGAGGACGCCGCCGCCAATCAATTAAAATGCAGGTTTTTATAGGAAAACATTATTTAATCGGAATATTTTTATTCCCGCCGGTTTAACTACATGATGAACCCGGGCGGACGATGAATATATTCGTTATTATTGAGGAGCGTAATTAAGCTTAAGTTAAACAATTTGAGCCCTGTGCAGTGCAGAGCTCAATATTTAACTTAGGTTTAGCTATCCGCCCGGTTTTAGGGTGCAGATCAAAGTGCGATCATTTTCTTTTCTCATTTCTTTTTTCTTGATTAAATCAATTCTAACTATTAATAGTTATTTTTTAATCGTCAAATATTATTAATTTAATAGATTGAGAACTTTTTGCTGGTATTTATAATTTCTGTCAATCGATCATCCATAAAAATGGATATTTTTTACTCTACTACTCAAAAATTAAGGAAACGATTATGT  >  1/816671‑817132

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