Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,350 A→T 7.3% L175Q (CTG→CAG)  gatD ← galactitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,3500AT 7.3% 143.4 / 4.5 56L175Q (CTG→CAG) gatDgalactitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base T (2/2);  total (28/28)
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.68e-01

CTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAC  >  NC_000913/2172209‑2172492
                                                                                                                                             |                                                                                                                                              
ctctGCATTCGCGGCGCGCTCCTTTCGCTACTGTTAGCTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAAGGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcag                                                                                                                                               <  2:250013/143‑1 (MQ=255)
  ctGCATTTGCGGCTCGCTCATCTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagca                                                                                                                                             <  2:350123/143‑1 (MQ=255)
          gcggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                                                                     <  2:928314/143‑4 (MQ=255)
          gcggcgCGCTCATTTCGCGACTGTTAAATGGTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATg                                                                                                                                     <  1:981219/143‑1 (MQ=255)
             gcgcgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagt                                                                                                                                  <  2:468885/143‑7 (MQ=255)
                cgcTCATTTCGCTACTGTTAAACGCTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCg                                                                                                                               <  1:384941/143‑1 (MQ=255)
                  cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGc                                                                                                                             >  2:226732/1‑143 (MQ=255)
                  cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTCGCTCCCAGCGCGACAGCGCACTGCATGGCCAGCAGCCCAACGGCTCCGGc                                                                                                                             >  2:239014/1‑143 (MQ=255)
                   tCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCa                                                                                                                            >  1:501390/1‑143 (MQ=255)
                      tttCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATCGCCAGCAGGCCAATGGTTCCGGCAccc                                                                                                                         >  2:643655/1‑142 (MQ=255)
                          gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGCCACCcataa                                                                                                                     >  2:248886/1‑143 (MQ=255)
                           ccaCTGGTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACACATCGCTCCTAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAACGGATCCGGCACCactaat                                                                                                                    >  2:1263624/3‑143 (MQ=255)
                           cTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATCGCAAGCCGGCCAATGGTTCCGGCACAactaat                                                                                                                    >  2:804497/1‑143 (MQ=255)
                                 ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                              >  1:836416/1‑143 (MQ=255)
                                       gTTTGCATCGCACCGAAAGATTTTGCCAGTGCAAGTTATCTGGTATTAATGTGGATCGCCGTCACACTCTTTGTGGTCATCGCGATAGTGCACTGAATGTCGAGCAGTCTAATGTTTCCGGCAGCAATAAAAATAACGttttt                                                                                                        >  2:153408/1‑143 (MQ=255)
                                        tttGCATGGCACCGAATGATTTTGCCAGTGCCAGTGTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                       <  2:707977/143‑1 (MQ=255)
                                           gCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAAGAACGTTTTTATtc                                                                                                    >  1:1165444/1‑143 (MQ=255)
                                            cATCGCACCGAACGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                   <  1:220376/143‑1 (MQ=255)
                                            cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATgct                                                                                                   >  1:334570/1‑143 (MQ=255)
                                             aTCGCACCGAAAGAGTTTGCCAGTGCCAGTTGTTCTCAACTAATGTCGACCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATcctc                                                                                                  <  2:222957/143‑1 (MQ=255)
                                              tCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTAGATCGCCGTCATGTTCTGTGGTCACACACTTACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                 >  1:295471/1‑143 (MQ=255)
                                               cGCACCGATAGAGTGTGCCCGTGCCAGTCTTTCTGCACTAATGTCGATCGCCGTCCCACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                                                                <  2:613430/143‑1 (MQ=255)
                                                gCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCCCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGCTCCGGCACCACTAATAATCACGTTTTTATACCcaca                                                                                               >  2:84844/1‑143 (MQ=255)
                                                 cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  1:818806/143‑1 (MQ=255)
                                                 cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  1:431544/143‑1 (MQ=255)
                                                 cACCGAAAGATTCTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  2:475556/143‑1 (MQ=255)
                                                        aGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcg                                                                                       <  1:766865/143‑1 (MQ=255)
                                                           ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTa                                                                                    >  1:447552/1‑143 (MQ=255)
                                                                  agtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaa                                                                             <  1:261164/143‑1 (MQ=255)
                                                                   gtgccaggTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaaa                                                                            <  1:338736/143‑1 (MQ=255)
                                                                        ggtggtgcTGAACTAGTGGTGTGCGCAGTCACACGCGTTGCTGCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCGATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                       <  2:507297/136‑1 (MQ=255)
                                                                             tgcTGAACTAATGCCGATCGCCGTCACACTCTTTGCCCCCAGCGGGTCAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAc                                                                  <  2:836416/141‑1 (MQ=255)
                                                                                tGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCaa                                                               >  2:279219/1‑143 (MQ=255)
                                                                                 gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCAAAACCTTGCGCTAAATGAAACGCATGAAGACCCAc                                                              >  1:494903/1‑143 (MQ=255)
                                                                                        tGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATCCCCACAACCTTCGCCTAAATGAAACGCAGGAAGCCCAACGGaaacg                                                       >  2:869461/1‑141 (MQ=255)
                                                                                         gTCGATCGTCGTCGGACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATCGTTCCGGCACCAATAATAATAACGTTTTTAGTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg                                                      <  2:917684/143‑1 (MQ=255)
                                                                                           cGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGATAATAGGc                                                    >  2:396508/1‑143 (MQ=255)
                                                                                            gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCCCCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTCACTGAACACCTGCAGACCAACGGTAATCCGccc                                                  >  1:146134/1‑141 (MQ=255)
                                                                                                      acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                         <  1:696175/143‑1 (MQ=255)
                                                                                                      acacTCCTTGCCCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                         <  2:800856/143‑1 (MQ=255)
                                                                                                       cacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACGAACGGTAATCGGCTCAATAAAAGcc                                        >  1:1129467/1‑143 (MQ=255)
                                                                                                        cctctTTGCTCCCAGCGCGACACCACACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                       <  2:501390/142‑1 (MQ=255)
                                                                                                        acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                       >  2:749738/1‑143 (MQ=255)
                                                                                                         ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGcccc                                      >  1:16512/1‑143 (MQ=255)
                                                                                                            tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCa                                                                     <  2:675617/109‑1 (MQ=255)
                                                                                                            tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCa                                                                     >  1:675617/1‑109 (MQ=255)
                                                                                                              tGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATACAAGCCCCATCCt                                 >  1:1182296/1‑143 (MQ=255)
                                                                                                                   ccggcgcGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                            <  2:295471/140‑1 (MQ=255)
                                                                                                                      gcgcGACAGCGCACTGAATGGCCAGCAGGCCAGTGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGc                         <  2:885582/143‑1 (MQ=255)
                                                                                                                              gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGt                 <  1:1143989/143‑1 (MQ=255)
                                                                                                                               ccaaCTGGATGCCCAGCATGTCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCATGCGCTAAATCAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCGCAATAGGCATATCCGTg                <  2:1109004/140‑1 (MQ=255)
                                                                                                                                gcACTGGATGGCCCGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg               <  2:164562/143‑1 (MQ=255)
                                                                                                                                gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg               <  1:298486/143‑1 (MQ=255)
                                                                                                                                 cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATATTAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGGCACATCCTCAATAGGCATATCCGTggg              >  1:105894/1‑143 (MQ=255)
                                                                                                                                     cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATGCGGGGGTaga          >  1:792492/4‑143 (MQ=255)
                                                                                                                                      attGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag         >  2:725312/3‑143 (MQ=255)
                                                                                                                                        tGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCa       <  2:930813/143‑1 (MQ=255)
                                                                                                                                           ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCACCCTCAATAGGCATATCCGCGCGTAGAGCAAAg    >  2:937599/1‑143 (MQ=255)
                                                                                                                                            cagcagGCCACTGGGTCCCGCACCAATAATAATAACGCTTTTACTCTCACAGCCTTGCGCTAAATGAAAAGCATGCAGCCCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa   <  2:937236/143‑1 (MQ=255)
                                                                                                                                            cagcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATACGCATATCCCTCGGCAGAGCAAAGa   >  2:256562/1‑143 (MQ=255)
                                                                                                                                             agcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGCGGGTAGAGCAAAGAc  >  1:150358/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
CTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAC  >  NC_000913/2172209‑2172492

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