Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,621,1180AC45.5% 74.6 / 12.6 56V7G (GTG→GGG) eamAcysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (21/9);  new base C (0/25);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 1.89e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620968‑1621268
                                                                                                                                                      |                                                                                                                                                      
ggTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAaca                                                                                                                                                        <  2:573935/151‑1 (MQ=255)
  tACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcc                                                                                                                                                      >  1:1105345/1‑151 (MQ=255)
        cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                                <  1:839697/151‑1 (MQ=255)
         ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                               <  2:557946/151‑1 (MQ=255)
          gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                              >  1:231910/1‑151 (MQ=255)
            cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                            >  2:1110750/1‑151 (MQ=255)
            cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                            <  2:568454/151‑1 (MQ=255)
             gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                                           <  1:244639/151‑1 (MQ=255)
                          gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCg                                                                                                                              >  2:309323/1‑151 (MQ=255)
                                 ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                        >  1:238350/1‑150 (MQ=255)
                                    aaaGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAAtt                                                                                                                    <  2:672662/151‑1 (MQ=255)
                                        cGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                                >  2:384566/1‑151 (MQ=255)
                                          aCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                              >  2:220219/1‑151 (MQ=255)
                                            cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCttt                                                                                                            >  1:139977/1‑151 (MQ=255)
                                             aGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                           >  1:100055/1‑151 (MQ=255)
                                                aTAAAGGGCAACCCGGCCCGTATCAGCGGGTCCATTTTTTAAACCCCCACTTTGATCACCAAAAAATTTAGCCCCCATACATCCCCTCCCATTAGCACCACCGCCCCATCTTTTCGCGACATTTTACCGCCTCTGAATTTCATCtt                                                                                                             <  1:952819/146‑1 (MQ=255)
                                                aTAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                           <  2:260995/148‑1 (MQ=255)
                                                  aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                        >  1:144093/1‑149 (MQ=255)
                                                  aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                        >  2:357128/1‑149 (MQ=255)
                                                  aaaGCCCAAACCGGCAAGCATCAGCGGTGGCATTTTATGAAGCCCCATTTTGATGACCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCGCCAACACCCCATTTTTTCGCCACATTCTACCGCCTCTGAATTTCATCTTTTGTAAg                                                                                                      <  1:372187/151‑1 (MQ=255)
                                                      cgcAACCCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                                   <  1:1048077/150‑1 (MQ=255)
                                                         aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                                 <  2:267516/146‑1 (MQ=255)
                                                          aaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAAc                                                                                              >  1:402385/1‑151 (MQ=255)
                                                           aCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                             >  2:716915/1‑151 (MQ=255)
                                                              gcccAGCATCGGCGGTGGCATGTTATGAACCCCCTCTTTGATGACCACAAAATTTAGCCCCCATCCGCCCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                            <  2:384767/147‑1 (MQ=255)
                                                                ccAGCATCAGCGGTGCCATGTTATGACGCCCCCTTTTGATGACCACAAAATTTAGCCCCCAACCGCCCACTACCAGTAGCGCCACCCCCCCTTCTTTTCGCACCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                          <  1:737453/149‑1 (MQ=255)
                                                                  aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                         <  2:1078357/148‑1 (MQ=255)
                                                                     aTCAGCGGTGGCATGTTATGACGCCCCTTTTTGATGACCAAAAAATTTCCCCCCCATACGCCCACTCCCAGTAGCGCAAACCCCCTATTTTTGCGCGACATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                       <  1:334714/147‑1 (MQ=255)
                                                                       cAGCGGTGGCATGTTATGAAGCCCCTCTTTGATGACCACAAAATTTACCCCCCATCCGCCCACTCCCATTAGCCCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGa                                                                                      <  2:156371/146‑1 (MQ=255)
                                                                              ggCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                          >  2:1025881/1‑151 (MQ=255)
                                                                              ggCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                          >  2:608428/1‑151 (MQ=255)
                                                                              ggCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                          <  2:1035986/151‑1 (MQ=255)
                                                                                 attttATGAACCCCCACTTGGATGCCCACAAAATTTGCCCCCCATACCACCACTACCATTACCGCAAACCCCCCATCTTTTCGCTACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                        <  1:1048699/147‑1 (MQ=255)
                                                                                    ttATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                     >  1:172727/1‑150 (MQ=255)
                                                                                      aTGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  <  2:1014010/151‑1 (MQ=255)
                                                                                         aaGCCCCATTTTGATGACCAAAAAATTTACCCCCCATCCGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                   <  1:956008/147‑1 (MQ=255)
                                                                                            ccccTCTTGAATGCCCACAAAATTTACCCCCCACCCGACCACTGCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                             <  1:173223/150‑1 (MQ=255)
                                                                                            ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                             <  1:798408/150‑1 (MQ=255)
                                                                                            ccccACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                             <  1:62749/150‑1 (MQ=255)
                                                                                             cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCAATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           >  1:625774/1‑151 (MQ=255)
                                                                                              ccACTTTGATGACCAAAAAATTTACCCCCCATACGCCCACTCCCAGTAGCACCAACCCCCCACTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                                          <  1:1001868/151‑1 (MQ=255)
                                                                                               cACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                         >  1:623277/1‑151 (MQ=255)
                                                                                                   ttGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTc                                                     >  1:818573/1‑151 (MQ=255)
                                                                                                       tgaCCACAAAATTTACCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                                  <  2:445453/150‑1 (MQ=255)
                                                                                                            acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                            <  2:1076273/151‑1 (MQ=255)
                                                                                                             caAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                           >  2:151226/1‑151 (MQ=255)
                                                                                                               aaaTTTAGCCCCCATACGACCCCTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                         <  2:745560/151‑1 (MQ=255)
                                                                                                                      gCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   >  2:1096958/1‑150 (MQ=255)
                                                                                                                       cccccACCCGACCACTCCCAGTAGCGCCACCCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                    <  1:452521/148‑1 (MQ=255)
                                                                                                                          ccATACGCCCACTCCCAGTAGCGCCACCACCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                              <  1:671634/151‑1 (MQ=255)
                                                                                                                                     cTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt                    <  2:1056243/150‑1 (MQ=255)
                                                                                                                                        ccAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg                  <  1:319150/149‑1 (MQ=255)
                                                                                                                                           gTAGCGCCACCCCCCCATCTTTTCGCAATATTCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtata             <  1:1094672/151‑1 (MQ=255)
                                                                                                                                              gcgcCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt          <  1:129838/151‑1 (MQ=255)
                                                                                                                                              gcgcCAACACCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt           <  2:224089/150‑1 (MQ=255)
                                                                                                                                               cgcCACCCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGACATTAATGCGCAAATAAGTCGCTATACtt          <  1:574595/150‑1 (MQ=255)
                                                                                                                                               cgcCAACCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt          <  1:309312/150‑1 (MQ=255)
                                                                                                                                                    acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtt    >  1:458483/1‑151 (MQ=255)
                                                                                                                                                      acCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt  >  2:870919/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620968‑1621268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: