Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,759,680 (C)5→6 coding (844/1221 nt) sufS ← cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,759,6751.C100.0% 240.1 / NA 68G283G (GGC→GGGsufScysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (33/35);  total (33/35)

TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAAC  >  NC_000913/1759537‑1759798
                                                                                                                                           |                                                                                                                           
ttGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                              <  2:1567824/139‑1 (MQ=255)
   tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccc                                                                                                                          >  2:1190034/1‑139 (MQ=255)
    gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                         <  1:255268/139‑1 (MQ=255)
     ggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc                                                                                                                        >  2:760777/1‑139 (MQ=255)
          tAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccccggta                                                                                                                   <  2:414039/139‑6 (MQ=255)
              gagTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTgg                                                                                                                <  2:387972/139‑1 (MQ=255)
               agTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTg                                                                                                                 >  2:24147/1‑137 (MQ=255)
               agTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTg                                                                                                                 <  1:24147/137‑1 (MQ=255)
                         ggTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCt                                                                                                     >  2:1483557/1‑139 (MQ=255)
                         ggTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                         <  1:1128186/118‑1 (MQ=255)
                         ggTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                         >  2:1128186/1‑118 (MQ=255)
                              aGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGTCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCaaa                                                                                                <  1:1190034/139‑1 (MQ=255)
                                    cAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                          >  2:849292/1‑139 (MQ=255)
                                          tgATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                                                                                    >  1:1315059/1‑139 (MQ=255)
                                          tgATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                                                                                    >  1:1534698/1‑139 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  2:1534698/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  2:730787/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  1:1106761/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  1:1044949/139‑1 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                                  >  1:771508/1‑139 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                                  >  2:540191/1‑139 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                                  <  2:1490312/139‑1 (MQ=255)
                                             tAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                 <  1:760777/139‑1 (MQ=255)
                                                 gcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc                                                                             <  1:141186/139‑1 (MQ=255)
                                                        gCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCCCCATGGTGCTTTGGTCCAGGTAg                                                                      <  1:1483557/139‑1 (MQ=255)
                                                               atccTGTTCATGCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCTGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   >  1:934065/4‑119 (MQ=255)
                                                               aTTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  2:934065/119‑1 (MQ=255)
                                                                 tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCAc                                                             >  2:472427/1‑139 (MQ=255)
                                                                   tGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                               >  2:1594724/1‑103 (MQ=255)
                                                                   tGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                               <  1:1594724/103‑1 (MQ=255)
                                                                          aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                                  >  2:1107407/1‑125 (MQ=255)
                                                                          aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                                  <  1:1107407/125‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                      >  1:1453396/1‑116 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                      <  2:1453396/116‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCtt                                                                <  1:413025/122‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCtt                                                                >  2:413025/1‑122 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCt                                                                 <  2:467818/121‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCt                                                                 >  1:467818/1‑121 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTg                                               >  1:725959/1‑139 (MQ=255)
                                                                                   tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGa                                           >  2:60010/1‑139 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGATACATACTCCAGCGCCGCGCCAAGACCAATGATGccccccggta                                                                                                                   <  2:1077842/61‑6 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccccggta                                                                                                                   >  1:1077842/1‑56 (MQ=255)
                                                                                          aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGa                                    >  1:371377/1‑139 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  2:79407/89‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  1:826414/89‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   >  2:558061/1‑89 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   >  2:826414/1‑89 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   >  1:79407/1‑89 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                   <  1:558061/89‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                             >  2:1479097/1‑127 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                             <  1:1479097/127‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                 >  1:575653/1‑139 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                           >  1:908181/1‑48 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                           <  2:908181/48‑1 (MQ=255)
                                                                                              cccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                                <  2:575653/139‑1 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                          <  1:326433/80‑1 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                          >  2:326433/1‑80 (MQ=255)
                                                                                                 aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc                                                                             <  2:406943/91‑1 (MQ=255)
                                                                                                 aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc                                                                             >  1:406943/1‑91 (MQ=255)
                                                                                                 aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGccc                             >  2:255973/1‑139 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                            >  1:1156857/1‑139 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                              <  1:203736/83‑1 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                              >  2:203736/1‑83 (MQ=255)
                                                                                                               cTCCAGCGCCGCGCCAAGACCAGTGATGCCCCCCGGTATTGGGTGTACCGGCTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc              <  1:456534/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc              <  2:771508/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc              <  1:60010/139‑1 (MQ=255)
                                                                                                                   aGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcg                   <  1:233226/131‑1 (MQ=255)
                                                                                                                   aGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcg                   >  2:233226/1‑131 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTcc                                                                            <  2:493128/71‑1 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTcc                                                                            >  1:493128/1‑71 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTcc                       <  1:346549/124‑1 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTcc                       >  2:346549/1‑124 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc              <  2:954582/131‑1 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc              >  1:954582/1‑131 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                            <  2:1238128/115‑1 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                            >  1:1238128/1‑115 (MQ=255)
                                                                                                                            ccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGcaac  <  1:472427/139‑1 (MQ=255)
                                                                                                                                gACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGt                                                <  2:359518/89‑1 (MQ=255)
                                                                                                                                gACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGt                                                >  1:359518/1‑89 (MQ=255)
                                                                                                                                           |                                                                                                                           
TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAAC  >  NC_000913/1759537‑1759798

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