Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,223,369 G→C 10.4% P393R (CCG→CGG)  PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,223,3690GC10.4% 162.6 / 6.0 57P393R (CCG→CGG) PP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base C (3/3);  total (28/29)
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.51e-01

ATAGGAATGGTGGCCGAACAGCAGACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAGAGG  >  NC_002947/2223221‑2223505
                                                                                                                                                       |                                                                                                                                         
aTAGGAATGGTGGCCGAACAGCAGACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCg                                                                                                                                           >  6:291623/1‑149 (MQ=255)
                   agcagACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCCGCAGCTGCCGAacacaatc                                                                                                                        <  1:242643/149‑5 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTg                                                                                                                       <  5:81259/147‑1 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTg                                                                                                                       >  6:81259/1‑147 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCCGCAGCTGCCGAacacaatcatt                                                                                                                     <  6:180442/149‑8 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCCGCAGCTGCCGAacacaatca                                                                                                                       <  1:115936/147‑6 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCCGCAGCTGCCGAacacaatca                                                                                                                       >  2:115936/1‑142 (MQ=255)
                      agACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGAGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCCGCAGCTGCCGAacacaatcatt                                                                                                                     >  5:180442/1‑142 (MQ=255)
                       gACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGt                                                                                                                     >  2:193761/1‑148 (MQ=255)
                        aCGGATCAGTTCAATATTTCTTTCGCTCCGATTGATCAGGATTTCTTCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTa                                                                                                                   <  2:340669/149‑1 (MQ=255)
                           gATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCt                                                                                                                <  3:173765/149‑1 (MQ=255)
                                 gTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCTATGCGGATGGAGGCGCTCGATCACCGTCAC‑‑TTCTTAGAGATAATCAATGATTGTGTTCGGAAGCTGCGGAAAACCTTTGGTTAACTTGATc                                                                                                           >  5:317414/1‑148 (MQ=255)
                                 gTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGAAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATc                                                                                                           >  6:380650/1‑148 (MQ=255)
                                  tCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACATTTGGTTACCTTGATCga                                                                                                         >  4:360974/1‑149 (MQ=255)
                                        tttCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGatgatg                                                                                                    >  3:42649/1‑148 (MQ=255)
                                        cttcgttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGatgatg                                                                                                    >  3:21479/2‑148 (MQ=255)
                                           cgttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGACTTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGaa                                                                                                  <  1:94186/148‑1 (MQ=255)
                                             ttcgCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCAc                                                                                              >  4:91527/1‑149 (MQ=255)
                                                      aTTGATCAGGATTTGTGCAATCGTTTTCCGAAAGCGTTGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTAGAGATAATCAATGATTGTTTTCGGCAGCTGCGGAACACCTTTGTTTACCTTGATCGATGATGAAGCACTCTGaaaaa                                                                                    <  3:360974/149‑1 (MQ=255)
                                                      aTTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGaaaaa                                                                                    <  2:300651/149‑1 (MQ=255)
                                                      aTTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGaaaaa                                                                                    <  8:62070/149‑1 (MQ=255)
                                                           tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCACAGGTCTTCGAGAAAAACAATGATTTTTTACGGCAGATGCGGAACACATTTGTTTACCTTGATCCATAATGAAGAAATATGAAAAAtc                                                                                  >  1:109049/1‑147 (MQ=255)
                                                           tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTGGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCAGTCTGAAAAATGaa                                                                                <  4:42649/148‑1 (MQ=255)
                                                           tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATg                                                                                  >  8:243807/1‑146 (MQ=255)
                                                           tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGaa                                                                                <  7:243807/148‑1 (MQ=255)
                                                           tCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAg                                                                               >  1:112349/1‑149 (MQ=255)
                                                            gatgATTTCTGCAATTGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGTTCCCCGTCTC‑‑TTCTTCGAGATTATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATTGATGATGTAGCACTCTGAAAAATGAAg                                                                               <  2:109049/145‑1 (MQ=255)
                                                             aGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGcc                                                                             >  2:62301/1‑149 (MQ=255)
                                                               gATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGccct                                                                           >  1:242720/1‑149 (MQ=255)
                                                                aTTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGccctc                                                                          <  2:242720/149‑1 (MQ=255)
                                                                     tGCAATCGTTTTCCGAAAGCGATGCGTATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGATATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGt                                                                     <  3:91527/149‑1 (MQ=255)
                                                                       cAATCGTTTTCCGAAAGCGATGCCGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCTAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTa                                                                    <  6:364697/148‑1 (MQ=255)
                                                                        aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGaaaaa                                                                                    <  4:153671/131‑1 (MQ=255)
                                                                        aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGaaaaa                                                                                    >  3:153671/1‑131 (MQ=255)
                                                                        aaTCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAg                                                                   >  3:128637/1‑148 (MQ=255)
                                                                          tCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAg                                                                <  5:291623/149‑1 (MQ=255)
                                                                             ttttCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTc                                                             >  3:307806/1‑149 (MQ=255)
                                                                             ttttCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTc                                                             >  8:347894/1‑149 (MQ=255)
                                                                                   gAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACg                                                         >  6:295433/1‑147 (MQ=255)
                                                                                    aaaGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTa                                                       <  1:62301/149‑1 (MQ=255)
                                                                                    aaaGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACg                                                         >  2:162069/1‑146 (MQ=255)
                                                                                             cGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAgagcggagcg                                             <  3:295199/149‑1 (MQ=255)
                                                                                                  ggAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTtct                                         <  5:295433/148‑1 (MQ=255)
                                                                                                          tCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGaa                                                                                <  3:221567/101‑1 (MQ=255)
                                                                                                          tCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGaa                                                                                >  4:221567/1‑101 (MQ=255)
                                                                                                            gATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTaa                              >  3:236412/1‑149 (MQ=255)
                                                                                                                      aC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATTAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCAc                     <  5:380650/148‑1 (MQ=255)
                                                                                                                      aC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCAc                     <  6:317414/148‑1 (MQ=255)
                                                                                                                                  gaTAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTa          <  4:21479/149‑1 (MQ=255)
                                                                                                                                   aTAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAg         <  4:128637/149‑1 (MQ=255)
                                                                                                                                     aaTCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAgg        >  6:58903/1‑148 (MQ=255)
                                                                                                                                      aTCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTtct                                         <  4:286856/114‑1 (MQ=255)
                                                                                                                                      aTCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTa      <  2:112349/149‑1 (MQ=255)
                                                                                                                                            gATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg  >  2:301552/1‑147 (MQ=255)
                                                                                                                                            gATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg  >  7:225470/1‑147 (MQ=255)
                                                                                                                                            gATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg  >  7:250136/1‑147 (MQ=255)
                                                                                                                                             aaaggtgtTCCGCAGCTGCCGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTCTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAgagg  >  2:288916/5‑147 (MQ=255)
                                                                                                                                               ttGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAgagcgga                                                <  6:362560/99‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                         
ATAGGAATGGTGGCCGAACAGCAGACGGATCAGGTCAATATTTCGTTCGCTCCGATTGATCAGGATTTCTGCAATCGTTTTCCGAAAGCGATGCGGATGGAGGCGCTCGATCCCCGTCAC‑‑TTCTTCGAGATAATCAATGATTGTGTTCGGCAGCTGCGGAACACCTTTGGTTACCTTGATCGATGATGAAGCACTCTGAAAAATGAAGCCCTCCCTGTAGAAGGTCCACGTAGAGCGGAGCGTTCTCAGAGATTTAATGGTTTCACCAACAAATTTAGGTAGAGG  >  NC_002947/2223221‑2223505

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