Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,245 A→G 7.1% intergenic (‑265/+29) PP_5290 ← / ← pyrE transposase/orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,039,2450AG 7.1% 159.9 / 3.0 57intergenic (‑265/+29)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base G (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCG  >  NC_002947/6039099‑6039386
                                                                                                                                                  |                                                                                                                                             
ggggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCATCTgaagaatg                                                                                                                                             >  5:144876/1‑147 (MQ=255)
  ggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTc                                                                                                                                           <  3:245226/149‑1 (MQ=255)
      aTTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCt                                                                                                                                       >  5:101538/1‑149 (MQ=255)
          ctcCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCtttt                                                                                                                                    <  5:303859/148‑1 (MQ=255)
             cTAGGGTCTGTACTGCTAGAACAGTCTGAAGCGAGTTATGCCCCTCCTCCTGCATTCAAGTCCTACCATCCAAGCGGGTTTACCCGCGAAAGCGCCAGATCACGCGTCACATAAAAAGGCGACCCTGTGAAGAAGGTCGCCTTTTTc                                                                                                                                  <  8:137412/147‑1 (MQ=18)
             cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttttt                                                                                                                                   <  8:227830/146‑1 (MQ=255)
             cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                <  7:229755/149‑1 (MQ=255)
             cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                >  5:219336/1‑149 (MQ=255)
                            cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTAAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                  >  4:26493/1‑148 (MQ=255)
                               gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                          >  4:194144/1‑137 (MQ=255)
                               gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                          <  3:194144/137‑1 (MQ=255)
                                aaCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAg                                                                                                               <  3:89706/147‑1 (MQ=255)
                                aaCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAg                                                                                                               >  4:89706/1‑147 (MQ=255)
                                 aCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATc                                                                                                            <  3:26493/149‑1 (MQ=255)
                                        gAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACt                                                                                                     >  2:268799/1‑149 (MQ=255)
                                                   tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAg                                                                                          <  6:219336/149‑1 (MQ=255)
                                                       cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                      >  7:4850/1‑149 (MQ=255)
                                                           ctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCg                                                                                    <  6:101538/147‑1 (MQ=255)
                                                              ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGc                                                                               >  6:210173/1‑149 (MQ=255)
                                                              ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGc                                                                               >  2:133527/1‑149 (MQ=255)
                                                                       aGTACTACCATACAAGCGGGTTTACCCGCGAAAGGTCCAGATCAGGCGTCACAAAAAAAGGCGACCCTCTGAAGAAGGTCTCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATg                                                                      <  8:166511/149‑1 (MQ=255)
                                                                        gTCCTACCATACAGGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGc                                                                     >  4:261259/1‑149 (MQ=255)
                                                                         tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGAt                                                                       <  3:194430/146‑1 (MQ=255)
                                                                         tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGAt                                                                       >  4:194430/1‑146 (MQ=255)
                                                                         tCCTACCATACAAGCGGCTTTACCCGCCAAAGGCCCAGATCAGCCGGCACAAAAAAAGGCGACACTCTGAAGAACGTCGCCTTTTTCATTCCTGGGCCAGCATCAGAACCCCTACTGCGCCCGGTAgcc                                                                                        >  5:27763/1‑127 (MQ=255)
                                                                           cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGcc                                                                                                <  4:131533/119‑1 (MQ=255)
                                                                           cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGcc                                                                                                >  3:131533/1‑119 (MQ=255)
                                                                           cTACCATACAAGCGGGTTGACCCGCGAAATGGCCAGAGCAGGCGTCACAAAAAAAGCCGACCCTCTGAAGCAGGTCGCCTTTTTCATTCCTGGGCCAGGATAAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGTCAGATGCTGt                                                                  <  7:151294/149‑1 (MQ=255)
                                                                            tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGAt                                                                       >  6:58951/1‑143 (MQ=255)
                                                                            tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGAt                                                                       <  5:58951/143‑1 (MQ=255)
                                                                            tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGtt                                                                 >  6:315198/1‑149 (MQ=255)
                                                                              ccATACCAGCGGTTTTACCCGCTAAAGGGCCAGATCAGGCGGCACAAAAAACGGCGACCCTCTGAAGACGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGtt                                                                 <  3:250097/147‑1 (MQ=255)
                                                                                aTACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTg                                                                   <  2:244334/143‑1 (MQ=255)
                                                                                aTACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTg                                                                   >  1:244334/1‑143 (MQ=255)
                                                                                 tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                            >  4:54645/1‑149 (MQ=255)
                                                                                    aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                                  <  2:173327/140‑1 (MQ=255)
                                                                                    aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                                  >  1:173327/1‑140 (MQ=255)
                                                                                    aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                         <  1:268799/149‑1 (MQ=255)
                                                                                    aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                         <  5:315198/149‑1 (MQ=255)
                                                                                    aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAACGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                         <  5:210173/149‑1 (MQ=255)
                                                                                             tACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGt                                                  >  4:250641/1‑147 (MQ=255)
                                                                                                  gcgAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAg                                           >  1:127877/1‑149 (MQ=255)
                                                                                                           ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCa                                                >  4:23854/1‑135 (MQ=255)
                                                                                                           ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCa                                                <  3:23854/135‑1 (MQ=255)
                                                                                                               aTCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCGCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCAcc                              >  5:320184/1‑149 (MQ=255)
                                                                                                                  atgcggcACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                           <  6:144876/147‑1 (MQ=255)
                                                                                                                  aGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                           <  7:254156/149‑1 (MQ=255)
                                                                                                                           aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                            <  1:163632/91‑1 (MQ=255)
                                                                                                                           aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                            >  2:163632/1‑91 (MQ=255)
                                                                                                                            aaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCTGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                 <  3:261259/149‑1 (MQ=255)
                                                                                                                              aaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCg                                                                             <  7:163896/87‑1 (MQ=37)
                                                                                                                              aaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCg                                                                             >  8:163896/1‑87 (MQ=37)
                                                                                                                                 aGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg            >  2:67118/1‑149 (MQ=255)
                                                                                                                                 aGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg            <  7:32569/149‑1 (MQ=255)
                                                                                                                                 aGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg            <  2:127877/149‑1 (MQ=255)
                                                                                                                                  ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc           <  8:4850/149‑1 (MQ=255)
                                                                                                                                   gCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt          <  6:320184/149‑1 (MQ=255)
                                                                                                                                           ctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg  >  4:290299/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCG  >  NC_002947/6039099‑6039386

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: ‑