Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,506 C→G 7.6% intergenic (+23/‑109) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5060CG 7.6% 233.0 / 8.1 66intergenic (+23/‑109)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (3/2);  total (34/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 = 1.00e+00

CCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTT  >  NC_000913/3891362‑3891652
                                                                                                                                                |                                                                                                                                                  
ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                <  4:562588/149‑1 (MQ=255)
ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                >  3:562588/1‑149 (MQ=255)
  aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAAtggctggc                                                                                                                                              <  2:281474/149‑1 (MQ=255)
  aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAAtggctggc                                                                                                                                              <  2:876068/149‑1 (MQ=255)
   aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATggctggct                                                                                                                                             <  4:71280/149‑1 (MQ=255)
   aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATggctggct                                                                                                                                             <  1:986618/149‑1 (MQ=255)
                                 ttATTCCGGCGATTGTCATTCTCTTTGGTGTAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                                                                    <  4:1067680/128‑1 (MQ=255)
                                 ttATTCCGGCGATTGTCATTCTCTTTGGTGTAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAAGCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                                                                    >  3:1067680/1‑128 (MQ=255)
                                                   ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                                                                           >  4:798896/1‑103 (MQ=255)
                                                   ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                                                                           <  3:798896/103‑1 (MQ=39)
                                                        tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                        >  1:964869/1‑149 (MQ=255)
                                                              ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                         <  4:1080661/110‑1 (MQ=255)
                                                              ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                         >  3:1080661/1‑110 (MQ=255)
                                                               taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacca                                                                                 >  3:661704/1‑149 (MQ=255)
                                                                   cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                                                                             >  2:483138/1‑149 (MQ=255)
                                                                               cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                                 <  1:483138/149‑1 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTgt                                                                     <  4:71943/142‑2 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     >  3:71943/1‑142 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     <  4:1121726/142‑1 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     >  3:428470/1‑142 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     <  4:428470/142‑1 (MQ=255)
                                                                                  gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     >  3:1121726/1‑142 (MQ=255)
                                                                                      ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                                                                       <  3:1421974/104‑1 (MQ=255)
                                                                                      ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                                                                       >  4:1421974/1‑104 (MQ=255)
                                                                                          aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                            <  3:201610/79‑1 (MQ=255)
                                                                                          aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                            >  4:201610/1‑79 (MQ=255)
                                                                                            gTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                          <  4:613372/79‑1 (MQ=255)
                                                                                            gTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                          >  3:613372/1‑79 (MQ=255)
                                                                                               tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                >  2:1207863/1‑54 (MQ=255)
                                                                                               tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                <  1:1207863/54‑1 (MQ=255)
                                                                                               tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                             >  3:599643/1‑105 (MQ=255)
                                                                                               tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                             <  4:599643/105‑1 (MQ=255)
                                                                                                ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                                >  1:61611/1‑133 (MQ=255)
                                                                                                ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                                <  2:61611/133‑1 (MQ=255)
                                                                                                  aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCtt                                              <  2:964869/149‑1 (MQ=255)
                                                                                                     gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                                  <  4:1233274/62‑1 (MQ=255)
                                                                                                     gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                                  >  3:1233274/1‑62 (MQ=255)
                                                                                                       gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                                 >  4:460518/1‑61 (MQ=255)
                                                                                                       gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                                 <  3:460518/61‑1 (MQ=255)
                                                                                                         tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTGTGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                       <  1:546775/133‑1 (MQ=255)
                                                                                                         tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                       >  2:546775/1‑133 (MQ=255)
                                                                                                            cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCtt                                                                        >  4:241531/1‑113 (MQ=255)
                                                                                                            cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCtt                                                                        <  3:241531/113‑1 (MQ=255)
                                                                                                             tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                   >  1:122265/1‑149 (MQ=255)
                                                                                                                         aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     >  3:404117/1‑103 (MQ=255)
                                                                                                                         aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     <  4:404117/103‑1 (MQ=255)
                                                                                                                         aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                       <  3:489838/117‑1 (MQ=255)
                                                                                                                         aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                       >  4:489838/1‑117 (MQ=255)
                                                                                                                           tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                                                                                <  1:723063/74‑1 (MQ=255)
                                                                                                                           tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                                                                                >  2:723063/1‑74 (MQ=255)
                                                                                                                              ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                                <  1:1337267/103‑1 (MQ=255)
                                                                                                                              ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                                >  2:1337267/1‑103 (MQ=255)
                                                                                                                              ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                  <  2:122265/149‑1 (MQ=255)
                                                                                                                                  agaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     >  2:150164/1‑94 (MQ=255)
                                                                                                                                  agaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                     <  1:150164/94‑1 (MQ=255)
                                                                                                                                  agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                                                                   <  1:1264136/80‑1 (MQ=255)
                                                                                                                                  agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                                                                   >  2:1264136/1‑80 (MQ=255)
                                                                                                                                   gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt             >  1:59792/1‑149 (MQ=255)
                                                                                                                                       ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGccc         <  4:661704/149‑1 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                        <  2:254428/68‑1 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                        >  1:254428/1‑68 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCAt                                      >  2:809480/1‑118 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCAt                                      <  1:809480/118‑1 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTccc                                 >  1:235942/1‑123 (MQ=255)
                                                                                                                                         agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTccc                                 <  2:235942/123‑1 (MQ=255)
                                                                                                                                             aTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                                                                      <  2:724495/66‑1 (MQ=255)
                                                                                                                                             aTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                                                                      >  1:724495/1‑66 (MQ=255)
                                                                                                                                              ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTa                                            <  3:1061250/107‑1 (MQ=255)
                                                                                                                                              ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTa                                            >  4:1061250/1‑107 (MQ=255)
                                                                                                                                              ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtttt  >  4:535970/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTT  >  NC_000913/3891362‑3891652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: