Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,299,574 C→A 7.1% P114P (CCG→CCT yhbP ← UPF0306 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,299,5740CA 7.1% 191.4 / 5.1 56P114P (CCG→CCTyhbPUPF0306 family protein
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base A (2/2);  total (28/28)
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.17e-01

CCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTC  >  NC_000913/3299428‑3299722
                                                                                                                                                   |                                                                                                                                                    
ccAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCgg                                                                                                                                                    <  4:441790/149‑1 (MQ=255)
      ggACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGa                                                                                                                                              <  1:158228/149‑1 (MQ=255)
          acCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTAGGTGAGTTTGATTTCATCGAGCCGCATGTCCCGCACCGGGGCCGACAGc                                                                                                                                          >  4:37916/1‑149 (MQ=255)
            cGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACAGGGGCCGACAGCAt                                                                                                                                        >  2:272089/1‑149 (MQ=255)
            cGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTATACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACAGGGGCCGACAGCAt                                                                                                                                        >  4:371416/1‑149 (MQ=255)
                   tCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCa                                                                                                                                 >  4:65856/1‑149 (MQ=255)
                           cGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGTAACCGGaaa                                                                                                                         <  4:217411/149‑1 (MQ=255)
                           cGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaaa                                                                                                                         <  4:84523/149‑1 (MQ=255)
                            gTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAg                                                                                                                        <  2:307340/149‑1 (MQ=255)
                             tttAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTTGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGc                                                                                                                       >  3:26009/1‑149 (MQ=255)
                                 aCCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgacg                                                                                                                   >  1:278612/1‑149 (MQ=255)
                                 aCCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgacg                                                                                                                   <  3:480001/149‑1 (MQ=255)
                                 aCCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgacg                                                                                                                   <  3:67649/149‑1 (MQ=255)
                                  ccGCCAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgacga                                                                                                                  >  4:172245/1‑149 (MQ=255)
                                      gAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGt                                                                                                              <  3:126031/149‑1 (MQ=255)
                                          cGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGc                                                                                                          <  4:412568/149‑1 (MQ=255)
                                           gTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGcc                                                                                                         <  3:172245/149‑1 (MQ=255)
                                              aCGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTg                                                                                                      >  4:181181/1‑149 (MQ=255)
                                                   tGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTgcgcgc                                                                                                 >  3:194253/1‑149 (MQ=255)
                                                               gAGTCGCGTAGCCAAATCATTTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCttt                                                                                      >  4:271718/1‑149 (MQ=255)
                                                               gAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTc                                                                                     <  1:474835/149‑1 (MQ=255)
                                                                    gcgTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGTCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTc                                                                                     >  1:520894/1‑144 (MQ=255)
                                                                    gcgTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTc                                                                                     <  2:520894/144‑1 (MQ=255)
                                                                    gcgTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCAc                                                                                >  1:192228/1‑149 (MQ=255)
                                                                           cAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCTAg                                                                         >  3:500580/1‑149 (MQ=255)
                                                                               tCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTg                                                                     <  4:236961/149‑1 (MQ=255)
                                                                                   tttttttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCgg                                                                  >  3:252913/1‑149 (MQ=255)
                                                                                     tttttACCAAAGCCCAGCGTGTTGTTGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACAGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGAt                                                                <  3:371416/149‑1 (MQ=255)
                                                                                     tttttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACAGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGAt                                                                <  1:272089/149‑1 (MQ=255)
                                                                                       tttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGAtct                                                              <  3:524646/149‑1 (MQ=255)
                                                                                         tACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCa                                                            <  3:226591/149‑1 (MQ=255)
                                                                                            cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCa                                                            <  4:77088/146‑1 (MQ=255)
                                                                                            cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCa                                                            >  3:77088/1‑146 (MQ=255)
                                                                                            cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCt                                                         >  3:522333/1‑149 (MQ=255)
                                                                                             aaaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCtt                                                        >  1:521895/1‑149 (MQ=255)
                                                                                              aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCttt                                                       >  3:30190/1‑149 (MQ=255)
                                                                                               aGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCTCACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                                      >  2:9480/1‑149 (MQ=255)
                                                                                                    aGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCAACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTg                                                 <  4:32120/149‑1 (MQ=255)
                                                                                                         gttgtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCacac                                            <  4:26009/149‑1 (MQ=255)
                                                                                                            gtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACcgc                                         <  2:498069/149‑1 (MQ=255)
                                                                                                              cGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGa                                       <  4:30190/149‑1 (MQ=255)
                                                                                                                gTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAAt                                     <  2:556368/149‑1 (MQ=255)
                                                                                                                   aaTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTaa                                  >  4:79841/1‑149 (MQ=255)
                                                                                                                   aaTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTaa                                  <  1:358149/149‑1 (MQ=255)
                                                                                                                    aTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc                                 <  2:370345/149‑1 (MQ=255)
                                                                                                                     tttGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc                                 >  3:12008/1‑148 (MQ=255)
                                                                                                                     tttGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc                                 <  4:12008/148‑1 (MQ=255)
                                                                                                                     tttGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACg                                >  1:476344/1‑149 (MQ=255)
                                                                                                                          tttCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTAcc                           >  1:47407/1‑149 (MQ=255)
                                                                                                                               tCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGtttt                      >  2:472702/1‑149 (MQ=255)
                                                                                                                                  aGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCgg                   >  1:269196/1‑149 (MQ=255)
                                                                                                                                    ccGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCt                 <  4:194253/149‑1 (MQ=255)
                                                                                                                                          ttCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCAcc                                                                               >  1:149575/1‑81 (MQ=255)
                                                                                                                                          ttCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCAcc                                                                               <  2:149575/81‑1 (MQ=255)
                                                                                                                                             ccACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTa        >  1:17090/1‑149 (MQ=255)
                                                                                                                                               acacCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTAcc      >  1:288558/1‑149 (MQ=255)
                                                                                                                                                   cGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTc  >  3:413796/1‑149 (MQ=255)
                                                                                                                                                   cGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTc  >  1:137576/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCGGCTGACCGTTTACCGTTC  >  NC_000913/3299428‑3299722

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