Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,377,769 G→T 6.7% intergenic (‑349/+46) sspA ← / ← rpsI stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,7690GT 6.7% 193.2 / 2.1 57intergenic (‑349/+46)sspA/rpsIstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base T (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCA  >  NC_000913/3377627‑3377915
                                                                                                                                              |                                                                                                                                                  
agagATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCgg                                                                                                                                              <  3:134373/149‑1 (MQ=255)
     tGCGCTTTAAACTGAAATAGTTGACTTAGTCCCTTATTGGTTATGTGGTTTTTGTTTTACCTGTCTGTCATGTAGCTGCAAAAAGCAACTTTGCATTTTTTTCGCTGATTCAGGTTTTAGCTATAAAAAAACCCGGCGAAGCGGGtttt                                                                                                                                         <  4:287051/149‑1 (MQ=255)
       cgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                       <  3:57420/149‑1 (MQ=255)
               tCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAAtt                                                                                                                               >  3:39519/1‑149 (MQ=255)
                  gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttt                                                                                                                                          <  2:484750/135‑1 (MQ=255)
                  gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttt                                                                                                                                          >  1:484750/1‑135 (MQ=255)
                  gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGt                                                                                                                             <  1:763689/148‑1 (MQ=255)
                    aaTAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                           >  2:310573/1‑148 (MQ=255)
                        gTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGc                                                                                                                      >  4:494627/1‑149 (MQ=255)
                            aCTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGa                                                                                                                  <  3:518479/149‑1 (MQ=255)
                              ttAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                           >  3:176021/1‑138 (MQ=255)
                              ttAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                           <  4:176021/138‑1 (MQ=255)
                              ttAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGc                                                                                                                >  4:208133/1‑149 (MQ=255)
                                  tCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCTAAAATTGTTTTCTGCCGGAGCAGaa                                                                                                            >  4:679445/1‑149 (MQ=255)
                                  tCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGaa                                                                                                            >  4:898376/1‑149 (MQ=255)
                                    ccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttttt                                                                                                                                        <  3:599061/119‑1 (MQ=255)
                                    ccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttttt                                                                                                                                        >  4:599061/1‑119 (MQ=255)
                                    ccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                          >  4:546402/1‑149 (MQ=255)
                                    ccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                          >  3:482798/1‑149 (MQ=255)
                                      ttATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCa                                                                                                        >  3:256308/1‑149 (MQ=255)
                                       tATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCa                                                                                                        <  1:609231/148‑1 (MQ=255)
                                                tgtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGttt                                                                                              <  4:199124/149‑1 (MQ=255)
                                                  tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTg                                                                                             >  1:685881/1‑148 (MQ=255)
                                                      ttttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgaga                                                                                         >  1:563918/1‑148 (MQ=255)
                                                          gTTTTACCTGTGTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTagc                                                                                                  >  4:201568/1‑133 (MQ=255)
                                                          gTTTTACCTGTCTGTCAGGTGGCAACAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAAc                                                                                                  <  3:201568/135‑1 (MQ=255)
                                                           ttttACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCg                                                                                   >  4:539487/1‑149 (MQ=255)
                                                                ccTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacga                                                                              >  3:285524/1‑149 (MQ=255)
                                                                    tctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTg                                                                                     >  2:40805/1‑138 (MQ=255)
                                                                    tctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTg                                                                                     <  1:40805/138‑1 (MQ=255)
                                                                        tcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGc                                                                      >  3:803037/1‑149 (MQ=255)
                                                                             tGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTAc                                                                 <  4:482798/149‑1 (MQ=255)
                                                                                cagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCa                                                              <  3:539487/149‑1 (MQ=255)
                                                                                  gcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAg                                                             >  2:50857/1‑148 (MQ=255)
                                                                                  gcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGa                                                            <  3:208133/149‑1 (MQ=255)
                                                                                   caAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGa                                                            <  1:310573/148‑1 (MQ=255)
                                                                                   caAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGAc                                                           <  4:39519/149‑1 (MQ=255)
                                                                                    aaaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGAcc                                                          >  3:887821/1‑149 (MQ=255)
                                                                                           aaCTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTc                                                   >  4:309864/1‑149 (MQ=255)
                                                                                              tttCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTa                                                >  4:293593/1‑149 (MQ=255)
                                                                                              tttCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTa                                                <  4:285524/149‑1 (MQ=255)
                                                                                               ttCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTa                                                <  2:563918/148‑1 (MQ=255)
                                                                                               ttCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTAc                                               <  3:494627/149‑1 (MQ=255)
                                                                                                 ccAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACg                                              <  2:685881/148‑1 (MQ=255)
                                                                                                   aGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCgg                                                                                                                                              <  4:351558/50‑1 (MQ=255)
                                                                                                   aGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCgg                                                                                                                                              >  3:351558/1‑50 (MQ=255)
                                                                                                                   cAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                            >  1:699665/1‑148 (MQ=255)
                                                                                                                    aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCTAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCAc                          <  3:679445/149‑1 (MQ=255)
                                                                                                                       ttttAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCGACCTGACGAGCGTCACGAg                       <  3:898376/149‑1 (MQ=255)
                                                                                                                          tAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTa                     <  2:387909/148‑1 (MQ=255)
                                                                                                                               ataAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                       >  4:447683/1‑93 (MQ=255)
                                                                                                                               ataAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                       <  3:447683/93‑1 (MQ=255)
                                                                                                                               ataAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGaa                >  1:146830/1‑148 (MQ=255)
                                                                                                                                taAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGc              >  3:843234/1‑149 (MQ=255)
                                                                                                                                   aaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAg           >  3:340616/1‑149 (MQ=255)
                                                                                                                                     aaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCt         >  3:492805/1‑149 (MQ=255)
                                                                                                                                        cccGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTa      <  3:293593/149‑1 (MQ=255)
                                                                                                                                            ccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCa  <  4:256308/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCA  >  NC_000913/3377627‑3377915

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