Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,180 C→T 25.5% intergenic (+238/‑896) USA300HOU_0072 → / → USA300HOU_0075 universal stress protein/lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,1800CT25.5% 130.6 / 33.5 64intergenic (+238/‑896)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (8/8);  total (32/32)
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.25e-01

ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGA  >  CP000730/76053‑76312
                                                                                                                               |                                                                                                                                    
aCTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTACAGTATAAttt                                                                                                                         <  2:269923/141‑1 (MQ=255)
   aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                      <  2:389144/141‑1 (MQ=255)
   aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                      <  2:22567/141‑1 (MQ=255)
       ttGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATGTAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaa                                                                                                                  >  2:205133/1‑141 (MQ=255)
        tgcttcACCGTTATCACTTATAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAAGAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                 <  1:104934/141‑1 (MQ=39)
        tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                 >  1:42667/1‑141 (MQ=255)
         gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                                >  1:27326/1‑141 (MQ=255)
         gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                                >  1:300480/1‑141 (MQ=255)
             cactggTATCACTTAGAACAGCTAAACAAAGAGTTAGGAATGCAGGGAAGACTAAAACAAGGGCAGTTAATTTAATACGTATAGCAAAATAAAGGTGCTGTGGCAAGGTTGAATCAATAGTATAATTTTAACAGAAAGGAg                                                                                                            >  2:255948/6‑141 (MQ=255)
             cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAg                                                                                                            <  1:17802/141‑1 (MQ=255)
                cTGTACCATATAGAGCAAATAACAAAAGTGTTAGTAATCCAGGGAATACTTACAGAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTATGTTGCTAAGTTGAATCTATAGTATCATTTTAACAACAAGAAGtct                                                                                                         >  2:104935/1‑141 (MQ=255)
                  gTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATCTAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                       >  1:308542/1‑141 (MQ=255)
                  gTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTCGAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                       >  2:170610/1‑141 (MQ=255)
                  atATCACTAAAAACAGCTAACCAAAATTTTAGAAATCCAGGGAATACTAACACAAAGACAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTCGCAAAGTTGAATCTATAGTATAATATTAACAAAAAAGAAtcttc                                                                                                       <  2:52992/140‑1 (MQ=255)
                    aTCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTg                                                                                                     >  2:151423/1‑141 (MQ=255)
                        tttAGAACAGCTAACCAAAATGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATg                                                                                                 <  1:382727/140‑1 (MQ=255)
                        cTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATg                                                                                                 <  1:43621/141‑1 (MQ=255)
                                   tAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAg                                                                                      >  1:367886/1‑141 (MQ=255)
                                           gtgtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAc                                                                              >  2:288862/1‑141 (MQ=255)
                                                aGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAACCAAttt                                                                         >  1:169217/1‑141 (MQ=255)
                                                   aaTCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAc                                                                      <  2:978/141‑1 (MQ=255)
                                                   aaTCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAc                                                                      <  1:19403/141‑1 (MQ=255)
                                                       cAGGGAATACTAACACAAAGGCAATTAATTTAATAATGATAGCAAAATATAGGGTATGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAAGAAATTCACAAgg                                                                  >  2:237099/1‑141 (MQ=21)
                                                                        aaGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGCCTTCTGTATGAACTATTTCAGCTATAAACAATTTACCAAGGATGTTATCACTGTAGcc                                                 >  2:305220/1‑141 (MQ=255)
                                                                         aGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCg                                                <  1:14287/141‑1 (MQ=255)
                                                                          ggCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGt                                               <  1:365392/141‑1 (MQ=255)
                                                                            cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGAAATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                             >  2:20760/1‑141 (MQ=255)
                                                                             aattaattTAATAATTATAGCAAAATATAGGTTCTGTTGAAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTgg                                            <  2:79275/141‑1 (MQ=255)
                                                                              attaattTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGc                                           >  2:258549/1‑141 (MQ=255)
                                                                                 aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                        <  2:42667/141‑1 (MQ=255)
                                                                                   tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                      >  2:33105/1‑141 (MQ=255)
                                                                                   tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                      <  1:135040/141‑1 (MQ=255)
                                                                                     taataaGTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATc                                    <  1:204806/141‑1 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTactaacta                                  >  2:84955/1‑141 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                  >  1:360595/1‑141 (MQ=255)
                                                                                       ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                  >  1:145166/1‑141 (MQ=255)
                                                                                         aatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTTTAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAg                                <  1:192098/141‑1 (MQ=255)
                                                                                             ataGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGat                              <  2:332479/139‑1 (MQ=255)
                                                                                             ataGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGat                              >  1:332479/1‑139 (MQ=255)
                                                                                               aGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGc                          >  1:327654/1‑141 (MQ=255)
                                                                                                 cagttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAt                        >  1:316209/9‑141 (MQ=11)
                                                                                                 cAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAt                        <  2:300480/141‑1 (MQ=255)
                                                                                                 agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                        >  1:245167/9‑141 (MQ=1)
                                                                                                    ttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                     <  2:108980/136‑1 (MQ=2)
                                                                                                    aaattcggTTTTGGTGCAATGTTGAATATATAATATAATTTTAACAAAAAGGAGTCTGCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                     <  2:315574/135‑1 (MQ=11)
                                                                                                    aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                     <  2:128291/136‑1 (MQ=9)
                                                                                                     taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                    >  2:194928/5‑141 (MQ=0)
                                                                                                     atatAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                    <  1:367491/141‑1 (MQ=40)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                    >  2:271317/5‑141 (MQ=2)
                                                                                                     aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                    <  2:351642/137‑1 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                   <  1:56602/138‑1 (MQ=2)
                                                                                                      atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                   >  2:251580/4‑141 (MQ=2)
                                                                                                      aggaGGTTCTGTTGCAACGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                   >  1:312135/4‑140 (MQ=255)
                                                                                                      aggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                   >  2:238224/4‑140 (MQ=0)
                                                                                                        gaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGtaa                 >  1:71858/2‑139 (MQ=0)
                                                                                                        gaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgcta                 <  2:56276/140‑4 (MQ=2)
                                                                                                           gTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg              <  1:63373/141‑1 (MQ=35)
                                                                                                            ttatGTTGCAAAGGTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGt             <  2:249975/138‑1 (MQ=11)
                                                                                                                gTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTCTTTCAGATATAACCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata         <  2:567/141‑1 (MQ=25)
                                                                                                                 tttcAAAGTTCAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACACGCATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatat        <  1:332073/138‑1 (MQ=11)
                                                                                                                  tGCAAAGCTGAACGGATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCGGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatata       >  2:232916/1‑141 (MQ=11)
                                                                                                                    cAAAGATGAATCTATAGTATAATTTTAACAAGAAGGAGTCCTCTGTATGAACTATTTCAGATATAAGCAATTTAACAAGGATGTAATCACTGCAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATc     >  1:293818/1‑141 (MQ=11)
                                                                                                                      aaGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTGCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTg   <  1:269560/141‑1 (MQ=255)
                                                                                                                       aGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCGTCTGTATGAACTAGTGCAGATATAAACGATTTAACAAGGATGTTATCACTGTAGTCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGa  <  1:151423/141‑1 (MQ=11)
                                                                                                                               |                                                                                                                                    
ACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGA  >  CP000730/76053‑76312

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