Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,397,8300TC44.2% 69.1 / 12.6 55I265V (ATT→GTT) uspEstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (21/8);  new base C (0/23);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 4.14e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGCT  >  NC_000913/1397680‑1397980
                                                                                                                                                       |                                                                                                                                                      
cttcttCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAAt                                                                                                                                                        >  2:979756/1‑151 (MQ=255)
     tcGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATcccc                                                                                                                                                    <  2:473732/150‑1 (MQ=255)
       gTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGcc                                                                                                                                                 >  2:877759/1‑151 (MQ=255)
        tCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCt                                                                                                                                                >  2:373708/1‑151 (MQ=255)
         cATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCt                                                                                                                                                <  1:271666/150‑1 (MQ=255)
         cATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTg                                                                                                                                               <  1:455448/151‑1 (MQ=255)
          aTCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGt                                                                                                                                              >  1:665948/1‑151 (MQ=255)
           tCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGt                                                                                                                                              >  1:583337/1‑150 (MQ=255)
            ccAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTa                                                                                                                                             >  1:694963/1‑150 (MQ=255)
                ttaaaggggTGTCTGATACTGGTCAGGTTTAATGCCCAGCAGGTCGCAGCAAAGATGACAAATCCCCTGTCCCCCCGTTTTGCCGAGGAATGCTGCTAAAATACCGGTGCGTCTTCCCGTGCCAAGACCCACAACCCCCGCCTGTAAGTg                                                                                                                                         <  2:329716/144‑1 (MQ=255)
                 tCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCt                                                                                                                                       >  2:1076093/1‑151 (MQ=255)
                 tCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCt                                                                                                                                       >  2:246514/1‑151 (MQ=255)
                      ggttgtCTGATACTGGTCAGGTTAAATGACCAGCAGGTCGCAGCGAAGATGATAAATCCCCTGTTCCGCCGTGTTGCGGAGGAATGCTGCTAAAATACCGGTGCGTCTTCCCGTCCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGcc                                                                                                                                  <  1:959384/148‑1 (MQ=255)
                           tCTGATACTGGTCAGTTTTAATGACCAGCAGGTCGCAGGGAAGATAATCAACCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGCCCTCCCGTGCCCAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCa                                                                                                                                 <  1:542/147‑1 (MQ=255)
                                 acgggTCATGTTAAAGGACCGGCGGGCCGCAGCAAAGATGATCAACCACTTGTTCCCCCGTGTGGCGGAGGAATGCTGCTGAAATCCGGTGGCGTCCTACCGGGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGaa                                                                                                                        <  1:891658/147‑1 (MQ=255)
                                 aCTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATc                                                                                                                             <  1:592284/145‑1 (MQ=255)
                                     gTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCAcc                                                                                                                   >  2:617667/1‑151 (MQ=255)
                                      tCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACct                                                                                                                  >  1:155301/1‑151 (MQ=255)
                                       cAGGTTAAATAACCGGCAGGTCGCAGCGAAAATGATAAATCCCCTGTCCCGCCGTGTTGCCGGGGAATGCTGCTGAAATCCCGGTGCGCCCTACCGTGCCCAAAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACct                                                                                                                  <  1:140916/150‑1 (MQ=255)
                                           tttAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtcttct                                                                                                             >  1:424964/1‑151 (MQ=255)
                                              aaTGCCCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAAGGCTGCTAAAATACCGGTGCGCCCTACCGTCCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                           <  1:892207/150‑1 (MQ=255)
                                                    cagcagGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCt                                                                                                    >  2:363513/1‑151 (MQ=255)
                                                        agGCCGCAGCGAAGATGATCAACCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTAAAATACCGGTGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCt                                                                                                    <  1:456593/147‑1 (MQ=255)
                                                         ggtcgcggagagattgatCAACCCCCTTTCCCCCCGTGTTGCCGGGGATTCTTGTTGAAACCCCGGTGCACCCTACCGCCCCCACACCCCCACGCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCttt                                                                                                  <  2:509796/134‑1 (MQ=255)
                                                               aGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTg                                                                                         >  1:557949/1‑151 (MQ=255)
                                                                   agaatgaccatcccccTGTCCCGCCGGGTTCCCGAGAAATGCTCTTAAAATCCGGGTGCTCCCTCCCGTCCCCAAACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACg                                                                                           <  2:571115/131‑1 (MQ=255)
                                                                   aagtctaacaatcgCTTGTCCCCCGGTTTTGCCCAGGATTCCTGCTGAAATCCCGGTGCGTCCTACCGTCCCTAGAACCCAACGCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTc                                                                                     <  2:1074031/143‑1 (MQ=255)
                                                                   aaGATGATCAATCCCCTTTTCCGCCGTGTTGCCGAGGAATGCTGCTAAAATCCCGGTGCGTCCTACCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGc                                                                                        <  1:4094/148‑1 (MQ=255)
                                                                   aaGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTc                                                                                     >  1:453245/1‑151 (MQ=255)
                                                                    aGATGATCACTCCCCTGTCCCGCCGTGTTGCCGAGAAATGCTGCTAAAATCCCGGGCCGCCCTACCGTCCCAAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCa                                                                                    <  2:39372/151‑1 (MQ=255)
                                                                    aGATAAACAATCCCCTGTTCCGCCGTTTTGCCGAGGAATGCTGCTAAAATCCCGGTGCTTCCTCCCGTGCCCAAACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCg                                                                                       <  2:332367/148‑1 (MQ=255)
                                                                     gatgatCAACCCCCTGTCCCGCCGTGTGGCCGAGGAATGCTGCTAAAATCCCGGTGCGTCCTCCCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCAt                                                                                   <  1:779604/151‑1 (MQ=255)
                                                                         atcaatccCCTGTTCCGCGGTGTGGCCGAGGAATGCTGCTAAAATCCCGGTGCGTCCTCCCGTGCCCAGACCCCCAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGttt                                                                               <  2:827177/151‑1 (MQ=255)
                                                                                    gTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGcc                                                                    >  2:1001060/1‑151 (MQ=255)
                                                                                          ccgTGTTGCCGAGGAATGCTGCTAAAATCCCGGTGCGTCCTCCCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAtt                                                               <  2:433648/150‑1 (MQ=255)
                                                                                            gttttGCCAGGGATTCCTGCTAAAATACCGGTGCGTCCTACCGTCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTc                                                             <  2:89179/147‑1 (MQ=255)
                                                                                            gtgtTGCCGAGAAATGCTGCTAAAATACCGGTGCTTCCTCCCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTc                                                             <  2:225325/150‑1 (MQ=255)
                                                                                              gtTGCCGAGGAATGCTGCTGAAATCCCGGTGCGTCCTACCGTGCCAAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCt                                                            <  1:708314/149‑1 (MQ=255)
                                                                                                tGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCg                                                        >  1:802601/1‑151 (MQ=255)
                                                                                                    gAGGAATGCTGCTGAAATCCCGGTGCGTCCTACCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAgg                                                    <  1:479017/151‑1 (MQ=255)
                                                                                                    gAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGAGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAgg                                                    >  2:126767/1‑151 (MQ=255)
                                                                                                      ggAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGc                                                  >  1:606717/1‑151 (MQ=255)
                                                                                                                  aaaTACGGTTGCGTCCTACCGTCCCAAAAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCa                                       <  1:626077/150‑1 (MQ=255)
                                                                                                                   aaTACCGGTGCGTCCTACCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGc                                     <  1:382225/151‑1 (MQ=255)
                                                                                                                               gccttACCGTGCCCAGAACCAAATTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCaa                                   <  1:172946/138‑1 (MQ=255)
                                                                                                                               gTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGccc                          >  1:187570/1‑151 (MQ=255)
                                                                                                                                    aCCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGa                      >  1:366372/1‑150 (MQ=255)
                                                                                                                                      cGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATg                   >  2:396960/1‑151 (MQ=255)
                                                                                                                                         gCCCACACCCACATGCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAAt                    <  1:962728/147‑1 (MQ=255)
                                                                                                                                          cccAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCa                <  1:218089/150‑1 (MQ=255)
                                                                                                                                          cccAAACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCAt               <  1:589465/151‑1 (MQ=255)
                                                                                                                                          cccAAAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATg                   <  1:402652/147‑1 (MQ=255)
                                                                                                                                           ccAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATc              <  1:519701/151‑1 (MQ=255)
                                                                                                                                               acccACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCg             <  2:899628/146‑1 (MQ=255)
                                                                                                                                                 ccACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTtat         <  1:345441/150‑1 (MQ=255)
                                                                                                                                                  cacaATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTtata        >  2:407880/1‑150 (MQ=255)
                                                                                                                                                       tCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGCt  >  1:543730/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTATAAACGCT  >  NC_000913/1397680‑1397980

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