Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,636,826 T→G L135L (CTA→CTC)  unknown ← D‑mannonate oxidoreductase (EC 1.1.1.57)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,636,8260TG100.0% 129.1 / NA 43L135L (CTA→CTC) unknownD‑mannonate oxidoreductase (EC 1.1.1.57)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/16);  total (27/16)

CACGATTTTGGTTACTGACTAAAAAGCCAGGTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTTAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGACCGCCCGCTTGATTAAT  >  1/1636582‑1637068
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
cACGATTTTGGTTACTGACTAAAAAGCCAGGTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGAtt                                                                                                                                                                                                                                               <  1:596448/250‑1 (MQ=255)
  cGATTTTGGTTACTGACTAAAAAGCCAGGTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTaaa                                                                                                                                                                                                                                            >  1:713045/1‑251 (MQ=255)
    aTTTTGGTTACTGACTAAAAAGCCAGGTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATa                                                                                                                                                                                                                                          >  1:891354/1‑251 (MQ=255)
                              gTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTaaa                                                                                                                                                                                                                >  1:860731/1‑251 (MQ=255)
                              gTGCAATGGCATT‑CAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGcc                                                                                                                                                                                                                             >  1:240722/1‑237 (MQ=255)
                                 cAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATc                                                                                                                                                                                                              <  1:288729/250‑1 (MQ=255)
                                      gCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATc                                                                                                                                                                                                              >  1:169205/1‑245 (MQ=255)
                                           aCAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATc                                                                                                                                                                                                              <  1:618788/240‑1 (MQ=255)
                                                   gACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATc                                                                                                                                                                                                              >  1:763806/1‑232 (MQ=255)
                                                                                    gCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCg                                                                                                                                                          >  1:316890/1‑251 (MQ=255)
                                                                                    gCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCg                                                                                                                                                          >  1:361643/1‑251 (MQ=255)
                                                                                       tGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGt                                                                                                                                                        <  1:466471/250‑1 (MQ=255)
                                                                                       tGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGt                                                                                                                                                        <  1:418419/250‑1 (MQ=255)
                                                                                         aaGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAACTTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGt                                                                                                                                                     >  1:425567/1‑251 (MQ=255)
                                                                                         aaGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGt                                                                                                                                                     >  1:83437/1‑251 (MQ=255)
                                                                                              aCTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCtt                                                                                                                                                >  1:207963/1‑251 (MQ=255)
                                                                                                               cACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTaccacc                                                                                                                                <  1:769253/250‑1 (MQ=255)
                                                                                                                   ggAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCa                                                                                                                            >  1:514727/1‑250 (MQ=255)
                                                                                                                             gCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGAc                                                                                                                                                                                                                         >  1:311287/1‑147 (MQ=255)
                                                                                                                                          gtggtgTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTacc                                                                                                                                   <  1:227721/220‑1 (MQ=255)
                                                                                                                                             gtgtAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGAtt                                                                                                                                                                                                                                               <  1:226488/109‑1 (MQ=255)
                                                                                                                                             gtgtAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCAcc                                                                                                                          >  1:852232/1‑226 (MQ=255)
                                                                                                                                                    ggCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGatta                                                                                                                                        <  1:874765/205‑1 (MQ=255)
                                                                                                                                                                              tatTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAg                                                                                                      >  1:411241/1‑213 (MQ=255)
                                                                                                                                                                              tatTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTCTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTa                                                                >  1:48506/1‑251 (MQ=255)
                                                                                                                                                                                      ccTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGttt                                                                                                                                                                                          <  1:301257/121‑1 (MQ=255)
                                                                                                                                                                                      ccTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGAt                                                                              >  1:261898/1‑229 (MQ=255)
                                                                                                                                                                                                  aTCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGAtt                                                                                                                                                                                                                                               >  1:153416/1‑56 (MQ=255)
                                                                                                                                                                                                  aTCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAacac                                            <  1:312376/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaaaCCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGc                             <  1:89470/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaCCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAg                               <  1:109887/247‑1 (MQ=255)
                                                                                                                                                                                                                   aaCCTGTGTAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCtt                           >  1:746095/1‑251 (MQ=255)
                                                                                                                                                                                                                        gtgtAGGTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGcc                                                                                                                                                                                                                             >  1:7376/1‑52 (MQ=255)
                                                                                                                                                                                                                             ggTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATg                                                                                                             <  1:776569/159‑1 (MQ=255)
                                                                                                                                                                                                                             ggTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAatat                                                            <  1:41781/208‑1 (MQ=255)
                                                                                                                                                                                                                             ggTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAaca                                             >  1:601411/1‑223 (MQ=255)
                                                                                                                                                                                                                             ggTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGAc                  >  1:571418/1‑250 (MQ=255)
                                                                                                                                                                                                                             ggTAATAATGTACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGAc                  >  1:893241/1‑250 (MQ=255)
                                                                                                                                                                                                                                       tACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAatat                                                            >  1:807349/1‑198 (MQ=255)
                                                                                                                                                                                                                                       tACCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGAc                  >  1:548796/1‑240 (MQ=255)
                                                                                                                                                                                                                                        aCCTAAATAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAat                                                              <  1:499518/195‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaTAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAg                                                                                                      >  1:45282/1‑151 (MQ=255)
                                                                                                                                                                                                                                            aaaTAGTTGAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGACCGCCCGCTTGATTAAt  >  1:276667/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
CACGATTTTGGTTACTGACTAAAAAGCCAGGTGCAATGGCATTACAACGAAGACCTACTTTAGAAAAGTAAACGGCAAGCCATTGCGTGAAGTTACTAATTGCAGCTTTAGCACCGGAATAAGCTGCAATTTTTGTTAGTGGTGTAAAGGCATTCATACTTGAAATATTGATGATATTTGCACCTTTTTTACCAATCATATCTTTCGCAAAAACCTGTGTAGGTAATAATGTACCTAAATAGTTTAGATTAAATACAAATTCAATGCCTGACTTATCTAAATCAAAGAACGTTTTTGTTGTTTCATTCAAGTCAAATTGATGAAATTCATTATCGGTAGTCGCTTTTGGATTATTACCACCGGCACCATTAATTAAAATGTCACAAGTGCCGAAATCGGCTTCAATTTGATTTCTGACTTCTTTAATATTTTCAAGTTCTAACACATTTGTTTTATAGGCTTTTGCACGACCGCCCGCTTGATTAAT  >  1/1636582‑1637068

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