Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC98.3% 178.2 / NA 60I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (36/23);  total (37/23)
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

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925
                                                                                                                                            |                                                                                                                                          
ggATCTGGTAGAGGGTAAACAGCTGGTATTTTGTGCTGGCGCGGTCGAATTCGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                          <  2:1161653/143‑1 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCCGTCGAATTTGTCTTTGTGCTTGCCGTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGCTGATCTGCATTGCGCCGTCATTGATCAGCAGTCTCTGAGTc                                                                                                                                    >  1:203960/1‑143 (MQ=255)
       gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGACg                                                                                                                                   >  1:195877/1‑143 (MQ=255)
          gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAGTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCct                                                                                                                                >  1:128770/1‑143 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAGTTTTTTCGGTGACGAGGTGATCTGGATTTGAGCGTGGTTGATATGCAGGATGTGCGGCGccgc                                                                                                                               >  2:339665/1‑141 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                               >  1:709213/1‑143 (MQ=255)
            gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCccct                                                                                                                              >  2:785635/1‑143 (MQ=255)
                  aCAGCTGATATTATTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACt                                                                                                                        >  1:787123/1‑143 (MQ=255)
                       tgttgttttttGCTTGGGCGGGCGAATTTGTGTTTGTGCTTGCGGTAAAAGTGGTCGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   <  1:770365/138‑1 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   <  1:820162/143‑1 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   >  2:215868/1‑143 (MQ=255)
                        gATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTGCTTTcccc                                                                                                                  >  1:422663/1‑141 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACg                                                                                                                 >  2:1215590/1‑143 (MQ=255)
                              ttttTCTGGCTGTGTCGAATTTGTCTTTGTGTTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGTTGACGATGTGATCTGCATTTCACCGTGTTTGATATGTAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                            <  1:433611/143‑1 (MQ=255)
                               tttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                           <  2:617079/143‑1 (MQ=255)
                                 tGCTGGTGCGGTCGAATTTGTCTTTGTGCTTGTGGTAATAGTTTAAGTTTTTGGATAATTTTTTCGGTGACGAGGTGATCTGCATGGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         <  2:891791/143‑1 (MQ=255)
                                 tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGGAGGCTCTGCGTCGCGTCTTCGACTGTCCACGGTTCGCCg                                                                                                         >  2:231391/1‑143 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTTCTAGCAGTACACACAGTATAGCACTTAGATTTTTTTCGGTGACGACATGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                    <  2:1277142/143‑1 (MQ=255)
                                        gcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTCCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGCCTGTCCACGCTTCGCCGTACAccc                                                                                                  >  1:655825/1‑142 (MQ=255)
                                            tCGATTTTGTCTTTGTGCTTGCGGTATATGTGGTTTTTTCCTGTGATTTTTTTTGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                              <  2:865008/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCTCTGCATTTCAACGTGATTGATATGCAGGCGCTGCGTCGCCTCTTATACTTTCCACGGTTCGACGTCATCCACCATCAACCgaacga                                                                                   >  1:308904/1‑138 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                   >  2:794560/1‑143 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTAATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAATTACCACAGTTCGCACTACTACACCCGCATCACCAGGa                                                                                   >  1:1005266/1‑143 (MQ=255)
                                                           tgctctcGGTAAAAGGGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                               <  2:282581/137‑1 (MQ=255)
                                                           tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                               <  2:316209/143‑1 (MQ=255)
                                                           tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                               <  1:468676/143‑1 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGagaacga                                                                             >  1:1123629/1‑137 (MQ=255)
                                                                    tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                      >  2:911490/1‑143 (MQ=255)
                                                                     aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTacgcc                                                                     >  1:641106/1‑139 (MQ=255)
                                                                        aggggAAGTACCGGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:826632/140‑1 (MQ=255)
                                                                        agcggAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:226409/140‑1 (MQ=255)
                                                                         gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGCACTCCACCACCATCCGCCGGCAGATGTCGTAATCCtgcg                                                                 >  1:542515/1‑141 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGCAGATGTCGTAATCCTCTGc                                                                >  1:280238/1‑143 (MQ=255)
                                                                           ggAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCAGGCTCTGCGTCGCGTGTTTTCGGCCCCACCGCTCGCCGTAATCCACCATCATCCTCAGGAAGATGTCGTAATCCTGTGcc                                                               >  2:1283690/1‑143 (MQ=255)
                                                                              aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCCTCCGCAGGAAGATGTCGTCATCCTGTGCCGCt                                                            >  1:44064/1‑143 (MQ=255)
                                                                               gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGGAGCTGGCGTACTCCTGTGCCGCtc                                                           >  1:1198788/1‑142 (MQ=255)
                                                                               gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                           >  2:945135/1‑143 (MQ=255)
                                                                                   ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCGCCGGACGAGGTCGTAATCCTGTGCCGCTTTcgg                                                       >  2:1016044/1‑141 (MQ=255)
                                                                                   ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGATCGCAGAACTAAAACAACAACAGCAGGAAGAAATCGTAATGCTGTGCCGCTTTCAg                                                       >  2:561793/1‑143 (MQ=255)
                                                                                     ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTAATCCTGTGCCGCATTCAGtt                                                     >  2:726167/1‑143 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTCGCCTCTTCTAATTTCCACGGTTCGCCGTACTCCACCCCCATCCGCACTATGATGTCGTAATCCTGTGCCGCTCTCAGTTc                                                    >  1:151808/1‑143 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   >  2:1004732/1‑143 (MQ=255)
                                                                                         cgtgttttCTGTGACGTGGGCATCTGCGTTTTTGAGTGACTTATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                 <  1:445054/139‑1 (MQ=255)
                                                                                          aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTCATTGCTCTGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                                <  1:566248/143‑1 (MQ=255)
                                                                                             ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                             <  2:1017549/143‑1 (MQ=255)
                                                                                                tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                          <  2:709213/143‑1 (MQ=255)
                                                                                                tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCACGACGCTGTCGTAATCCTGTGCCGCTCTCCCCTCGGCATCGcac                                          >  1:620640/1‑143 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCCACCACCATCCGCAGGCAGATGTCTTAATCCCCTGCCGCTTTCAGTTCTGTATCGaacaaa                                       >  2:766023/1‑143 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                           <  2:392674/107‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                           >  1:392674/1‑107 (MQ=255)
                                                                                                    tGACGAGGTGATCTGCATTTCACCGTGAGTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaac                                      >  2:1060238/1‑143 (MQ=255)
                                                                                                          ggTGTTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCGTACTCCTCCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  1:320261/143‑1 (MQ=255)
                                                                                                            tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAGGATGTCGTAATCCTGTGCCGCCTGCAGTGCGGCATCGACCAAACACCCTTTg                              >  1:7903/1‑143 (MQ=255)
                                                                                                                     tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTAATCCTGTGCCGCATTCAGTTCGGTATCGAACAAACACACTTTGAAACGCCAt                     >  1:541172/1‑143 (MQ=255)
                                                                                                                          ccGTGATGGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCAGGGTGCGCCGTACGCCACCGCCATCCGCAGTAAGATGTCGTAATCCGGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa                <  1:520296/143‑1 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCAATGCCCAg               >  1:74247/1‑143 (MQ=255)
                                                                                                                               attgatATGCCGGCTCTGCGTCGCCTCTTCTACTTCCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGCCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGACCAGGtct           <  2:1158521/143‑3 (MQ=255)
                                                                                                                                  gactggCAGGCTCTGGCGCGCCTCTTCTACTGTCCACGGTTCGGCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACGCTTTGAAACGCCATGCCCAGGTAAAGa        <  1:215868/138‑1 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  >  2:1106966/1‑143 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACCTCCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  <  2:936502/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925

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