Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,714,436 2 bp→GC coding (261‑262/660 nt) yiaD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,4360CG99.0% 255.4 / ‑4.9 100N87K (AAC→AAGyiaDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (52/47);  minor base A (1/0);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
*NC_0009133,714,4370GC98.0% 244.4 / NA 100V88L (GTG→CTG) yiaDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (50/47);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GCTGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAACGTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATAT  >  NC_000913/3714364‑3714504
                                                                        ||                                                                   
gcTGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGGAGGGAGCGGAGGTg                                                                   >  1:1621173/1‑76 (MQ=255)
  tGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGGTATTACATGGATGTGGGGGAAGCGAAGGTggg                                                                 >  1:279891/1‑74 (MQ=255)
   gATTGGCGCAGCAGCAGGCGCGGCGCTGGGCGGCGGCGTTGGGTATTACATGGCTGGCCAGGAAGCGAAGCTgcgc                                                                >  2:1632178/1‑76 (MQ=255)
   gATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTgcgc                                                                >  2:2483328/1‑76 (MQ=255)
    aTTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGcgcg                                                               >  1:1708381/1‑76 (MQ=255)
     ttGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGa                                                              <  1:1148667/76‑1 (MQ=255)
     ttGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGa                                                              <  1:2642542/76‑1 (MQ=255)
       gtcgcAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGCAGCTGCGCGACa                                                            <  1:3872364/74‑1 (MQ=255)
        gcgcAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCCGGAAGCGCAGCTGCGCGACaa                                                           >  2:1976061/1‑76 (MQ=255)
        gcgcAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACaa                                                           <  1:3231701/76‑1 (MQ=255)
         cgcAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCCCGACaaa                                                          >  2:3355854/1‑76 (MQ=255)
           cannngcagGCGCAGCTCTGNGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGANGCTGCGCGACAAAAt                                                        <  1:4597361/76‑1 (MQ=255)
           cagcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAAt                                                        <  1:3655119/76‑1 (MQ=255)
           cagcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAAt                                                        <  1:3227243/76‑1 (MQ=255)
           cagcagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAAt                                                        <  1:802827/76‑1 (MQ=255)
              cagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcg                                                     >  1:971693/1‑76 (MQ=255)
              cagcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcg                                                     <  1:4718943/76‑1 (MQ=255)
               agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                    <  1:518284/76‑1 (MQ=255)
               agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                    >  1:504773/1‑76 (MQ=255)
               agcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATgcgc                                                    >  1:1045751/1‑76 (MQ=255)
                gcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGcgcg                                                   >  1:2146939/1‑76 (MQ=255)
                gcagGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCCGCGCGACAAAATGcgcg                                                   >  2:4127507/1‑76 (MQ=255)
                 cagGCGCGGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCgg                                                  <  2:1318732/76‑1 (MQ=255)
                 cagGCGCATCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCgg                                                  <  2:1362809/76‑1 (MQ=255)
                  agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                 >  1:602042/1‑76 (MQ=255)
                  agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGc                                                 >  1:2111598/1‑76 (MQ=255)
                   ggCGCAGCTCGGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCa                                                <  2:4310551/76‑1 (MQ=255)
                    gcgcCGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCAc                                               <  2:3597544/76‑1 (MQ=255)
                     cgcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAACGCGCGGCAcc                                              >  2:2994882/1‑75 (MQ=255)
                      gcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTg                                             <  1:1598701/76‑1 (MQ=255)
                       cAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTgg                                            <  1:2314153/76‑1 (MQ=255)
                        aGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgt                                           >  1:2487342/1‑76 (MQ=255)
                         gCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACCAAATGCGCGGCACTCgtg                                          >  2:136013/1‑76 (MQ=255)
                          ctctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGGGCGGGCCTGgggt                                         >  2:21068/1‑76 (MQ=255)
                          ctctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGgtgt                                         >  1:3878886/1‑76 (MQ=255)
                           tctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGtt                                        <  1:2786053/76‑1 (MQ=255)
                            ctGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTa                                       <  2:2649808/76‑1 (MQ=255)
                               ggcggcggcGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCg                                    >  2:511980/1‑76 (MQ=255)
                               ggcggcggcGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCg                                    >  1:1811643/1‑76 (MQ=255)
                               ggcggcggcGTTGGGTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAAGGCGCGGCAGTGGTGTTAGCg                                    >  1:4677147/1‑76 (MQ=255)
                               ggcggcggcGTTGATTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCg                                    <  1:2169786/76‑1 (MQ=255)
                                 gggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTa                                  <  1:276221/75‑1 (MQ=255)
                                 cggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCgaa                                  <  1:3474679/76‑3 (MQ=255)
                                 cggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTAGTGTTAGCGTa                                  >  1:3440157/1‑76 (MQ=255)
                                  ggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTaa                                 >  1:2919628/1‑76 (MQ=255)
                                   gcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAc                                >  2:2858457/1‑76 (MQ=255)
                                    cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                               >  1:2407821/1‑76 (MQ=255)
                                    cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAcc                               <  1:2468110/76‑1 (MQ=255)
                                     ggggTTGGGTATTACATGGAGGGGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                              <  2:4750353/73‑1 (MQ=255)
                                     ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                              <  1:2330262/76‑1 (MQ=255)
                                     ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                              >  2:1793193/1‑76 (MQ=255)
                                     ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                              <  1:160993/76‑1 (MQ=255)
                                     ggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAAccc                              >  1:2037019/1‑76 (MQ=255)
                                      gcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                             >  2:2801900/1‑76 (MQ=255)
                                      gcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCg                             >  2:3324205/1‑76 (MQ=255)
                                       cgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGc                            >  2:1775855/1‑76 (MQ=255)
                                        gTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCa                           >  1:1077356/1‑76 (MQ=255)
                                        gTTGGTTATTACATGGATGTGCAGGAAGCGAAACTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCa                           >  1:3106942/1‑76 (MQ=255)
                                        gTTGGTTAGTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCa                           <  2:4844600/76‑1 (MQ=255)
                                         ttGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAg                          >  1:3329981/1‑76 (MQ=255)
                                          tGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGc                         <  1:21068/76‑1 (MQ=255)
                                           ggTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCg                        >  1:3445747/1‑76 (MQ=255)
                                            gTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCgg                       <  2:2002190/76‑1 (MQ=255)
                                             ttattaCATGGATGGGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCggg                      <  2:447471/76‑1 (MQ=255)
                                              tattaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCgggg                     >  2:2481056/1‑76 (MQ=255)
                                               attaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGa                    <  1:3137199/76‑1 (MQ=255)
                                               attaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGa                    >  2:1486695/1‑76 (MQ=255)
                                               attaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGa                    <  1:3525348/76‑1 (MQ=255)
                                               attaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGa                    <  2:525288/76‑1 (MQ=255)
                                                ttaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGAt                   >  1:36423/1‑76 (MQ=255)
                                                 taCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATa                  >  1:4408050/1‑76 (MQ=255)
                                                   cATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAGCCCGCAGCGGGGATAAc                <  1:4730826/76‑1 (MQ=255)
                                                   cATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAAc                >  1:81219/1‑76 (MQ=255)
                                                   cATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAAc                <  2:2884903/76‑1 (MQ=255)
                                                   cATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAAc                >  1:1642573/1‑76 (MQ=255)
                                                    aTGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACa               <  2:1685713/76‑1 (MQ=255)
                                                    aTGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACa               >  2:2301412/1‑76 (MQ=255)
                                                    aTGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGAGAACa               >  2:597281/1‑76 (MQ=255)
                                                     tGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACAt              <  1:4397008/76‑1 (MQ=255)
                                                     tGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAAAAt              >  1:4801940/1‑76 (MQ=255)
                                                        aTGTGCAGGCAGCGCAGCTGCGCTCGAAAATGCGCGGCATTGGTGTTAGCGTCACCCGCAGCGGGGATAACATTAt           <  1:315053/76‑1 (MQ=255)
                                                         tgtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATc          <  1:4037567/76‑1 (MQ=255)
                                                          gtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATcc         <  2:2416090/76‑1 (MQ=255)
                                                          gtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATcc         >  1:1031914/1‑76 (MQ=255)
                                                           tgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCt        <  1:63467/76‑1 (MQ=255)
                                                            gCAGGAAGCGAAGCTGCGCGCCAAAATGCGCGGCACTGGTGTTAGCGTAGCCCGCAGCGGGGATAACATTATGCTc       <  2:374014/76‑1 (MQ=255)
                                                            gCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTc       <  1:2811949/76‑1 (MQ=255)
                                                            gCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGAACCCCCCAGCGGGGAAACAATTACCCTc       >  2:3437443/1‑76 (MQ=255)
                                                             cAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCa      >  1:3584960/1‑76 (MQ=255)
                                                              aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCCGCGGGGATAACATTATCCTCaa     <  1:665093/76‑1 (MQ=255)
                                                              aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCaa     >  2:2951490/1‑76 (MQ=255)
                                                              aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCaa     >  1:314295/1‑76 (MQ=255)
                                                              aGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGCGTTAGCGTAACCCGCAGCGGGGACAACATTATCCTCaa     >  1:4206760/1‑76 (MQ=255)
                                                               ggAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAat    >  2:2894780/1‑76 (MQ=255)
                                                               ggAAGCGAAGATGCGCGACAAAAAGCGCGGCACTGGTGATAGCGTAACCCGCAGCAGCGATAACATTATCGTcgat    >  1:3466382/1‑73 (MQ=255)
                                                                gaagcgaagcTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAata   <  1:584739/76‑1 (MQ=255)
                                                                gaagcgaagcTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAata   >  1:884985/1‑76 (MQ=255)
                                                                 aaGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAatat  <  1:3265884/76‑1 (MQ=255)
                                                                 aaGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAatat  <  1:3613413/76‑1 (MQ=255)
                                                                 aaGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAatat  <  2:1862621/76‑1 (MQ=255)
                                                                        ||                                                                   
GCTGATTGGCGCAGCAGCAGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAACGTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATAT  >  NC_000913/3714364‑3714504

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