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.C98.4% 222.5 / ‑3.4 63coding (849/1221 nt)sufScysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (33/29);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTG  >  NC_000913/1759537‑1759794
                                                                                                                                           |                                                                                                                        
ttGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                           <  2:1340848/139‑1 (MQ=255)
   tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                        <  1:753134/138‑1 (MQ=255)
   tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                        >  2:753134/1‑138 (MQ=255)
    gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                      <  1:1005756/139‑1 (MQ=255)
            gagagTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTAtt                                                                                                              <  1:52099/139‑1 (MQ=255)
                    gaTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAc                                                                                                      >  1:1040654/1‑139 (MQ=255)
                      tCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCg                                                                                                    >  2:301308/1‑139 (MQ=255)
                       ccGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCgg                                                                                                   >  2:1294133/1‑139 (MQ=255)
                                     aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                           <  2:148205/133‑1 (MQ=255)
                                     aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                           >  1:148205/1‑133 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  2:939592/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  1:1294133/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  1:301308/139‑1 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                              >  2:628730/1‑139 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTAAAACCGCCATGGTGCtt                                                                              >  1:952859/1‑139 (MQ=255)
                                               gcgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTgg                                                                           <  2:952859/139‑1 (MQ=255)
                                                    tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAc                                                                                                      >  1:738776/1‑107 (MQ=255)
                                                    tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAc                                                                                                      <  2:738776/107‑1 (MQ=255)
                                                     aaTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCaaa                                                                                            >  1:712465/1‑116 (MQ=255)
                                                     aaTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCaaa                                                                                            <  2:712465/116‑1 (MQ=255)
                                                           tCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                           <  2:1188296/111‑1 (MQ=255)
                                                           tCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAc                                                                                           >  1:1188296/1‑111 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCg                                                                                         <  2:498599/112‑1 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCg                                                                                         >  1:498599/1‑112 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                              <  2:125196/139‑1 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                              >  2:954379/1‑139 (MQ=255)
                                                              gATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCtt                                                            >  1:1103401/1‑139 (MQ=255)
                                                                 tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                   >  2:236794/1‑129 (MQ=255)
                                                                 tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                   <  1:236794/129‑1 (MQ=255)
                                                                  cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:541271/139‑1 (MQ=255)
                                                                  cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:954379/139‑1 (MQ=255)
                                                                          aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGa                                                >  2:415820/1‑139 (MQ=255)
                                                                              ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                        <  2:1322580/63‑1 (MQ=255)
                                                                              ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcc                                                                                                                        >  1:1322580/1‑63 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:914685/126‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        >  2:914685/1‑126 (MQ=255)
                                                                                 tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCTGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         <  1:599576/139‑1 (MQ=255)
                                                                                   tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGa                                       >  2:663608/1‑139 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                  <  2:1200518/139‑1 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                  >  1:1200518/1‑139 (MQ=255)
                                                                                          aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc                                                                                                                     >  1:676593/1‑54 (MQ=255)
                                                                                          aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc                                                                                                                     <  2:676593/54‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAc                                                                                                      <  1:158382/66‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAc                                                                                                      >  2:158382/1‑66 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         >  2:709752/1‑127 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         <  1:709752/127‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                             >  1:603351/1‑139 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                             >  1:211091/1‑139 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                      <  1:867371/48‑1 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGcccc                                                                                                                      >  2:867371/1‑48 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        >  1:174489/1‑139 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        >  1:391184/1‑139 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        <  1:415820/139‑1 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                            <  1:1024114/129‑1 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                            >  2:1024114/1‑129 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcc                          <  1:715192/125‑1 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcc                          >  2:715192/1‑125 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcggcg            <  2:211091/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc          <  2:1103401/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc          <  2:391184/139‑1 (MQ=255)
                                                                                                                  cAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAt        >  2:44166/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:1284633/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:790906/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  1:718723/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCCGGTATTCGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:605460/1‑139 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTcc    <  2:1040654/139‑1 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                  >  2:1342282/1‑75 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                  <  1:1342282/75‑1 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTg  <  2:603351/139‑1 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  2:557444/60‑1 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               >  1:557444/1‑60 (MQ=255)
                                                                                                                                    aaTGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                          >  2:480185/1‑55 (MQ=255)
                                                                                                                                    aaTGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                          <  1:480185/55‑1 (MQ=255)
                                                                                                                                           |                                                                                                                        
TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTG  >  NC_000913/1759537‑1759794

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: