Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,925 +A coding (35/1431 nt) yaaJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,9221.A98.2% 168.6 / ‑3.7 54coding (38/1431 nt)yaaJputative transporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (26/27);  total (27/27)
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.97e-01

AACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7788‑8056
                                                                                                                                       |                                                                                                                                       
aaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaa                                                                                                                                    >  1:222920/1‑138 (MQ=255)
   cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:227868/139‑5 (MQ=255)
   cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:42347/139‑5 (MQ=255)
   cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  1:1133184/139‑5 (MQ=255)
     ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:1246290/1‑139 (MQ=255)
     ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:1546807/139‑1 (MQ=255)
     ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:509420/139‑1 (MQ=255)
      cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               >  1:751701/1‑139 (MQ=255)
       tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCt                                                                                                                              <  1:7651/139‑1 (MQ=255)
                    gtGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGt                                                                                                                            <  2:1045773/127‑1 (MQ=255)
                    aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGt                                                                                                                            >  1:1045773/1‑128 (MQ=255)
                    aTGCTATTCTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGgag                                                                                                                 >  2:404668/1‑139 (MQ=255)
                     tGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGgaga                                                                                                                >  1:1340108/1‑139 (MQ=255)
                         aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:1128101/1‑119 (MQ=255)
                         aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:1128101/119‑1 (MQ=255)
                          tttttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  1:1513899/1‑112 (MQ=255)
                          tttttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:1513899/116‑5 (MQ=255)
                              tAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTg                                                                                                       <  2:751701/139‑1 (MQ=255)
                               aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                      <  2:1337189/139‑1 (MQ=255)
                                  aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcat                                                                                                   <  1:404668/139‑1 (MQ=255)
                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  2:598834/1‑102 (MQ=255)
                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  1:598834/106‑5 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  1:342633/106‑1 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  2:342633/1‑106 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                               <  2:1340108/139‑1 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  1:481979/138‑1 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                >  2:481979/1‑138 (MQ=255)
                                      ttGCCAAACTGGCGGATGTAGCGAAACTACACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAAAAAAAATCTGgaatgaaa                                                                                               >  2:999417/1‑136 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTAt                                                                                    >  1:142625/1‑139 (MQ=255)
                                                        tAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGt                                                                             >  2:581888/1‑139 (MQ=255)
                                                               aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                            <  2:350062/117‑1 (MQ=255)
                                                               aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                            >  1:350062/1‑117 (MQ=255)
                                                                 tGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:265535/79‑1 (MQ=255)
                                                                 tGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:265535/1‑79 (MQ=255)
                                                                   cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                  >  1:1060722/1‑139 (MQ=255)
                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                >  2:875263/1‑139 (MQ=255)
                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  2:1280612/139‑1 (MQ=255)
                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  1:1087607/139‑1 (MQ=255)
                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  2:142625/139‑1 (MQ=255)
                                                                        aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtg                                                             >  1:267753/1‑139 (MQ=255)
                                                                            ggTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTa                                                         <  1:581888/139‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        >  1:478321/1‑139 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        <  2:1160720/139‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        <  1:875263/139‑1 (MQ=255)
                                                                                  aaaaGTGAACCAACAACCTTCGCGGAACAGCAGGTAAATCATTACCGATCCCCATAAGGACGCTGTTAATGAAGTAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGTTCGCGATGACTTTCGTGTGTAAACGtt                                                   >  1:490928/1‑139 (MQ=255)
                                                                                      gTGAACCAACAACCTGCGGCGAATAGCAGGTAAATCAGTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                               <  1:999417/139‑1 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:1252014/1‑42 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:1252014/42‑1 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                >  1:358702/1‑74 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  2:358702/74‑1 (MQ=255)
                                                                                                           aaGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACg                          <  2:1246290/139‑1 (MQ=255)
                                                                                                                cAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     >  2:1353452/1‑139 (MQ=255)
                                                                                                                   gTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAgg                  <  2:1060722/139‑1 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        >  2:1495907/1‑91 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        <  1:1495907/91‑1 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     <  2:1117373/126‑1 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     >  1:1117373/1‑126 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGt                 >  2:527435/1‑130 (MQ=255)
                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGt                 <  1:527435/130‑1 (MQ=255)
                                                                                                                                atccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                          <  1:121765/66‑1 (MQ=255)
                                                                                                                                atccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                      >  2:1118365/4‑105 (MQ=255)
                                                                                                                                atccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                      <  1:1118365/102‑1 (MQ=255)
                                                                                                                                 tccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTgtgt                                                          <  2:753518/82‑1 (MQ=255)
                                                                                                                                 tccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTgtgt                                                          >  1:753518/3‑84 (MQ=255)
                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                          >  2:121765/1‑69 (MQ=255)
                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc    >  2:689536/1‑139 (MQ=255)
                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc    >  1:897132/1‑139 (MQ=255)
                                                                                                                                  ccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTg  <  2:267753/138‑1 (MQ=255)
                                                                                                                                   cccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                     >  1:1447636/1‑103 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7788‑8056

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