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

TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTG  >  NC_000913/1759537‑1759794
                                                                                                                                           |                                                                                                                        
ttGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                           <  2:1340848/139‑1 (MQ=255)
   tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                        <  1:753134/138‑1 (MQ=255)
   tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                        >  2:753134/1‑138 (MQ=255)
    gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                      <  1:1005756/139‑1 (MQ=255)
            gagagTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTAtt                                                                                                              <  1:52099/139‑1 (MQ=255)
                    gaTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc                                                                                                      >  1:1040654/1‑139 (MQ=255)
                      tCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCg                                                                                                    >  2:301308/1‑139 (MQ=255)
                       ccGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCgg                                                                                                   >  2:1294133/1‑139 (MQ=255)
                                     aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAAc                                                                                           <  2:148205/133‑1 (MQ=255)
                                     aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAAc                                                                                           >  1:148205/1‑133 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  2:939592/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  1:1294133/139‑1 (MQ=255)
                                           gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  1:301308/139‑1 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                              >  2:628730/1‑139 (MQ=255)
                                            aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTAAAACCGCCATGGTGCtt                                                                              >  1:952859/1‑139 (MQ=255)
                                               gcgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTgg                                                                           <  2:952859/139‑1 (MQ=255)
                                                    tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc                                                                                                      >  1:738776/1‑107 (MQ=255)
                                                    tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc                                                                                                      <  2:738776/107‑1 (MQ=255)
                                                     aaTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCaaa                                                                                            >  1:712465/1‑116 (MQ=255)
                                                     aaTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCaaa                                                                                            <  2:712465/116‑1 (MQ=255)
                                                           tCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAAc                                                                                           <  2:1188296/111‑1 (MQ=255)
                                                           tCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAAc                                                                                           >  1:1188296/1‑111 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCg                                                                                         <  2:498599/112‑1 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCg                                                                                         >  1:498599/1‑112 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                              <  2:125196/139‑1 (MQ=255)
                                                            cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                              >  2:954379/1‑139 (MQ=255)
                                                              gATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCtt                                                            >  1:1103401/1‑139 (MQ=255)
                                                                 tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                   >  2:236794/1‑129 (MQ=255)
                                                                 tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                   <  1:236794/129‑1 (MQ=255)
                                                                  cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:541271/139‑1 (MQ=255)
                                                                  cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:954379/139‑1 (MQ=255)
                                                                          aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGa                                                >  2:415820/1‑139 (MQ=255)
                                                                              ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                        <  2:1322580/63‑1 (MQ=255)
                                                                              ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                        >  1:1322580/1‑63 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        <  1:914685/126‑1 (MQ=255)
                                                                               gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                        >  2:914685/1‑126 (MQ=255)
                                                                                 tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCCTGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         <  1:599576/139‑1 (MQ=255)
                                                                                   tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGa                                       >  2:663608/1‑139 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                  <  2:1200518/139‑1 (MQ=255)
                                                                                        ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                  >  1:1200518/1‑139 (MQ=255)
                                                                                          aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccc                                                                                                                     >  1:676593/1‑54 (MQ=255)
                                                                                          aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccc                                                                                                                     <  2:676593/54‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc                                                                                                      <  1:158382/66‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc                                                                                                      >  2:158382/1‑66 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         >  2:709752/1‑127 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                         <  1:709752/127‑1 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                             >  1:603351/1‑139 (MQ=255)
                                                                                             ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                             >  1:211091/1‑139 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                      <  1:867371/48‑1 (MQ=255)
                                                                                               ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                      >  2:867371/1‑48 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        >  1:174489/1‑139 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        >  1:391184/1‑139 (MQ=255)
                                                                                                  gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                        <  1:415820/139‑1 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                            <  1:1024114/129‑1 (MQ=255)
                                                                                                        aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                            >  2:1024114/1‑129 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcc                          <  1:715192/125‑1 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcc                          >  2:715192/1‑125 (MQ=255)
                                                                                                              aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcggcg            <  2:211091/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc          <  2:1103401/139‑1 (MQ=255)
                                                                                                                tCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc          <  2:391184/139‑1 (MQ=255)
                                                                                                                  cAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAt        >  2:44166/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:1284633/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:790906/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  1:718723/1‑139 (MQ=255)
                                                                                                                     cgccgcgccAAGACCAATGATGCCCCCC‑GGTATTCGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc     >  2:605460/1‑139 (MQ=255)
                                                                                                                      gccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTcc    <  2:1040654/139‑1 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                  >  2:1342282/1‑75 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                  <  1:1342282/75‑1 (MQ=255)
                                                                                                                        cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTg  <  2:603351/139‑1 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               <  2:557444/60‑1 (MQ=255)
                                                                                                                          cgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                               >  1:557444/1‑60 (MQ=255)
                                                                                                                                    aaTGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                          >  2:480185/1‑55 (MQ=255)
                                                                                                                                    aaTGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                          <  1:480185/55‑1 (MQ=255)
                                                                                                                                           |                                                                                                                        
TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTG  >  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: ‑