Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,921,549 T→A 9.5% intergenic (+24/+37) PP_1720 → / ← PP_1721 Alcohol dehydrogenase, zinc‑containing/putative phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,921,5490TA 9.5% 127.4 / 4.7 42intergenic (+24/+37)PP_1720/PP_1721Alcohol dehydrogenase, zinc‑containing/putative phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (15/23);  new base A (2/2);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CAACCGTGATGACGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCAGCCATTCCAGCAACCGCGTG  >  AE015451/1921327‑1921769
                                                                                                                                                                                                                              |                                                                                                                                                                                                                            
cAACCGTGATGACGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                                                                                                 >  1:30150/1‑236 (MQ=255)
 aaCCGTGATGACGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgc                                                                                                                                                                                                                       >  2:276900/1‑229 (MQ=255)
      tgatgaCGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCg                                                                                                                                                                                                             >  2:566925/1‑234 (MQ=255)
          gaCGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGcac                                                                                                                                                                                                         <  2:685673/234‑1 (MQ=255)
                     ggCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                                                                                                      >  2:542696/1‑210 (MQ=255)
                            cTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                                                                                                        >  2:791198/1‑233 (MQ=255)
                                tGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcag                                                                                                                                                                                      <  2:504565/231‑1 (MQ=255)
                                      aaCTGCAGCAGCAAGTGTGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGca                                                                                                                                                                                                          <  1:735027/205‑1 (MQ=255)
                                              cagcaAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCg                                                                                                                                                                      <  2:701875/233‑1 (MQ=255)
                                                     tATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAtt                                                                                                                                                            <  1:552044/236‑1 (MQ=255)
                                                       tGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgc                                                                                                                                                                                                                       <  2:437157/175‑1 (MQ=255)
                                                                ttCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCa                                                                                                                                                  <  1:276902/235‑1 (MQ=255)
                                                                 tCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGGCCCTTCGCGGGTTTACCCGCGAAGGGGCCg                                                                                                                                                                                                             >  2:375638/1‑175 (MQ=255)
                                                                        ggggCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAAcc                                                                                                                                                                                           <  1:566927/186‑1 (MQ=255)
                                                                                 gTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTc                                                                                                                                                                       <  2:293653/197‑1 (MQ=255)
                                                                                                    acCTATCCGGTGGAGTTTGCCCAGGTGGCGTATGCCGAGCTGGAGACCAACCAGGGGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTc                                                                                                                                                                       <  2:10391/178‑1 (MQ=255)
                                                                                                      cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGTTTATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                                                                                                      >  2:633143/1‑129 (MQ=255)
                                                                                                             gTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATcc                                                                                                     <  1:627474/235‑1 (MQ=255)
                                                                                                              tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAccc                                                                                                                                                                                                                         <  1:107612/118‑1 (MQ=255)
                                                                                                                            gcggcgTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg                                                                                     <  1:542698/236‑1 (MQ=255)
                                                                                                                               gcgTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCAGcccc                                                                                                                                                                                                        <  2:647262/118‑4 (MQ=255)
                                                                                                                                cggaTGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGt                                                                                                          <  2:465717/208‑1 (MQ=255)
                                                                                                                                           cTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcag                                                                                                                                                                                      <  2:502924/124‑1 (MQ=255)
                                                                                                                                                      ccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacc                                                                                                                    <  2:372326/179‑1 (MQ=255)
                                                                                                                                                            gtCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTc                                                     >  1:509986/1‑236 (MQ=255)
                                                                                                                                                             gcGGGCCAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCCTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCg                                                                      <  1:110530/217‑1 (MQ=255)
                                                                                                                                                                   aaGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATACAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAacg                                                                                                                                     >  2:524623/1‑148 (MQ=255)
                                                                                                                                                                            agggtgatCGACCCGAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGt                                                                                                          <  2:30150/165‑1 (MQ=255)
                                                                                                                                                                                           aGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCa                      <  1:791200/236‑1 (MQ=255)
                                                                                                                                                                                               tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTc                                                                       >  1:115745/1‑183 (MQ=255)
                                                                                                                                                                                                 gtgtAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGa                                                                     <  2:115744/183‑1 (MQ=255)
                                                                                                                                                                                                    tAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGCAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCCCGCGCAACACTGCGTCCGGGGTCACCCCATGCggg                                                                                                           >  2:238835/1‑142 (MQ=255)
                                                                                                                                                                                                    tAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCAGCCATTcc              >  1:706907/1‑235 (MQ=255)
                                                                                                                                                                                                    tAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCACTCCAgg                                         >  1:491708/1‑208 (MQ=255)
                                                                                                                                                                                                     aaGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCCACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCCGCCATTCCAg            >  1:19787/1‑236 (MQ=255)
                                                                                                                                                                                                         ggggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                                                                                                           >  2:696220/1‑137 (MQ=255)
                                                                                                                                                                                                              tGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCAGCCATTCCAGCAACcgcg    >  1:519137/1‑235 (MQ=255)
                                                                                                                                                                                                               ggCCCCTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCAGCCATTCCAGCAACCGCGTg  >  1:407404/1‑236 (MQ=255)
                                                                                                                                                                                                                   ccTTCGCGGGTTTACCCTCGAAGGGCCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCgg                                                                                                                           <  2:673601/111‑1 (MQ=255)
                                                                                                                                                                                                                    cTTCGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGt                                                                                                          <  2:54680/127‑1 (MQ=255)
                                                                                                                                                                                                                     tgcgcgGGTAAACCCGCGAAGGTGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCa                                                                                                                                                  <  2:15229/84‑1 (MQ=255)
                                                                                                                                                                                                                        gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTcc                                            <  2:680244/185‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                            
CAACCGTGATGACGGCTTCAAGGCCGAGCTGTTGCGCGAACTGCAGCAGCAAGTATGGCCGCTGTTCGCCGAGGGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCTCGAATCGGAATAGAAACTCGCGCCCTCCAGGTTTTCCTGCTCCTGATCCAGCCATTCCAGCAACCGCGTG  >  AE015451/1921327‑1921769

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