Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 469,786 A→G G118G (GGA→GGG)  USA300HOU_0445 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730469,7860AG95.2% 204.4 / ‑1.5 63G118G (GGA→GGG) USA300HOU_0445hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (43/17);  total (44/19)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGAGATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACTT  >  CP000730/469649‑469923
                                                                                                                                         |                                                                                                                                         
cGAGATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCat                                                                                                                                        >  2:275922/1‑141 (MQ=255)
  agaTGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAACGGGGCatta                                                                                                                                      <  1:766421/141‑1 (MQ=255)
         gAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGtt                                                                                                                               <  2:270000/141‑1 (MQ=255)
           aTTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTaa                                                                                                                             >  2:640256/1‑141 (MQ=255)
           aTTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTaa                                                                                                                             >  2:109007/1‑141 (MQ=255)
                  aaGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTc                                                                                                                      <  1:833354/141‑1 (MQ=255)
                   aGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCt                                                                                                                     <  1:85279/141‑1 (MQ=255)
                        gATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTataat                                                                                                                >  1:574361/1‑141 (MQ=255)
                          tAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTAtaataa                                                                                                              >  1:259444/1‑141 (MQ=255)
                             aGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAggg                                                                                                           <  2:701121/141‑1 (MQ=255)
                              gggATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGc                                                                                                          <  2:4493/141‑1 (MQ=255)
                                gATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCaa                                                                                                        <  2:398940/141‑1 (MQ=255)
                                    tgGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTc                                                                                                    >  1:3279/1‑141 (MQ=255)
                                        cGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAg                                                                                                >  1:330882/1‑141 (MQ=255)
                                           tatataCAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGCAGGTATGGACTTATACATAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATa                                                                                             <  1:171661/141‑1 (MQ=255)
                                             tataCAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGa                                                                                           >  2:305901/1‑141 (MQ=255)
                                                 cAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGaaaaa                                                                                       <  2:421144/141‑1 (MQ=255)
                                                        cGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATaa                                                                                >  2:566745/1‑141 (MQ=255)
                                                         gCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATaaa                                                                               >  2:720219/1‑141 (MQ=255)
                                                         gCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATaaa                                                                               >  1:381004/1‑141 (MQ=255)
                                                              aaGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtga                                                                          >  1:371465/1‑141 (MQ=255)
                                                                gTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTagaaa                                                                        >  2:640234/1‑141 (MQ=255)
                                                                  aatttaCAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatg                                                                      >  2:530219/6‑141 (MQ=255)
                                                                  aataaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatg                                                                      >  1:12673/1‑141 (MQ=255)
                                                                   ataaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatga                                                                     >  2:381682/1‑141 (MQ=255)
                                                                   ataaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatga                                                                     >  1:489112/1‑141 (MQ=255)
                                                                    taaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaa                                                                    >  1:675137/1‑141 (MQ=255)
                                                                    taaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaa                                                                    >  1:794101/1‑141 (MQ=255)
                                                                    taaACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaa                                                                    >  1:415105/1‑141 (MQ=255)
                                                                     aaaCAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACGCGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaaa                                                                   >  1:659672/1‑141 (MQ=255)
                                                                     aaaCAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaaa                                                                   >  1:233414/1‑141 (MQ=255)
                                                                     aaaCAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTtgaaatgaaa                                                                   >  2:548355/1‑141 (MQ=255)
                                                                        cAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAaat                                                                <  1:305157/141‑1 (MQ=255)
                                                                         aatgtgtaatGAGTAATGAAAGAATGATCATATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTAATTTGTTAAGACATTCTATAATAAGGGCAAAAAAAAAAAAAAAAAAAATTATAAAGTTGAAATGAAAaata                                                               <  2:728924/132‑1 (MQ=21)
                                                                                        aTGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATAAACAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCttat                                                <  2:584692/141‑1 (MQ=255)
                                                                                           aaGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAg                                             >  1:666046/1‑141 (MQ=255)
                                                                                             ggTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGAt                                           >  2:699537/1‑141 (MQ=255)
                                                                                             ggTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGAt                                           >  2:609697/1‑141 (MQ=255)
                                                                                                 tGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAg                                       >  2:718043/1‑141 (MQ=255)
                                                                                                  ggTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGt                                      <  2:369721/141‑1 (MQ=255)
                                                                                                   gTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTa                                     >  2:274374/1‑141 (MQ=255)
                                                                                                     cTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTaga                                   >  2:392199/1‑141 (MQ=255)
                                                                                                     cTTATACTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAAAAATCTTATTAGATAAAGTaga                                   >  2:392257/1‑141 (MQ=255)
                                                                                                         taCTTAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATGATCTTATTAGATAAAGTAGAAGAt                               >  1:406926/1‑141 (MQ=255)
                                                                                                            ttAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACa                            >  1:725618/1‑141 (MQ=255)
                                                                                                            ttAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATT‑CCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACaa                           >  1:805396/1‑141 (MQ=255)
                                                                                                            ttAGATAGAAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATT‑CCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACaa                           >  1:831274/1‑141 (MQ=255)
                                                                                                             tagatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACaa                           >  1:510800/1‑141 (MQ=255)
                                                                                                              agatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACaaa                          <  2:553742/141‑1 (MQ=255)
                                                                                                               gatagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAAt                         >  1:203321/1‑141 (MQ=255)
                                                                                                                atagaAAAACACGAAATGAAAAGAGGAATTAAATTGTTAAAACAATAAATAAAAAAGACAAATACCAAGATAGAAAAAATTATAAAATAAAAATAAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGAAACaaaac                        <  2:487934/141‑3 (MQ=12)
                                                                                                                 tagaAATACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCt                       >  2:446341/1‑141 (MQ=255)
                                                                                                                    aaaTACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTaaa                    >  2:601087/1‑141 (MQ=255)
                                                                                                                     aaTACACGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTaaaa                   >  2:825679/1‑141 (MQ=255)
                                                                                                                          acGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAagag              <  1:766365/141‑1 (MQ=255)
                                                                                                                          acGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAagag              >  1:743908/1‑141 (MQ=255)
                                                                                                                           cGGACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAgaga             >  1:857858/1‑141 (MQ=255)
                                                                                                                             gACTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAAt           >  2:840212/1‑141 (MQ=255)
                                                                                                                              aCTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATa          >  2:149482/1‑141 (MQ=255)
                                                                                                                               cTGCAAAGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAg         <  1:110824/141‑1 (MQ=255)
                                                                                                                                   aaaGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGaaaa     >  2:613861/1‑141 (MQ=255)
                                                                                                                                     aGGGGCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTAAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACt   <  2:605868/141‑1 (MQ=255)
                                                                                                                                      ggggCATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACtt  <  1:612814/141‑1 (MQ=255)
                                                                                                                                      gggACATTCTTTTGTTAAGACATTATATAAAAAATGAAAAATACAAGATAAAAAAAATTATAAAATTAAAATAAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACtt  <  2:735821/141‑1 (MQ=14)
                                                                                                                                         |                                                                                                                                         
CGAGATGGCGAATTTAAAAAGGGTGATAAAGGGATGTGGACGATATATACAGATTTCGCCAAAAGTAATAAACAAGGTGGATTGAGTAATGAAGGTATGGTCTTATACTTAGATAGAAATACACGGACTGCAAAGGGACATTATTTTGTTAAGACATTCTATAATAAGGGCAAATTCCCAGATAGAAAAAATTATAAAGTTGAAATGAAAAATAATAAAATTATCTTATTAGATAAAGTAGAAGATACAAATCTAAAAAAGAGAATAGAAAACTT  >  CP000730/469649‑469923

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