Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,634,642 C→T G82S (GGC→AGC)  tfaQ ← Qin prophage; putative tail fiber assembly protein TfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,6420CT95.2% 122.6 / 0.4 62G82S (GGC→AGC) tfaQQin prophage; putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (32/27);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCA  >  NC_000913/1634510‑1634773
                                                                                                                                    |                                                                                                                                   
tCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGcc                                                                                                                              <  1:289008/140‑1 (MQ=18)
  cGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCAc                                                                                                                            >  2:547267/1‑140 (MQ=17)
  cGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCAc                                                                                                                            >  2:465920/1‑140 (MQ=17)
         gTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATaa                                                                                                                     >  1:1015793/1‑140 (MQ=17)
          tttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATaaa                                                                                                                    >  1:529457/1‑140 (MQ=17)
          tttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATaaa                                                                                                                    <  2:304994/140‑1 (MQ=18)
                    tCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAcc                                                                                                          <  2:268942/140‑1 (MQ=18)
                     ccGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                         <  1:866324/140‑1 (MQ=18)
                     ccGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                         <  2:546761/140‑1 (MQ=17)
                             ttCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATc                                                                                                 >  2:372461/1‑140 (MQ=17)
                               cACCCAGGCTGTGCCGTTCCACTTCTGAAGCTCCCCTTCCGGCTATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                               >  2:1069228/1‑140 (MQ=12)
                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                               >  2:136191/1‑140 (MQ=17)
                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                               >  2:433914/1‑140 (MQ=17)
                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTCACCCCGATGATCtt                                                                                               >  2:644464/1‑140 (MQ=17)
                                     ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                         >  1:629878/1‑140 (MQ=17)
                                     ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                         >  1:1152860/1‑140 (MQ=17)
                                      gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACgag                                                                                        >  2:1179745/1‑140 (MQ=17)
                                        tgtgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGAt                                                                                      <  1:1179745/140‑1 (MQ=17)
                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                    >  2:15438/1‑140 (MQ=17)
                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                    >  2:84614/1‑140 (MQ=32)
                                                tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                              >  1:269135/1‑140 (MQ=17)
                                                tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                              >  2:204268/1‑140 (MQ=17)
                                                    cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                          >  2:802413/1‑140 (MQ=17)
                                                    cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                          >  1:678075/1‑140 (MQ=17)
                                                    cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                          >  2:949456/1‑140 (MQ=17)
                                                       cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                       <  1:547267/140‑1 (MQ=17)
                                                       cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                       <  1:465920/140‑1 (MQ=17)
                                                       cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                       <  2:1015793/140‑1 (MQ=17)
                                                         gAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                     <  1:1119234/140‑1 (MQ=17)
                                                                  cTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                       <  2:512267/129‑1 (MQ=17)
                                                                  cTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                       >  1:512267/1‑129 (MQ=17)
                                                                  cTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaaa                                                            >  2:829233/1‑140 (MQ=17)
                                                                              aCCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAag                                                <  1:15438/140‑1 (MQ=17)
                                                                              aCCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAag                                                <  2:64163/140‑1 (MQ=17)
                                                                               ccAGGTAACATTTACCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTCCCCCGATGATCTTCAACGAGATGCCACGATGACTCATAACTGTTGAAAACAGCCACGAagc                                               <  1:1069228/140‑1 (MQ=255)
                                                                                 aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                             <  1:829233/140‑1 (MQ=17)
                                                                                 aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                             <  1:1054285/140‑1 (MQ=32)
                                                                                  ggTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                            >  2:228240/1‑140 (MQ=17)
                                                                                     aaCATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCg                                         <  1:315598/140‑1 (MQ=17)
                                                                                     aaCATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCg                                         <  2:269135/140‑1 (MQ=17)
                                                                                          tttCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAata                                    >  1:386152/1‑140 (MQ=17)
                                                                                             ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                 >  2:926316/1‑140 (MQ=17)
                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                               <  2:678075/140‑1 (MQ=17)
                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                               <  1:233013/140‑1 (MQ=17)
                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTAATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                               <  1:204268/140‑1 (MQ=12)
                                                                                                 tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                               <  2:1063469/138‑1 (MQ=17)
                                                                                                 tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                               >  1:1063469/1‑138 (MQ=18)
                                                                                                  aaCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCgg                            >  2:123484/1‑140 (MQ=17)
                                                                                                  aaCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCgg                            >  1:630905/1‑140 (MQ=17)
                                                                                                    cGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTg                          <  1:136191/140‑1 (MQ=17)
                                                                                                     ggACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                         <  1:228240/140‑1 (MQ=17)
                                                                                                     ggACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                         <  1:84614/140‑1 (MQ=32)
                                                                                                      gACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCa                        >  2:1004943/1‑140 (MQ=18)
                                                                                                       aCCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCaa                       <  1:433914/140‑1 (MQ=17)
                                                                                                               cAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTa               >  2:735752/1‑140 (MQ=18)
                                                                                                                        ataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCt      >  1:1205773/1‑140 (MQ=18)
                                                                                                                          aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg    >  2:556769/1‑140 (MQ=18)
                                                                                                                          aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg    >  1:932941/1‑140 (MQ=18)
                                                                                                                           aCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGc   <  1:735752/140‑1 (MQ=17)
                                                                                                                           aCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGc   <  1:1004943/140‑1 (MQ=17)
                                                                                                                            cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa  <  2:1205773/140‑1 (MQ=17)
                                                                                                                            cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa  >  1:1005729/1‑140 (MQ=18)
                                                                                                                            cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa  >  1:533174/1‑140 (MQ=17)
                                                                                                                                    |                                                                                                                                   
TCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCA  >  NC_000913/1634510‑1634773

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