Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 139,326 C→T G634S (GGC→AGC)  gcd ← glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913139,3260CT97.9% 180.6 / NA 47G634S (GGC→AGC) gcdglucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/24);  total (22/25)
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

CGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCA  >  NC_000913/139153‑139529
                                                                                                                                                                              |                                                                                                                                                                                                            
cggcgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                                                                                                                                                                       >  1:150446/1‑212 (MQ=255)
cggcgTACCAATACGTTTCTTCCACAACACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                                                                                                                                                                       <  2:150446/212‑1 (MQ=255)
    gTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAg                                                                                                                                        <  3:153958/239‑1 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAt                                                                                                                                                             <  2:153250/217‑1 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAt                                                                                                                                                             >  1:153250/1‑217 (MQ=255)
     tACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGt                                                                                                                                       >  3:341431/1‑239 (MQ=255)
           tACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 <  3:141306/239‑1 (MQ=255)
                ttcttcCACACGCATTCATTGGTTTTCAGATCCAGCGCCGAGATATCACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACCCCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGt                                                                                                                             <  4:341431/238‑1 (MQ=255)
                   ttcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  4:117017/182‑1 (MQ=255)
                   ttcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  3:117017/1‑182 (MQ=255)
                                cATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCa                                                                                                            <  1:34381/239‑1 (MQ=255)
                                                             taACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACAACGCGGGATCAGTTTCGCAACAAACGGCAGTGCCATTGGGTTGGGAATCGTCACGTCACGATTTGGATCAACGGAAATCTCGtcc                                                                                 >  4:155009/1‑237 (MQ=255)
                                                                  ccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  3:230515/135‑1 (MQ=255)
                                                                  ccATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  4:230515/1‑135 (MQ=255)
                                                                    aTGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAAc                                                                        >  1:95466/1‑239 (MQ=255)
                                                                          gctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCGGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATTAACGGAAATCCCGCCCCAATCGAAAATtccc                                                                  >  4:317991/1‑239 (MQ=255)
                                                                          gctgTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGCAATCCCGCCCCATTCGAACATcccc                                                                  >  3:292323/1‑239 (MQ=255)
                                                                               ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 <  2:228688/171‑1 (MQ=255)
                                                                               ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAtt                                                                                                                                 >  1:228688/1‑171 (MQ=255)
                                                                                          aCCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAg                                                   >  3:162451/1‑238 (MQ=255)
                                                                                          aCCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGACCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTTACGATTTGGATCAACGGACATACAGCCGCATTCGTACAATCCCAGGTTACCCGGGAAgc                                                  >  2:67225/1‑238 (MQ=255)
                                                                                          aCCAAATGGTAGGAAGACGGGGTTGAGCGTAGACACAAT‑CGGTACACCGTACTGTGGATAAATCCCGGTCTCCGTACCCGTCCTTTTGG‑CCTCTTTCGGCTGCTCCATCGGGTTTCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGa                                                  <  4:292323/239‑1 (MQ=255)
                                                                                           ccAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCGGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGa                                                  <  3:317991/238‑1 (MQ=255)
                                                                                               aTGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTTAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGc                                             <  3:155009/239‑1 (MQ=255)
                                                                                                      gagGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGACCATCCCCAGGTTAACCGGGAAGACCAGCggaa                                         >  4:476304/1‑233 (MQ=255)
                                                                                                          aaCGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTAccc                                       <  3:239800/234‑1 (MQ=255)
                                                                                                                 tGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATgg                            <  3:476304/238‑1 (MQ=255)
                                                                                                                     catgtacgccgtACGGTAAGCGGGACTGGGGCTGAATGGGGGATCCGGTA‑CCGGGACTATGGACATCGTTCGGCTGCTCCAGCGGTGTGCCAGGACCACGCGGGATCAGTTGCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTGCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACAGCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGAggg                            <  4:325497/224‑1 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGgct                                                                                                                                                                                            <  4:355815/70‑1 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGgct                                                                                                                                                                                            >  3:355815/1‑70 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAc                                                                                                                                                                    >  3:290504/1‑94 (MQ=255)
                                                                                                                          acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAc                                                                                                                                                                    <  4:290504/94‑1 (MQ=255)
                                                                                                                                        acCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         >  3:185673/1‑107 (MQ=255)
                                                                                                                                        acCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         <  4:185673/107‑1 (MQ=255)
                                                                                                                                           gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  3:283090/239‑1 (MQ=255)
                                                                                                                                            tACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         <  1:224139/103‑1 (MQ=255)
                                                                                                                                            tACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCa                                                                                                                                         >  2:224139/1‑103 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  <  3:323377/52‑1 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCgg                                                                                                                                                                                  >  4:323377/1‑52 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  >  1:5852/1‑228 (MQ=255)
                                                                                                                                                      gAATGCCGGATTCCGTACCCGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  2:5852/228‑1 (MQ=255)
                                                                                                                                                                        ccGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                                                          >  3:318874/1‑138 (MQ=255)
                                                                                                                                                                        ccGTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                                                          <  4:318874/138‑1 (MQ=255)
                                                                                                                                                                          gTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  <  3:106574/208‑1 (MQ=255)
                                                                                                                                                                          gTGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCa  >  4:106574/1‑208 (MQ=255)
                                                                                                                                                                           tGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATc                                                                     >  4:328242/1‑140 (MQ=255)
                                                                                                                                                                           tGCTTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATc                                                                     <  3:328242/140‑1 (MQ=255)
                                                                                                                                                                              |                                                                                                                                                                                                            
CGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGT‑GACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGG‑CATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATGGCGGGGTGAAAATGCCTTCATAGCGCA  >  NC_000913/139153‑139529

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