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