Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,760,433 C→T 5.2% A31T (GCC→ACC)  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,760,4330CT 5.2% 260.0 / 2.3 77A31T (GCC→ACC) sufScysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCTTGCGCACGTTCTCCATTTTCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGTTTAAGC  >  NC_000913/1760298‑1760577
                                                                                                                                       |                                                                                                                                                
cccTCTTGCGCACGTTCTCCATTTTCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCgag                                                                                                                                      <  5:309024/148‑1 (MQ=255)
    cTTGCGCACGTTCTCCATTTTCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGTGGCGCTGTCGAGAt                                                                                                                                    >  1:88790/1‑146 (MQ=255)
         gcACGTTCTCCATTTTCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAg                                                                                                                             >  7:55863/1‑148 (MQ=255)
                  ccATTTTCTCGGTCGCCTGGGCGATTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAAcc                                                                                                                    >  8:120044/1‑148 (MQ=255)
                       ttCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGttt                                                                                                                <  7:377353/147‑1 (MQ=255)
                        tctcGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGTGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTa                                                                                                               <  2:88790/147‑1 (MQ=255)
                                 ccTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTTCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGa                                                                                                     <  2:280400/148‑1 (MQ=255)
                                      gcgcTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTGCTGCGCACTGGCGGCGCTGTCGAGATCAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCa                                                                                                <  8:174516/148‑1 (MQ=255)
                                            aaGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGc                                                                                                                         >  1:301157/1‑117 (MQ=255)
                                            aaGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGc                                                                                                                         <  2:301157/117‑1 (MQ=255)
                                            aaGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGa                                                                                           >  3:350003/1‑147 (MQ=255)
                                                   gAATACCACGATGTACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCaaa                                                                                                                      >  5:316879/1‑113 (MQ=255)
                                                   gAATACCACGATGTACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCaaa                                                                                                                      <  6:316879/113‑1 (MQ=255)
                                                   gAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATACGCCAGCGGCAACCCGTTTACCTCACGCGCAAGCACCGGAAAGTCGGc                                                                                   >  4:19151/1‑148 (MQ=255)
                                                   gAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGttt                                                                                                                >  7:261225/1‑119 (MQ=255)
                                                   gAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGttt                                                                                                                <  8:261225/119‑1 (MQ=255)
                                                   gAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGc                                                                                   >  4:14091/1‑148 (MQ=255)
                                                    aaTACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCaaa                                                                                                                      <  8:305627/112‑1 (MQ=255)
                                                    aaTACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCaaa                                                                                                                      >  7:305627/1‑112 (MQ=255)
                                                    aaTACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGcc                                                                                  >  7:251856/1‑148 (MQ=255)
                                                       aCCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTc                                                                                      <  5:13046/141‑1 (MQ=255)
                                                       aCCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTc                                                                                      >  6:13046/1‑141 (MQ=255)
                                                       aCCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCg                                                                                >  6:182789/1‑147 (MQ=255)
                                                          aCGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGTGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACt                                                                            >  4:361072/1‑148 (MQ=255)
                                                          aCGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACt                                                                            >  4:101865/1‑148 (MQ=255)
                                                           cGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACtt                                                                           <  5:182789/148‑1 (MQ=255)
                                                            gATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCa                                                                                                <  3:245263/126‑1 (MQ=255)
                                                            gATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCa                                                                                                >  4:245263/1‑126 (MQ=255)
                                                                cACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCg                                                                      <  8:55863/148‑1 (MQ=255)
                                                                        cgTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGaaaa                                                              <  8:251856/148‑1 (MQ=255)
                                                                              cATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatcat                                                        >  2:234255/1‑148 (MQ=255)
                                                                              cATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatcat                                                        <  7:120044/148‑1 (MQ=255)
                                                                                 gACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTg                                                                         >  1:13781/1‑128 (MQ=255)
                                                                                 gACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTg                                                                         <  2:13781/128‑1 (MQ=255)
                                                                                 gACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCtt                                                     >  7:103421/1‑148 (MQ=255)
                                                                                  aCGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCtt                                                     <  1:234255/147‑1 (MQ=255)
                                                                                   cGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGc                                                   <  8:103421/148‑1 (MQ=255)
                                                                                     aTAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCa                                                  <  1:195908/147‑1 (MQ=255)
                                                                                     aTAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCa                                                  >  2:195908/1‑147 (MQ=255)
                                                                                        aaaCTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGTGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcctc                                              <  3:361072/148‑1 (MQ=255)
                                                                                          aCTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcctcct                                            >  1:145451/1‑148 (MQ=255)
                                                                                                cctcggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGc                                                                                   <  8:365379/103‑1 (MQ=255)
                                                                                                cctcggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGc                                                                                   >  7:365379/1‑103 (MQ=255)
                                                                                                  tcggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCg                                    >  1:103359/1‑148 (MQ=255)
                                                                                                   cggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatca                                                         <  5:265771/126‑1 (MQ=255)
                                                                                                   cggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatca                                                         >  6:265771/1‑126 (MQ=255)
                                                                                                   cggcGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGt                                   >  6:83696/1‑148 (MQ=255)
                                                                                                         cAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCa                                                  >  5:368067/1‑127 (MQ=255)
                                                                                                         cAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCa                                                  <  6:368067/127‑1 (MQ=255)
                                                                                                            tCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcctcct                                            <  7:57859/130‑1 (MQ=255)
                                                                                                            tCACCTGGCTAGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcctcct                                            >  8:57859/1‑130 (MQ=255)
                                                                                                             cACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGAt                           >  8:181779/1‑146 (MQ=255)
                                                                                                             cACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGAt                           <  7:181779/146‑1 (MQ=255)
                                                                                                               ccTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCgg                       <  6:173865/148‑1 (MQ=255)
                                                                                                                cTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcc                                                >  1:11619/1‑122 (MQ=255)
                                                                                                                cTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCAcc                                                <  2:11619/122‑1 (MQ=255)
                                                                                                                cTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCggg                      >  6:141667/1‑148 (MQ=255)
                                                                                                                 tGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGc                     >  6:240021/1‑148 (MQ=255)
                                                                                                                 tGGCTCGGTTTCTGCGCACTGGCGGCGCAGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACACATTCTGGc                     >  5:19206/1‑148 (MQ=255)
                                                                                                                   gCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCa                   <  5:83696/148‑1 (MQ=255)
                                                                                                                    cTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCtt                                                     >  8:333376/1‑113 (MQ=255)
                                                                                                                    cTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCtt                                                     <  7:333376/113‑1 (MQ=255)
                                                                                                                    cTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAg                  >  5:72908/1‑148 (MQ=255)
                                                                                                                    cTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGGCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACTGATTCGGGc                     >  7:303617/1‑145 (MQ=255)
                                                                                                                     tCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGcc                    <  8:303617/145‑1 (MQ=255)
                                                                                                                       ggTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCa                <  4:350003/147‑1 (MQ=255)
                                                                                                                       ggTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCAc               <  3:14091/148‑1 (MQ=255)
                                                                                                                       ggTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCAc               <  2:145451/148‑1 (MQ=255)
                                                                                                                       ggTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCAc               <  3:19151/148‑1 (MQ=255)
                                                                                                                        gTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTg                                 <  3:103391/129‑1 (MQ=255)
                                                                                                                        gTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTg                                 >  4:103391/1‑129 (MQ=255)
                                                                                                                               cgcACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACctgctg         <  2:103359/146‑1 (MQ=255)
                                                                                                                               cgcACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGtt       <  6:72908/148‑1 (MQ=255)
                                                                                                                                gcACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGttt      <  5:141667/148‑1 (MQ=255)
                                                                                                                                gcACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGttt      <  5:240021/148‑1 (MQ=255)
                                                                                                                                    tGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGTTTAAGc  >  4:84172/1‑148 (MQ=255)
                                                                                                                                    tGGCGGCGCTGGCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGAt                           <  3:113184/123‑1 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCg                                                                      <  8:129371/79‑1 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCg                                                                      >  7:129371/1‑79 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatc                                                          >  5:190808/1‑91 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAatc                                                          <  6:190808/91‑1 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTc                         >  6:1876/1‑124 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTc                         <  5:1876/124‑1 (MQ=255)
                                                                                                                                     ggcggcGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGTTTaa    >  4:168251/1‑145 (MQ=255)
                                                                                                                                      gcggcgCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGTTTAAg   <  3:168251/145‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CCCGCTTGCGCACGTTCTCCATTTTCTCGGTCGCCTGGGCGCTTAAGGTATGAATACCACGATGCACCGCCGCGTAGCCATGACGATAAAACTCGGCCTCGGCGTCAATCACCTGGCTCGGTTTCTGCGCACTGGCGGCGCTGTCGAGATAAGCCAGCGGCAAACCGTTTACCTCACGCGAAAGCACCGGAAAGTCGGCCCGCACTTTGTCGACGGAAAAAATCATCTTGCACCTCCTGGCAGCCGTTGACCGATTCGGGCCAGCACCTGCTGTTTAAGC  >  NC_000913/1760298‑1760577

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