Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 554,812 A→G 6.8% T68A (ACC→GCC)  cysS → cysteinyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913554,8120AG 6.8% 167.9 / 3.4 59T68A (ACC→GCC) cysScysteinyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACC  >  NC_000913/554665‑554959
                                                                                                                                                   |                                                                                                                                                   
cACGCCGGGGAAGTCGGCTTGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCGTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTAcc                                                                                                                                                   >  1:78721/1‑150 (MQ=255)
cACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTAcc                                                                                                                                                   >  1:125466/1‑150 (MQ=255)
    ccGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTAACGATCTCTGGCATATCGGGCACGGGCGTACCTTTGTTGCGTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGata                                                                                                                                               <  2:807440/150‑1 (MQ=255)
      ggtgAAGTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCCGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTATCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATc                                                                                                                                             <  2:554126/147‑1 (MQ=255)
          aaGTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATcgacg                                                                                                                                         <  1:859066/150‑1 (MQ=255)
            gTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCgacgac                                                                                                                                       >  1:571119/1‑150 (MQ=255)
              cGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACaa                                                                                                                                     >  2:796212/1‑150 (MQ=255)
               ggCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACaaa                                                                                                                                    >  1:336117/1‑150 (MQ=255)
                  aTGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTAACGATATCGACGACAAAatc                                                                                                                                 >  1:327205/1‑150 (MQ=255)
                      aCGTGTGTGGAATGAACGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGTTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAatcatga                                                                                                                             <  1:716553/150‑3 (MQ=255)
                            gtgGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAAcgcg                                                                                                                       >  1:35553/1‑150 (MQ=255)
                             gggAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACgcgc                                                                                                                      <  1:993427/149‑1 (MQ=255)
                              ggAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCGTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGcc                                                                                                                     <  2:78721/150‑1 (MQ=255)
                                       gTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAAAGCGCCAATGAAAAt                                                                                                            >  2:347064/1‑150 (MQ=255)
                                             gATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATGACCGATATCGACGACAAAATCATCAAACGAGCCAATGAAAATGGCGaa                                                                                                      >  1:776304/1‑150 (MQ=255)
                                              aTCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGaaa                                                                                                     >  1:334338/1‑150 (MQ=255)
                                                  ctGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCtt                                                                                                 <  2:1105110/150‑1 (MQ=255)
                                                    gTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTtg                                                                                               >  1:967637/1‑150 (MQ=255)
                                                     tCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTtgt                                                                                              <  1:281624/150‑1 (MQ=255)
                                                     tCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTtgt                                                                                              <  1:626410/150‑1 (MQ=255)
                                                     tCATATCGGTCACGGGCGTACCATTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCTGCTATAAACTGAAGTATGAGCGCAACATTACCGATATCGACGACAAACTCATCAAACGCGCCAATGAAAATGGCGAAAGCTTtgt                                                                                              <  1:164132/150‑1 (MQ=255)
                                                      cATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACCGAAGTATGTGCGCAACATTACCGAAATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGCAAGCTTtgtg                                                                                             >  2:609474/1‑150 (MQ=255)
                                                       atatCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTgg                                                                                            <  2:1344071/150‑1 (MQ=255)
                                                       atatCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTgg                                                                                            <  1:194990/150‑1 (MQ=255)
                                                       atatCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTgg                                                                                            <  2:1133362/150‑1 (MQ=255)
                                                        tatCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCTATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGc                                                                                           <  2:1208008/150‑1 (MQ=255)
                                                           cGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATTAACTGAAGTATGTGCGCAACATTACCGATATCGACGCCAAAATCATCAAACGCGCCAATGAAAATGGCGATAGCTTTGTGGCGAt                                                                                        >  1:1163134/1‑150 (MQ=255)
                                                            ggTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAGCATTGCCGATATCGACGACTAAAGCATCAAACGCGCCAAGGAAAAGGTCGAAAGCTTTGTGGCGAtg                                                                                       >  2:103223/1‑150 (MQ=255)
                                                            ggTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTGCCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAtg                                                                                       >  2:1275382/1‑150 (MQ=255)
                                                             gTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCTGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACTCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAtgg                                                                                      <  1:819844/150‑1 (MQ=255)
                                                             gTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGTCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAtgg                                                                                      <  2:984390/150‑1 (MQ=255)
                                                             gTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAtgg                                                                                      >  1:1288081/1‑150 (MQ=255)
                                                                aCGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTACGCAACATTACCGATATCGACGACAAAATCATCAACCGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAtggtgg                                                                                   <  1:801257/150‑1 (MQ=255)
                                                                       tACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCAt                                                                            >  2:711189/1‑150 (MQ=255)
                                                                        aCCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATg                                                                           <  2:183627/150‑1 (MQ=255)
                                                                              gttgCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGcc                                                                     <  2:728799/150‑1 (MQ=255)
                                                                              gttgCAGTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGcc                                                                     <  1:202121/150‑1 (MQ=255)
                                                                                     ttGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGc                                                              <  1:419452/150‑1 (MQ=255)
                                                                                       gACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGcac                                                            <  1:1265599/150‑1 (MQ=255)
                                                                                        aCGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAATTGcaca                                                           >  2:996374/1‑150 (MQ=255)
                                                                                                 cgcgcTATCTGCGTTTCCTCGTCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTg                                                  <  2:742339/150‑1 (MQ=255)
                                                                                                 cgcgcTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGACTTTg                                                  >  2:1067751/1‑150 (MQ=255)
                                                                                                      tATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCACCATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCt                                             >  1:893216/1‑150 (MQ=255)
                                                                                                        tCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCttt                                           >  1:275421/1‑150 (MQ=255)
                                                                                                            cGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAACCGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGaaa                                       >  2:1281472/1‑149 (MQ=255)
                                                                                                                tCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCAGCAAACGCGCCAATGAAAATGGCGATAGCGTTGTGGAGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATGTTTATGCTTTGAACAtag                                   >  2:262639/1‑148 (MQ=255)
                                                                                                                 ccTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCt                                  <  2:571119/150‑1 (MQ=255)
                                                                                                                    cGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGAGGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTgcg                               <  1:837448/150‑1 (MQ=255)
                                                                                                                           aaaCTGTAGTATGTGCGCAACATTACCGATATCGACGTCAAAATCATCAAACGCGCCAATCATAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCTCCGAAATGCACAAAGAGTTTGATGCTTTGAACATTCTGCGCCCGGat                        <  2:275572/150‑1 (MQ=255)
                                                                                                                              cTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAACCGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATg                     <  1:1183237/150‑1 (MQ=255)
                                                                                                                                  aGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGc                 >  2:28638/1‑150 (MQ=255)
                                                                                                                                   gTATGTGCGCAACATTACCGATATCGACGACAAAATCCTCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGcc                <  2:846520/150‑1 (MQ=255)
                                                                                                                                         gcgcAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCcgcgcgc          >  2:417666/1‑150 (MQ=255)
                                                                                                                                           gcAACATTGCCGATATCGACGACAAAATCATCTCACGCGCCACTGAAATAGGCGAAAGCTTTGAGGCGATGGAGGATCGCAGGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAAGATTCTGCGCCCGGATATGGAGCCGCGCGCGa        <  1:103223/150‑1 (MQ=255)
                                                                                                                                           gcAACATTGCCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGa        <  1:1275382/150‑1 (MQ=255)
                                                                                                                                             aaCATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACg      >  2:604766/1‑150 (MQ=255)
                                                                                                                                                aTTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCAc   <  2:1231306/150‑1 (MQ=255)
                                                                                                                                                aTTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCAc   >  1:775204/1‑150 (MQ=255)
                                                                                                                                                 ttACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCAcc  <  2:1288081/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTTACGATCTCTGTCATATCGGTCACGGGCGTACCTTTGTTGCTTTTGACGTGGTTGCGCGCTATCTGCGTTTCCTCGGCTATAAACTGAAGTATGTGCGCAACATTACCGATATCGACGACAAAATCATCAAACGCGCCAATGAAAATGGCGAAAGCTTTGTGGCGATGGTGGATCGCATGATCGCCGAAATGCACAAAGATTTTGATGCTTTGAACATTCTGCGCCCGGATATGGAGCCGCGCGCGACGCACC  >  NC_000913/554665‑554959

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