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.2% 192.9 / ‑3.6 56coding (849/1221 nt)sufScysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (28/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑C‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTT  >  NC_000913/1759541‑1759808
                                                                                                                                       |                                                                                                                                      
gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑ccc                                                                                                                                   <  1:1295611/139‑1 (MQ=255)
      tAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccccggta                                                                                                                             <  1:1669843/139‑6 (MQ=255)
             tGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑c                                                                                                                                     <  1:450069/124‑1 (MQ=255)
             tGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑c                                                                                                                                     >  2:450069/1‑124 (MQ=255)
              gagaTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGgtgt                                                                                                                      >  2:693269/1‑139 (MQ=255)
                            aTTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGc                                                                                                                >  1:747792/1‑131 (MQ=255)
                            aTTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGc                                                                                                                <  2:747792/131‑1 (MQ=255)
                                cAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                                    <  1:1029742/139‑1 (MQ=255)
                                      tgATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                                                                                              >  2:558006/1‑139 (MQ=255)
                                       gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             <  1:558006/139‑1 (MQ=255)
                                       gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             <  2:554454/139‑1 (MQ=255)
                                              cATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTcc                                                                                      >  2:1302094/1‑139 (MQ=255)
                                              cATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTcc                                                                                      >  1:189834/1‑139 (MQ=255)
                                                tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                         >  2:1471804/1‑87 (MQ=255)
                                                tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                         <  1:1471804/87‑1 (MQ=255)
                                                                 ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATTCTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAg                                                                   >  1:221152/1‑139 (MQ=255)
                                                                 ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccc                                                                                                                                  >  2:421540/1‑75 (MQ=255)
                                                                 ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccc                                                                                                                                  <  1:421540/75‑1 (MQ=255)
                                                                    aTACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTAc                                                                                                                    >  2:1685066/1‑87 (MQ=255)
                                                                    aTACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CGCCCCAGTATTGGGTGTAc                                                                                                                    <  1:1685066/87‑1 (MQ=255)
                                                                         cGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccccggta                                                                                                                             <  1:1183139/72‑6 (MQ=255)
                                                                         cGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccccggta                                                                                                                             >  2:1183139/1‑67 (MQ=255)
                                                                           gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                                        >  2:504588/1‑108 (MQ=255)
                                                                           gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                                        <  1:504588/108‑1 (MQ=255)
                                                                           gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      <  2:1010429/126‑1 (MQ=255)
                                                                           gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      >  1:1010429/1‑126 (MQ=255)
                                                                             tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                       <  1:693269/139‑1 (MQ=255)
                                                                               tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                       <  2:1266325/137‑1 (MQ=255)
                                                                               tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                       >  1:1266325/1‑137 (MQ=255)
                                                                                    ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                                        >  2:84416/1‑99 (MQ=255)
                                                                                    ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                                        <  1:84416/99‑1 (MQ=255)
                                                                                    ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑ccc                                                                                                                                   >  1:1258895/1‑55 (MQ=255)
                                                                                    ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑ccc                                                                                                                                   <  2:1258895/55‑1 (MQ=255)
                                                                                      aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑ccc                                                                                                                                   <  1:231571/53‑1 (MQ=255)
                                                                                      aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑ccc                                                                                                                                   >  2:231571/1‑53 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCACTGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             >  1:306525/1‑89 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                           >  2:434426/1‑91 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                           <  1:434426/91‑1 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             <  2:306525/89‑1 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      >  1:1244337/1‑112 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      <  2:1244337/112‑1 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTc                                                                     >  1:453099/1‑113 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTc                                                                     <  2:453099/113‑1 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                           >  2:1464847/1‑139 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAc                                            <  1:101414/138‑1 (MQ=255)
                                                                                         ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAc                                            >  2:101414/1‑138 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                         >  2:1154643/1‑44 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                         <  1:1154643/44‑1 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             >  2:91802/1‑87 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                             <  1:91802/87‑1 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                                <  2:614028/100‑1 (MQ=255)
                                                                                           ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                                >  1:614028/1‑100 (MQ=255)
                                                                                             aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGccc                                       >  2:1557829/1‑139 (MQ=255)
                                                                                              gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAc                                            <  1:1192604/133‑1 (MQ=255)
                                                                                              gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAc                                            >  2:1192604/1‑133 (MQ=255)
                                                                                              gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                                          <  2:654267/135‑1 (MQ=255)
                                                                                              gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg                                          >  1:654267/1‑135 (MQ=255)
                                                                                                   gAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                                <  2:1075930/92‑1 (MQ=255)
                                                                                                   gAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                                >  1:1075930/1‑92 (MQ=255)
                                                                                                   gAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccc                                                                                                                                  >  1:183453/1‑41 (MQ=255)
                                                                                                   gAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑c‑cccc                                                                                                                                  <  2:183453/41‑1 (MQ=255)
                                                                                                    aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGc                                         <  1:1374498/130‑1 (MQ=255)
                                                                                                    aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGc                                         >  2:1374498/1‑130 (MQ=255)
                                                                                                    aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTCCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTccc                                <  2:221152/139‑1 (MQ=255)
                                                                                                            tCCAGCGCCGCGCCAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc                        <  1:1302094/139‑1 (MQ=255)
                                                                                                                 cgccgcgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtctc                   >  2:989383/1‑139 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      >  1:59707/1‑83 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                      <  2:59707/83‑1 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTccc                                >  1:457721/1‑121 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTccc                                <  2:457721/121‑1 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc                        <  2:1061422/129‑1 (MQ=255)
                                                                                                                      cgccAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc                        >  1:1061422/1‑129 (MQ=255)
                                                                                                                         cAAGACCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGcaaca           <  2:189834/139‑1 (MQ=255)
                                                                                                                             aCCAATGATGCC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGc       <  1:1464847/139‑1 (MQ=255)
                                                                                                                                  tgatgcC‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttctt  <  2:935968/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑C‑CCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTT  >  NC_000913/1759541‑1759808

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