Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,677,867 C→T K165K (AAG→AAA yphE ← putative ABC transporter ATP‑binding protein YphE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,677,8670CT100.0% 1420.3 / NA 334K165K (AAG→AAAyphEputative ABC transporter ATP‑binding protein YphE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (169/165);  total (169/165)

Only 100 of 340 total aligned reads displayed.

AGTGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCCTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGC  >  NC_000913/2677769‑2677964
                                                                                                  |                                                                                                 
aGTGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCttt                                                                                                 <  1:5130285/101‑1 (MQ=255)
aGTGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCttt                                                                                                 <  2:5130374/101‑1 (MQ=255)
 gTGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTtc                                                                                                >  2:5130929/1‑101 (MQ=255)
  tGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTtca                                                                                               >  1:5130930/1‑101 (MQ=255)
    ccGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTtcatc                                                                                             >  1:5130936/1‑101 (MQ=255)
      gACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCAc                                                                                           >  2:5130939/1‑101 (MQ=255)
       aCATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCAcc                                                                                          >  1:5130945/1‑101 (MQ=255)
        cATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCAccc                                                                                         >  2:5130949/1‑101 (MQ=255)
        cATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCAccc                                                                                         <  2:5130497/101‑1 (MQ=255)
          tCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCcgc                                                                                       >  1:5130952/1‑101 (MQ=255)
          tCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCcgc                                                                                       <  1:5130506/101‑1 (MQ=255)
           cTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCcgcg                                                                                      <  2:5130348/101‑1 (MQ=255)
            tttttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCgcgc                                                                                     <  2:5130705/101‑1 (MQ=255)
            tttttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACAACGCGCGGCTCGCCTTTCATCACCCgcgc                                                                                     >  1:5130956/1‑101 (MQ=255)
             ttttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCa                                                                                    >  2:5130959/1‑101 (MQ=255)
              tttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCAccc                                                                                         >  2:5130965/1‑95 (MQ=255)
              tttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCaa                                                                                   >  1:5130964/1‑101 (MQ=255)
               ttCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAAt                                                                                  <  2:5130552/101‑1 (MQ=255)
                 cACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAAttt                                                                                >  1:5130970/1‑101 (MQ=255)
                 cACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAAttt                                                                                <  1:5130256/101‑1 (MQ=255)
                      cgcTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTccacc                                                                           >  2:5130972/1‑101 (MQ=255)
                      cgcTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTccacc                                                                           >  1:5130975/1‑101 (MQ=255)
                       gcTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTccacca                                                                          <  2:5130562/101‑1 (MQ=255)
                        cTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAg                                                                         >  2:5130978/1‑101 (MQ=255)
                         tGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAgc                                                                        >  1:5130982/1‑101 (MQ=255)
                          gATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAgct                                                                       <  1:5130455/101‑1 (MQ=255)
                           aTCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAgctg                                                                      >  1:5130986/1‑101 (MQ=255)
                             cACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAgctgct                                                                    >  2:5130989/1‑101 (MQ=255)
                              aCCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCtt                                                                   <  1:5130428/101‑1 (MQ=255)
                               ccAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCttt                                                                  >  2:5130993/1‑101 (MQ=255)
                                cAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCtttt                                                                 <  2:5130441/101‑1 (MQ=255)
                                cAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCtttt                                                                 <  1:5130373/101‑1 (MQ=255)
                                cAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCtttt                                                                 <  1:5130477/101‑1 (MQ=255)
                                cAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCtttt                                                                 >  1:5130996/1‑101 (MQ=255)
                                  gTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTgc                                                               >  2:5131002/1‑101 (MQ=255)
                                   ttCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTgcg                                                              <  2:5130289/101‑1 (MQ=255)
                                     cAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAACTGCTTTTGCGcc                                                            >  2:5131006/1‑101 (MQ=255)
                                      aaCTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCg                                                           <  2:5130638/101‑1 (MQ=255)
                                       aCTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCgg                                                          >  1:5131009/1‑101 (MQ=255)
                                         ttCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTgcg                                                              >  2:5131011/1‑95 (MQ=255)
                                          tCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCt                                                       <  1:5130487/101‑1 (MQ=255)
                                            cGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTa                                                     <  1:5130573/101‑1 (MQ=255)
                                             gCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAg                                                    >  2:5131017/1‑101 (MQ=255)
                                             gCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAg                                                    >  1:5131018/1‑101 (MQ=255)
                                              cACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGc                                                   <  1:5130570/101‑1 (MQ=255)
                                              cACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGc                                                   <  1:5130526/101‑1 (MQ=255)
                                                cTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCg                                                           >  1:5131021/1‑91 (MQ=255)
                                                 tCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCgg                                                          <  2:5130098/91‑1 (MQ=255)
                                                  cGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCCCCCGCGCAATTTCCACCAGCTGCTTTTgc                                                               <  2:5130328/85‑1 (MQ=255)
                                                   gCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGtt                                                >  2:5131025/1‑99 (MQ=255)
                                                   gCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGa                                              <  2:5130529/101‑1 (MQ=255)
                                                    cAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTc                                                                               >  2:5131030/1‑67 (MQ=255)
                                                    cAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGaa                                             <  1:5130880/101‑1 (MQ=255)
                                                     aaGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGaaa                                            >  1:5131032/1‑101 (MQ=255)
                                                      aGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAAc                                           <  1:5130401/101‑1 (MQ=255)
                                                       gCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACa                                          >  1:5131037/1‑101 (MQ=255)
                                                        cGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGACTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACaa                                         >  2:5131040/1‑101 (MQ=255)
                                                         gAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCa                                                                                    >  2:5131048/1‑57 (MQ=255)
                                                         gAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAg                                        >  1:5131045/1‑101 (MQ=255)
                                                         gAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAg                                        <  1:5130589/101‑1 (MQ=255)
                                                           aCTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAgtt                                      <  2:5130411/101‑1 (MQ=255)
                                                            cTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAgttg                                     <  1:5130075/101‑1 (MQ=255)
                                                            cTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAgttg                                     >  1:5131051/1‑101 (MQ=255)
                                                            cTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAgttg                                     <  1:5130674/101‑1 (MQ=255)
                                                            cTGGTAGGTTCATCAAGAATGACCACGCGCGGCCCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAgttg                                     <  2:5130458/101‑1 (MQ=255)
                                                              ggTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAAc                                           <  1:5130418/93‑1 (MQ=255)
                                                               gTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTc                                  <  1:5130356/101‑1 (MQ=255)
                                                                  ggTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAgg                               >  2:5131057/1‑101 (MQ=255)
                                                                  ggTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAgg                               >  1:5131056/1‑101 (MQ=255)
                                                                  ggTTCATCAAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAgg                               <  2:5130720/101‑1 (MQ=255)
                                                                     tgatcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACt                            >  1:5131061/3‑101 (MQ=255)
                                                                     tcatcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACt                            <  1:5130471/101‑1 (MQ=255)
                                                                     tcatcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACt                            <  1:5130742/101‑1 (MQ=255)
                                                                     tcatcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACt                            <  1:5130555/101‑1 (MQ=255)
                                                                      catcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTa                           <  2:5130504/101‑1 (MQ=255)
                                                                       atcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTaa                          <  1:5130813/101‑1 (MQ=255)
                                                                        tcaAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAAc                         <  1:5130606/101‑1 (MQ=255)
                                                                         caAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTaa                          <  1:5130847/99‑1 (MQ=255)
                                                                         caAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACg                        <  1:5130652/101‑1 (MQ=255)
                                                                         caAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACg                        <  2:5130597/101‑1 (MQ=255)
                                                                         caAGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACg                        >  2:5131065/1‑101 (MQ=255)
                                                                          aaGAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGt                       <  1:5130623/101‑1 (MQ=255)
                                                                            gAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCa                     <  2:5130830/101‑1 (MQ=255)
                                                                            gAATGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCa                     >  1:5131067/1‑101 (MQ=255)
                                                                              aTGACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAAc                   <  2:5130574/101‑1 (MQ=255)
                                                                               tGACCACGCGCGGCTCGCCTTTCATGACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACg                  >  1:5131073/1‑101 (MQ=255)
                                                                                gACCACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGc                 >  2:5131074/1‑101 (MQ=255)
                                                                                   cACGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCa              >  2:5131080/1‑101 (MQ=255)
                                                                                    aCGCGCGGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAg             >  1:5131081/1‑101 (MQ=255)
                                                                                      gcgcgGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAgcg           <  1:5130436/101‑1 (MQ=255)
                                                                                       cgcgGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAgcgc          >  2:5131084/1‑101 (MQ=255)
                                                                                         cgGCTCGCCTTTCATCCCCCGCGCAATTTCCACCAgctg                                                                      <  2:5130543/39‑1 (MQ=37)
                                                                                         cgGCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCt        >  2:5131089/1‑101 (MQ=255)
                                                                                          ggCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTg       >  1:5131091/1‑101 (MQ=255)
                                                                                           gCTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGATTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGt      >  1:5131095/1‑101 (MQ=255)
                                                                                            cTCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTa     <  2:5130647/101‑1 (MQ=255)
                                                                                             tCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTaa    <  2:5130607/101‑1 (MQ=255)
                                                                                             tCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTaa    <  1:5130703/101‑1 (MQ=255)
                                                                                             tCGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTaa    >  2:5131100/1‑101 (MQ=255)
                                                                                              cGCCTTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGGTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCaa                    >  1:5131103/1‑84 (MQ=255)
                                                                                                  |                                                                                                 
AGTGCCGACATCTTTTTCACCGCGCTGATCACCAGTTCAACTTCCGCACTCGCAAGCGAACTGGTAGGTTCATCAAGAATGACCACGCGCGGCTCGCCCTTCATCACCCGCGCAATTTCCACCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGC  >  NC_000913/2677769‑2677964

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