Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,980,512 G→T 6.5% Q83K (CAA→AAA)  USA300HOU_1839 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,980,5120GT 6.5% 211.3 / 3.2 61Q83K (CAA→AAA) USA300HOU_1839hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base T (2/2);  total (26/35)
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.44e-01

TTGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAATATAT  >  CP000730/1980378‑1980649
                                                                                                                                      |                                                                                                                                         
ttGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGACGTAAATTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCt                                                                                                                                     <  1:18980/141‑1 (MQ=255)
  gATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCttt                                                                                                                                   >  1:225140/1‑141 (MQ=255)
  gATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCttt                                                                                                                                   >  1:250385/1‑141 (MQ=255)
     ttCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCa                                                                                                                                >  2:127047/1‑141 (MQ=255)
       ccATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAAt                                                                                                                              >  1:234334/1‑141 (MQ=255)
       ccATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAAt                                                                                                                              >  1:234365/1‑141 (MQ=255)
        cATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAAtt                                                                                                                             >  2:77720/1‑141 (MQ=255)
          ttAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCATCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAAtttt                                                                                                                           <  2:442653/141‑1 (MQ=255)
              tttGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCat                                                                                                                       >  2:90929/1‑141 (MQ=255)
              tttGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTAACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCat                                                                                                                       >  1:221507/1‑141 (MQ=255)
                           ttCTAAATCAAGATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTGTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTATTCAATTTTTCATTAATTAACGTTTc                                                                                                          >  2:268284/1‑141 (MQ=255)
                                  tCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGc                                                                                                   <  2:20010/141‑1 (MQ=255)
                                  tCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGc                                                                                                   <  2:286942/141‑1 (MQ=255)
                                   cAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCa                                                                                                  >  2:66794/1‑141 (MQ=255)
                                   cAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCa                                                                                                  >  1:325007/1‑141 (MQ=255)
                                      tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:175239/1‑141 (MQ=255)
                                      tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:36389/1‑141 (MQ=255)
                                      tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:36388/1‑141 (MQ=255)
                                              tGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCa                                                                                       <  1:271141/141‑1 (MQ=255)
                                                agaAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAAt                                                                                     <  2:181048/141‑1 (MQ=255)
                                                    gTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCTCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACt                                                                                 >  1:445986/1‑141 (MQ=255)
                                                        tttGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGc                                                                             >  2:246745/1‑141 (MQ=255)
                                                           gTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCaa                                                                          <  1:44914/141‑1 (MQ=255)
                                                               tATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAAtttt                                                                      <  1:328406/141‑1 (MQ=255)
                                                                          ggtCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:180299/139‑1 (MQ=255)
                                                                          gATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTCAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:14340/141‑1 (MQ=255)
                                                                          gATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:139425/141‑1 (MQ=255)
                                                                            tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:84843/141‑1 (MQ=255)
                                                                            tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  2:368886/141‑1 (MQ=255)
                                                                            tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  2:49058/141‑1 (MQ=255)
                                                                            tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:428226/141‑1 (MQ=255)
                                                                            tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:152077/141‑1 (MQ=255)
                                                                             cTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttctt                                                        >  1:162117/1‑141 (MQ=255)
                                                                                  gTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGt                                                   >  1:227149/1‑141 (MQ=255)
                                                                                   tCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTc                                                  <  1:246522/141‑1 (MQ=255)
                                                                                    cAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCt                                                 >  1:147380/1‑141 (MQ=255)
                                                                                     aaCTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCtt                                                <  1:61208/141‑1 (MQ=255)
                                                                                              aaCTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTAtt                                       <  1:443696/141‑1 (MQ=255)
                                                                                                     gTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTg                                <  2:234334/141‑1 (MQ=255)
                                                                                                     gTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTg                                <  2:234365/141‑1 (MQ=255)
                                                                                                        gAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGgtc                             <  1:436350/141‑2 (MQ=255)
                                                                                                           tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:36389/141‑1 (MQ=255)
                                                                                                           tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:36388/141‑1 (MQ=255)
                                                                                                           tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:175239/141‑1 (MQ=255)
                                                                                                             aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:225140/141‑1 (MQ=255)
                                                                                                             aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:250385/141‑1 (MQ=255)
                                                                                                             aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:399984/141‑1 (MQ=255)
                                                                                                             aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:305760/141‑1 (MQ=255)
                                                                                                              tAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATg                       >  2:174466/1‑141 (MQ=255)
                                                                                                                 tttCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGa                    <  2:369635/140‑1 (MQ=255)
                                                                                                                 aTTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGa                    <  2:172193/141‑1 (MQ=255)
                                                                                                                   tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGGAAATGCTTTTTAGTGACGCGGCGAATCTTGGTTGCAATTTCTTCGGGGCTTTTCGGTATTAAACCTGGTGTAAATGGGAcg                  >  1:344940/1‑140 (MQ=255)
                                                                                                                   tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAACTGCTCTTCAATAACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGAct                  >  1:80319/1‑141 (MQ=255)
                                                                                                                   tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGAct                  >  1:424834/1‑141 (MQ=255)
                                                                                                                     tATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGActct                <  1:249119/141‑1 (MQ=255)
                                                                                                                        tgcttgctGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTaaa             <  2:227149/141‑1 (MQ=255)
                                                                                                                        tgcttgctGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTaaa             <  1:198708/141‑1 (MQ=255)
                                                                                                                               tgtgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCCATCCCTTGGCCACTTCTAGTTGCAATTTCTTCGAGCCTAATAGCTCACAAACAAACTGCAAACGGGACTCTAAATTTAAat      >  1:92672/1‑141 (MQ=255)
                                                                                                                               tgtgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAat      >  2:177039/1‑141 (MQ=255)
                                                                                                                                 tgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAatat    >  1:272533/1‑141 (MQ=255)
                                                                                                                                   aTTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAatatat  <  2:371727/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                         
TTGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAATATAT  >  CP000730/1980378‑1980649

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: