Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 852,869 C→T *43* (TAG→TAA)  mntS ← Mn(2)‑response protein, MntR‑repressed

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913852,8690CT98.6% 214.9 / ‑4.3 70*43* (TAG→TAA) mntSMn(2)‑response protein, MntR‑repressed
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/35);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

CGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAA  >  NC_000913/852734‑852999
                                                                                                                                       |                                                                                                                                  
cgccgTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTAt                                                                                                                                 >  2:738674/1‑139 (MQ=255)
      tCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCg                                                                                                                           <  1:508723/139‑1 (MQ=255)
       cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAGAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCgg                                                                                                                          <  2:360310/139‑1 (MQ=255)
       cAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCgg                                                                                                                          <  2:542599/139‑1 (MQ=255)
               ctactGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTa                                                                                                                  >  1:1223173/1‑139 (MQ=255)
                     aCGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgc                                                                                                            >  1:53664/1‑139 (MQ=255)
                      cGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgct                                                                                                           <  2:351467/139‑1 (MQ=255)
                          cAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCgg                                                                                                       >  1:635332/1‑139 (MQ=255)
                           aaaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      <  2:635332/139‑1 (MQ=255)
                           aaaTCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      <  2:1247986/139‑1 (MQ=255)
                               cTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTttgt                                                                                                  >  2:314291/1‑139 (MQ=255)
                                ttttGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTttgtt                                                                                                 <  1:335954/139‑1 (MQ=255)
                                   tGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgtt                                                                                              <  2:53664/139‑1 (MQ=255)
                                     cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgctg                                                                                                          >  2:1177985/1‑125 (MQ=255)
                                     cGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgctg                                                                                                          <  1:1177985/125‑1 (MQ=255)
                                              ttGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                   <  1:314291/139‑1 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgct                                                                                                           <  2:1024028/109‑1 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgct                                                                                                           >  1:1024028/1‑109 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCg                                                                                                        >  2:810631/1‑112 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCg                                                                                                        <  1:810631/112‑1 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCacaaca                                                                             >  1:1023977/1‑139 (MQ=255)
                                                    aTCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCacaaca                                                                             >  2:584607/1‑139 (MQ=255)
                                                      ccGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATa                                                                           <  1:968396/139‑1 (MQ=255)
                                                         aCAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGa                                                                                                                         <  1:52601/90‑1 (MQ=255)
                                                         aCAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGa                                                                                                                         >  2:52601/1‑90 (MQ=255)
                                                         aCAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATagag                                                                        >  1:432873/1‑139 (MQ=255)
                                                              aaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCg                                                                                                                           >  2:1376156/1‑83 (MQ=255)
                                                              aaGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCg                                                                                                                           <  1:1376156/83‑1 (MQ=255)
                                                                   gggTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTa                                                              >  1:347719/1‑139 (MQ=255)
                                                                        tttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCTTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGt                                                         <  2:1023977/139‑1 (MQ=255)
                                                                        tttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGt                                                         <  2:297363/139‑1 (MQ=255)
                                                                        tttCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGt                                                         <  1:738674/139‑1 (MQ=255)
                                                                              aCCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTa                                                   >  2:1087876/1‑139 (MQ=255)
                                                                                cAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaaa                                                 <  2:1386783/139‑1 (MQ=255)
                                                                                  aaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGt                                                         <  1:270650/129‑1 (MQ=255)
                                                                                  aaaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGt                                                         >  2:270650/1‑129 (MQ=255)
                                                                                   aaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      <  2:623265/83‑1 (MQ=255)
                                                                                   aaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      >  1:623265/1‑83 (MQ=255)
                                                                                   aaCGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAggg                                              >  2:741688/1‑139 (MQ=255)
                                                                                          tATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGa                                       <  1:1127714/139‑1 (MQ=255)
                                                                                                tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgtt                                                                                              <  1:1096189/78‑1 (MQ=255)
                                                                                                tGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTttgttgtt                                                                                              >  2:1096189/1‑78 (MQ=255)
                                                                                                 gCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAaca                                <  1:1087876/139‑1 (MQ=255)
                                                                                                 gCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAaca                                <  2:347719/139‑1 (MQ=255)
                                                                                                  cATGTAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                   >  2:67865/1‑87 (MQ=255)
                                                                                                  cATGTAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGc                                                                                   <  1:67865/87‑1 (MQ=255)
                                                                                                  cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGa                                                                                                                         >  2:787078/1‑49 (MQ=255)
                                                                                                  cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGa                                                                                                                         <  1:787078/49‑1 (MQ=255)
                                                                                                  cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaaa                                                 <  2:132099/121‑1 (MQ=255)
                                                                                                  cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTaaa                                                 >  1:132099/1‑121 (MQ=255)
                                                                                                  cATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAacac                               >  2:222111/1‑139 (MQ=255)
                                                                                                    tGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACAcgc                             >  2:580722/1‑139 (MQ=255)
                                                                                                      aaaGCAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACgcgc                           >  2:1377523/1‑139 (MQ=255)
                                                                                                          cAGGAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATac                       <  1:222111/139‑1 (MQ=255)
                                                                                                            ggAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGca                                                                                  <  2:1278153/78‑1 (MQ=255)
                                                                                                            ggAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGca                                                                                  >  1:1278153/1‑78 (MQ=255)
                                                                                                            ggAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgaga                                           <  2:1250483/117‑1 (MQ=255)
                                                                                                            ggAGTGAAGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGgaga                                           >  1:1250483/1‑117 (MQ=255)
                                                                                                                  aaGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                >  2:783642/1‑74 (MQ=255)
                                                                                                                  aaGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGcaca                                                                                <  1:783642/74‑1 (MQ=255)
                                                                                                                  aaGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTg               >  1:316881/1‑139 (MQ=255)
                                                                                                                   aGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      <  2:111541/51‑1 (MQ=255)
                                                                                                                   aGCGGCGTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGt                                                                                                      >  1:111541/1‑51 (MQ=255)
                                                                                                                    gcggcgTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgct                                                                                                           >  2:1379874/1‑45 (MQ=255)
                                                                                                                    gcggcgTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTtgct                                                                                                           <  1:1379874/45‑1 (MQ=255)
                                                                                                                     cggcgTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAAc            <  1:741688/139‑1 (MQ=255)
                                                                                                                       gcgTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGca                                                                                  >  1:797446/1‑67 (MQ=255)
                                                                                                                       gcgTACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGca                                                                                  <  2:797446/67‑1 (MQ=255)
                                                                                                                          tACCGCGACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACtcatt       >  1:818405/1‑139 (MQ=255)
                                                                                                                               cgACGCCGTTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATaa  >  1:509655/1‑139 (MQ=255)
                                                                                                                                  cgccgtTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGa                                                                                            <  2:744770/46‑1 (MQ=255)
                                                                                                                                  cgccgtTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGa                                                                                            >  1:744770/1‑46 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CGCCGTTCAGGTCAACTACTGACGAACAAATCTTTTGCGAAGGGGATTGCAAATCCGACAAAAAGTGGGGTTTTTCAGACCAAAACGCGTTATGCCTGCATGAAAGCAGGAGTGAAGCGGCGTACCGCGACGCCGCTATTTATCGGAAGGTTTATCTTGCTGCGGTTTGTTGTTGACCATATCGCACAACATAGAGAGCAGCATTAACCGTACTTTAAAGGGAGAATGACTAAACACGCGCATACACCTCTTGAACTCATTCATAA  >  NC_000913/852734‑852999

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