Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,416,458 A→T 100% I216N (ATT→AAT)  USA300HOU_1314 ← oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,416,4580AT100.0% 232.1 / NA 67I216N (ATT→AAT) USA300HOU_1314oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (29/38);  total (29/38)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTCG  >  CP000730/1416325‑1416595
                                                                                                                                     |                                                                                                                                         
aaCACAATCATAGTAATCCCTCTTTTATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                                                                    <  2:259035/141‑1 (MQ=255)
  actactagttttactCCCTCGTTGATTGTGTTCTATATCTTATCAGTCGTTACGAATACGGATCCACAAGGTCTTACGTGAAATGATTAATATTCCCTCGTTCAATGATTTGACCGTTTGTTAAAACATAATTGTACTAAc                                                                                                                                  <  2:302402/127‑1 (MQ=255)
     tttcataaaacttccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAat                                                                                                                               <  2:603634/128‑1 (MQ=255)
     aaTCATAGTAATCCCTCTTTGATTATGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAat                                                                                                                               <  1:367272/141‑1 (MQ=255)
     aaTCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAat                                                                                                                               >  1:755800/1‑141 (MQ=255)
           aGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCa                                                                                                                         >  2:573830/1‑141 (MQ=255)
           aGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCa                                                                                                                         >  1:665103/1‑141 (MQ=255)
            gTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCaa                                                                                                                        >  1:827187/1‑141 (MQ=255)
              aaTCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATa                                                                                                                      >  2:73472/1‑141 (MQ=255)
                    tctTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                                                                                >  2:588269/1‑141 (MQ=255)
                      tttAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACtt                                                                                                              >  1:568635/1‑141 (MQ=255)
                        tAATTGTGTTCTATATTTAATTAGACGTTCAGTACACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTaa                                                                                                            >  2:179905/1‑141 (MQ=255)
                          atggggtTCTATATTGAAGTAGACGTTCAGTACACAGACTCAAACGCTCAAACCTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                                                                          <  1:763383/136‑1 (MQ=255)
                          aTTGTGCTCTATATGTAATTAGACGTTCGGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                                                                          <  2:59519/141‑1 (MQ=255)
                           ttttgtTCTATATTTAATGAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGtt                                                                                                         <  2:454620/138‑1 (MQ=255)
                              tgtTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtg                                                                                                      <  1:349724/141‑1 (MQ=255)
                              tgtTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtg                                                                                                      <  1:134628/141‑1 (MQ=255)
                               ggtCTATATTTAATTAGACGGTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgt                                                                                                     <  2:428427/139‑1 (MQ=255)
                                 tCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTAGTGACAATCTTTAAATACATGACTTAAGTTAtgtgg                                                                                                   >  2:466795/1‑140 (MQ=255)
                                 gcTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCGTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgtgt                                                                                                   >  2:399802/2‑141 (MQ=255)
                                     tatTTAATTATACGTTCAGTATACGGATGCAAATGTTCATACTTGAAATGATTAATATTACATCGTTTAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATAACTTAAGTTAtgtattgaa                                                                                               <  2:108410/141‑7 (MQ=255)
                                           aTTAGACGTTCCGGATACGGACGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                                                                         <  1:457231/141‑1 (MQ=255)
                                           aTTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                                                                         >  2:599379/1‑141 (MQ=255)
                                            ttAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTATTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaata                                                                                        <  1:600858/141‑1 (MQ=255)
                                                aCGTTCAGTATACGGATGCAAATGCTCATCCTTGAAATGATTAATATTACCTCGTTCAACGATTTGACCTTCTTTTAAAACCTAATTGTCCTGACAATATTTCAATACATGAATTACTACAGGGGGGATAATAAATAAGGt                                                                                    >  2:826574/1‑141 (MQ=255)
                                                aCGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGt                                                                                    <  1:657243/141‑1 (MQ=255)
                                                 cGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTCCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGtt                                                                                   <  2:485057/141‑1 (MQ=255)
                                                 cGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCATTACATGACTTAAGTTATGTGTGATAATAAATAATGtt                                                                                   <  1:387007/141‑1 (MQ=255)
                                                 cGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGtt                                                                                   <  1:228624/141‑1 (MQ=255)
                                                  gTTCAATATACGGATTCAAATGCTCATTCGTGAAATTATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGttt                                                                                  <  2:804893/141‑1 (MQ=255)
                                                     cAGTATACGGATGCAAATGCTCATACTTGAAATGCTTAATCTTACATCGTTCTAATATTTGGCCTTATTTTAGATCATAATTGTACTGAGAATATGTCAATACAGGACTTAAGTTATGGGTGATAATAAATAATGTTTGAc                                                                               >  1:414194/1‑141 (MQ=255)
                                                     cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGAc                                                                               >  1:192012/1‑141 (MQ=255)
                                                       gTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                                                                             >  1:635492/1‑141 (MQ=255)
                                                       gTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                                                                             >  2:285434/1‑141 (MQ=255)
                                                         ataCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATg                                                                           >  1:825222/1‑141 (MQ=255)
                                                            cGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGCCAATATTTCAATACATGACTTAAGTTATGTTTGATAATAAATAATGTTTGAcctttttt                                                                        >  2:457071/1‑135 (MQ=255)
                                                                tGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                                                                    <  1:115476/141‑1 (MQ=255)
                                                                  caaactcTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATAc                                                                  <  2:79234/135‑1 (MQ=255)
                                                                  cAAATGCTCATACTTGAAATGATTAATATTACTTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGATAATATTTCAATACATGACTTAAGTTATGTTTAATAATAAATAATGTTTGACCATGTTCTAATAc                                                                  <  2:357955/141‑1 (MQ=255)
                                                                   aaaTGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGTCCATGTTCTAATACa                                                                 >  1:441097/1‑141 (MQ=255)
                                                                       gCTCATACTTGAAATGATTAATATTACCTAGTTCAATGGGTTGTCCTTCTTTTAAAACATACTTTTACAGACACAATTTCAATACTTGGCTTACGAGTTTTGTGATAAGAGTTCATGTTTGACCGTGTTCTAATaaaactt                                                             >  2:284249/1‑138 (MQ=255)
                                                                            ttcTTGAATTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGt                                                        <  2:584720/139‑1 (MQ=255)
                                                                                 gAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGtaataa                                                   >  1:756401/1‑141 (MQ=255)
                                                                                    aTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACAATGTTCTAATACAATATGATGTAAAAAATc                                                >  1:298973/1‑141 (MQ=255)
                                                                                    aTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAACACATTACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATc                                                >  2:3797/1‑141 (MQ=255)
                                                                                      gATGAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCa                                              <  2:274073/141‑1 (MQ=255)
                                                                                       aTTAATCTCACCTCGTTCAATGATTTTACCTTCTTTTACAACATAATTGTACTCACAACATTTCAATACATGACTTATGTTATCTGTGATAATAAATAATGTTTGACCACGTTCTAATAAAATATGCTGTAATACATcctt                                             <  2:72138/141‑3 (MQ=255)
                                                                                         tAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATTTGCTTTAATAAATCCATCa                                           <  1:475620/141‑1 (MQ=255)
                                                                                         tAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATCTGCTTTAATAAATCCATCa                                           <  1:475434/141‑1 (MQ=255)
                                                                                         tAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCa                                           <  1:869341/141‑1 (MQ=255)
                                                                                              ttACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGa                                      <  2:563395/141‑1 (MQ=255)
                                                                                                 ccTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTa                                   <  2:827187/141‑1 (MQ=255)
                                                                                                    cGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCg                                >  1:553189/1‑141 (MQ=255)
                                                                                                     gTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGt                               <  1:107583/141‑1 (MQ=255)
                                                                                                     gTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAAAGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGt                               <  1:813259/141‑1 (MQ=255)
                                                                                                       tCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCCTGTTCTAATACAATATGCTGTAATAACTCCATCACTTGCTTCACGTTc                             >  2:170653/1‑141 (MQ=255)
                                                                                                       tCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCGCTTGATTACCGTTc                             >  1:230578/1‑141 (MQ=255)
                                                                                                         aaTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAGGTTCTAATACAATATGCGGTAAAAAATCCGTCACATGATTACCCTTCaa                           >  2:737172/1‑141 (MQ=255)
                                                                                                            gATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGc                        <  1:341076/141‑1 (MQ=255)
                                                                                                                     ttcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatga               <  1:142522/141‑1 (MQ=255)
                                                                                                                     ttcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatga               <  1:142823/141‑1 (MQ=255)
                                                                                                                     ttcttTTAAAACATAATTGTACTGACAATATTTAAATCCATCACTTAGGTGATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatga               <  2:390556/141‑1 (MQ=255)
                                                                                                                        ttttAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGc            <  1:73472/141‑1 (MQ=255)
                                                                                                                        ttttAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGc            <  2:635492/141‑1 (MQ=255)
                                                                                                                              aaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTgg      >  2:664358/1‑141 (MQ=255)
                                                                                                                              aaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTgg      >  2:450339/1‑141 (MQ=255)
                                                                                                                                  tAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTCg  <  2:568635/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                         
AACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTCG  >  CP000730/1416325‑1416595

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