Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,199,222 C→A 10.7% F80L (TTC→TTA)  USA300HOU_1127 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,199,2220CA10.7% 169.7 / 2.2 63F80L (TTC→TTA) USA300HOU_1127hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base A (3/5);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTATTAGAGAAACTAAGGTAGGTTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTGAATATAAGGGCAAT  >  CP000730/1199083‑1199353
                                                                                                                                           |                                                                                                                                   
cAAGTATTAGAGAAACTAAGGTAGGTTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTc                                                                                                                                     <  2:524455/140‑1 (MQ=255)
 aaatagtaaggagaCTATGGTAGGTTGTTTTTTAGGCAAAATATATGAACCTTCGTTACAACCATTTAGAAAAGATATTCCACCTATTTTAAGTATCGACATATCATCAATCGCTGCAATGATCGTTTTAGTATTATTcc                                                                                                                                    <  2:120860/127‑1 (MQ=255)
       tAGAGAAACTAAGGTAGGTTATTTTTTAGCGAAGATATAGGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAg                                                                                                                              >  2:424681/1‑140 (MQ=255)
       tAGAGAAACTAAAGTAGGTTATTATTTAGCGAAAATATAGGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAg                                                                                                                              >  2:424643/1‑140 (MQ=255)
           gaAACTGAGGTAGGTTATTTGTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGtt                                                                                                                          <  2:55573/140‑1 (MQ=255)
           gaAACTAAGGTAGGTTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGtt                                                                                                                          <  2:347819/140‑1 (MQ=255)
               cTAAGGTAGGTTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTc                                                                                                                      >  2:369794/1‑140 (MQ=255)
                  aggtaggtTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCaa                                                                                                                   >  1:712186/1‑140 (MQ=255)
                         ttATTTTTTAGCGAAAATATATGAACCTTTCTTACATCCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTa                                                                                                            <  1:53873/140‑1 (MQ=255)
                           aTTTTTTAGCGAAAATATATGATCCTTTCCTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATTTCCTCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTCTTCCAATTATTTCAt                                                                                                          >  2:6401/1‑140 (MQ=255)
                           aTTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTACACCTCTTGGAATTATCGACATATCATCAATCGCTGCAATTATCGCTTTAGTATTATTACACAAAGGGTTACACCAAATCTTTAAt                                                                                                          >  2:674551/1‑140 (MQ=255)
                             tttttAGCGAAAATATATTAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTg                                                                                                        <  1:291359/140‑1 (MQ=255)
                             tttttAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTg                                                                                                        <  2:594795/140‑1 (MQ=255)
                               tttAACGAAAATTTATAAACTTTTCTTATAACCAAAGAGAAAAGTAATTCCTCCTATTGCAATTACGAACATATCATCAATTGTTGCTATTTTCGTGTTAAAAAAAAAAAAAAAAGGGTTACTCCAAATCTTTAATTGGa                                                                                                      <  2:257428/140‑1 (MQ=255)
                               tctAGCAAAAATATATAAACCTTTCTTACAACCATTTAAAAAAGCAATTCCACCTATTAGAATTATCAACATATCATCAATCGCCGCAATTATCGTTTTAAAAAAAAAAAAAAAAGGGTTACTCCAAATCTTTAATTGGa                                                                                                      <  2:377175/138‑1 (MQ=255)
                                  aGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGCCCTATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAAGATTACTAAAAATCTTTAATTAGAttt                                                                                                   >  1:649424/1‑140 (MQ=255)
                                           atatGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAAt                                                                                          <  2:444163/140‑1 (MQ=255)
                                           atatGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATGGGTATTATCGACATGTCAGCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTTGATTTTAATTCAAt                                                                                          >  2:718908/1‑140 (MQ=255)
                                            tatGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTaacacaaga                                                                                         >  1:708787/1‑133 (MQ=255)
                                            tatGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGAAttaaaatccaatt                                                                                         <  1:322072/140‑12 (MQ=255)
                                             atGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTa                                                                                        <  2:211523/140‑1 (MQ=255)
                                              tGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGAGATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTAc                                                                                       >  2:259037/1‑140 (MQ=255)
                                              tGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTAc                                                                                       >  2:192735/1‑140 (MQ=255)
                                              tGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTAc                                                                                       >  2:206650/1‑140 (MQ=255)
                                              tGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGAGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTaa                                                                                       >  1:156570/1‑139 (MQ=255)
                                               gAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACa                                                                                      >  1:469625/1‑140 (MQ=255)
                                                   cTTTCTTACACCCCTTTAGATCAGCAATTCCACATATTGGAATTATCGACATATCATTAATCGCTGTAATTATCGTTTTAGTAATATTACAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATaa                                                                                  <  1:195459/140‑1 (MQ=255)
                                                   cTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATaa                                                                                  <  1:711572/140‑1 (MQ=255)
                                                   cTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATaa                                                                                  <  1:276714/140‑1 (MQ=255)
                                                   cTTTCTTACAACCATTTAGAAAACTTATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATaa                                                                                  <  2:616988/140‑1 (MQ=255)
                                                           cAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAATAATATAAAAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCg                                                                          <  2:437583/140‑1 (MQ=255)
                                                              ccATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATGAATTAAGCGACa                                                                       >  2:186100/1‑140 (MQ=255)
                                                              ccATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACa                                                                       >  2:685149/1‑140 (MQ=255)
                                                                 tttAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATaa                                                                    >  1:649547/1‑140 (MQ=255)
                                                                   tAGAAAAGTAATTCCACCTATTGGAATTATCGACATCTCATCAAACGCTGCAATTATGGTTTTAGTATTATTCCAAAAAGGGTTACCGAAAATCTTTAATTGGATTTTAATTAAATTACAATAAATTAAGCGACATAACg                                                                  >  2:184132/1‑140 (MQ=255)
                                                                    aGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGa                                                                 >  2:139645/1‑140 (MQ=255)
                                                                    aGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGa                                                                 >  2:419871/1‑140 (MQ=255)
                                                                    aGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGa                                                                 >  2:714643/1‑140 (MQ=255)
                                                                    aGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTACAAAAAGGATTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGAAATAACGa                                                                 >  2:139154/1‑140 (MQ=255)
                                                                      aaaaGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTTCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAgg                                                               <  1:323971/140‑1 (MQ=255)
                                                                      aaaaGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAgg                                                               >  2:100998/1‑140 (MQ=255)
                                                                        aaGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCa                                                             <  1:420706/140‑1 (MQ=255)
                                                                         attAATTCCACCTATTGGAATTATCTACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAATGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCaa                                                            <  2:31576/138‑1 (MQ=255)
                                                                         aGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCaa                                                            <  1:369794/140‑1 (MQ=255)
                                                                            aaTTCCACCTATTTGAATTATGGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAggg                                                         <  2:591815/140‑1 (MQ=255)
                                                                               tCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATACCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACTATAAATTAAGCGACATAACGAGGCAAGGGTAc                                                      >  2:406304/1‑140 (MQ=255)
                                                                                   ccTATTGGCATTATCGACATATCATCAATCGTTGCAATTATCGTTTTAGAATAAATAAAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGa                                                  <  2:572806/140‑1 (MQ=255)
                                                                                   ccTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGa                                                  <  1:146403/140‑1 (MQ=255)
                                                                                   ccTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGAGACATAACGAGGCAAGGGTACATGa                                                  >  1:230449/1‑140 (MQ=255)
                                                                                    cTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGAt                                                 >  1:4821/1‑140 (MQ=255)
                                                                                         ggAATTATCGACATATGATCAATCGCTGCAATTATCGTTTTAGTATTATTGCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGGAAAGGAAACttt                                            >  2:196043/1‑140 (MQ=255)
                                                                                           aaTTATCAACATATCATAAATCGCTGCAATAAGCGTGTTAGTATTATTCAAACTAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATAAAGCGACATAGCGAGGCAAGGGAACATGAGACTTTAg                                          >  2:261433/1‑140 (MQ=255)
                                                                                              tATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCt                                       <  1:521523/140‑1 (MQ=255)
                                                                                                acGACATATCATCAATCGCTGCAATTATCGTTTTGGTATTACTCCAACAAGCGATGATTCAAATAATTAAAGAGATTTTAATCCACTTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTCGCCtgg                                     >  2:180926/2‑138 (MQ=255)
                                                                                                   aCATATCATCAATCGCTGCAATTATCGCTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAACTTAAGCGACATAACGCGGCAAGGGTACATTACACTTTAGCCCCGttt                                  >  2:10353/1‑140 (MQ=255)
                                                                                                             aaTCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATATCGTGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCa                        <  1:192735/140‑1 (MQ=255)
                                                                                                                   tGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAATCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTc                  >  1:214160/1‑140 (MQ=255)
                                                                                                                     cAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTg                <  1:419871/140‑1 (MQ=255)
                                                                                                                     cAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTg                <  1:714643/140‑1 (MQ=255)
                                                                                                                     cAATTATCGTTTTAGTATTATTCCAAAAAGGGCTAATAGTAACCTTTAACTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTg                >  1:716831/1‑140 (MQ=255)
                                                                                                                       aTTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTATTCTGaa              >  2:572616/1‑140 (MQ=255)
                                                                                                                       aTTATCGTTTTAGTATTATTACAAAAAGGGTTACTCCATATCCTTAATTGGATTTTAAGTCAATTACAATAAATTAAGAGACATAACGAAGCAAGGGAACATGATACTTTAGCCTCGTTTTTGATATGCATTGTGCTGac              >  2:98033/1‑139 (MQ=255)
                                                                                                                                  aGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACCTGATACTGTAGCCTCGTTTTTGATATGCATTTTTCTGAATATAAGGGCaa   <  1:384464/140‑1 (MQ=255)
                                                                                                                                   gTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTGAATATAAGGGCAAt  <  1:629636/140‑1 (MQ=255)
                                                                                                                                   gTATTATTCCAAAAAGGGCTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATCCATTTTTCTGAATATAAGGGCAAt  <  1:463079/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CAAGTATTAGAGAAACTAAGGTAGGTTATTTTTTAGCGAAAATATATGAACCTTTCTTACAACCATTTAGAAAAGTAATTCCACCTATTGGAATTATCGACATATCATCAATCGCTGCAATTATCGTTTTAGTATTATTCCAAAAAGGGTTACTCCAAATCTTTAATTGGATTTTAATTCAATTACAATAAATTAAGCGACATAACGAGGCAAGGGTACATGATACTTTAGCCTCGTTTTTGATATGCATTTTTCTGAATATAAGGGCAAT  >  CP000730/1199083‑1199353

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