Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_hsa_tpiA 257,908 Δ776 bp [crl] [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_hsa_tpiA 257908–258675 258683 9–776 66 [0] [0] 65 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 2579070 (0.000)58 (1.140) 46/262 0.2 100% intergenic (+8/‑769) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913_3_hsa_tpiA 258684 = 0 (0.000)pseudogene (9/331 nt) crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

ACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_hsa_tpiA/257759‑258055
                                                                                                                                                    |                                                                                                                                                    
aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACAGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGACCGTATATTTGTGAAGGTA                                                                                                                                                      >  1:189932‑M1/1‑149 (MQ=255)
         ttAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                                                                                                             >  1:379051‑M1/1‑140 (MQ=255)
         ttAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                                                                                                             >  2:589276‑M1/1‑140 (MQ=255)
             aaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagatca                                                                                                                                         >  1:463371‑M1/1‑136 (MQ=255)
                   gCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAAAAACAGTCAATCGTGAAGGTAagtgcaaagataatcgatt                                                                                                                                   <  1:478885‑M1/149‑20 (MQ=255)
                     tAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                                                                                                                 >  1:923032‑M1/1‑128 (MQ=255)
                         tgttgcCGATTAGGTAAAACAGTTGCATAACAACAGGAGAGAGTAATGACGTTACCGAGTGCACATCCGAAGAGTAGATTTATCAAATTATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                                                                             <  2:189932‑M1/147‑26 (MQ=255)
                          tttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                                                                            <  2:826363‑M1/149‑27 (MQ=255)
                                      gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggct                                                                                                                <  1:851863‑M1/149‑39 (MQ=255)
                                       tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                                                                               >  1:279681‑M1/1‑110 (MQ=255)
                                       tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgaaaagataatcgattctttttcgattgtctggctg                                                                                                               >  2:719217‑M1/1‑110 (MQ=255)
                                          aaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcacagataatcgattctttttcgattgtctggctgtat                                                                                                            >  2:277189‑M1/1‑107 (MQ=255)
                                                 gCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                                                                     >  1:579960‑M1/1‑100 (MQ=255)
                                                    tCAAAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                  <  2:134188‑M1/149‑53 (MQ=255)
                                                    cgctaacaggggATAGGAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                  <  2:681092‑M1/145‑53 (MQ=255)
                                                         acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGTAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcgttttcgattgtctggctgtatgcgtcaacgtgaatc                                                                                             <  2:611899‑M1/149‑58 (MQ=255)
                                                         acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaacaaaaccgattccttttcgattgtctggctgtatgcgtcaacgtgaaac                                                                                             <  2:831189‑M1/149‑58 (MQ=255)
                                                               gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcac                                                                                       <  2:29915‑M1/149‑64 (MQ=255)
                                                                 tAGCAATGACGGTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccg                                                                                     <  1:231863‑M1/149‑66 (MQ=255)
                                                                  aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgg                                                                                    <  1:589276‑M1/149‑67 (MQ=255)
                                                                     aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaag                                                                                 >  1:553578‑M1/1‑80 (MQ=255)
                                                                      aTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagt                                                                                <  2:579960‑M1/149‑71 (MQ=255)
                                                                                 gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattct                                                                     >  2:374781‑M1/1‑68 (MQ=255)
                                                                                 gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagata                                                                                                                                          <  2:919663‑M1/80‑13 (MQ=255)
                                                                                 gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgaaaagata                                                                                                                                          >  1:919663‑M1/1‑68 (MQ=255)
                                                                                         acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaga                                                                                                                                            <  1:777733‑M1/70‑11 (MQ=255)
                                                                                         acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaga                                                                                                                                            >  2:777733‑M1/1‑60 (MQ=255)
                                                                                           ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                  >  2:324953‑M1/1‑58 (MQ=255)
                                                                                           ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                  <  1:324953‑M1/94‑37 (MQ=255)
                                                                                                  gCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggag                                                    <  1:374781‑M1/149‑99 (MQ=255)
                                                                                                   cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttttatttttttgttgtatgagtcgaagtgaaaccggcaccggaagttggtgaattctgaggctggtggatggagc                                                   >  2:54189‑M1/1‑50 (MQ=255)
                                                                                                      aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgg                                                >  1:30014‑M1/1‑47 (MQ=255)
                                                                                                        tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaa                                              <  2:572343‑M1/149‑105 (MQ=255)
                                                                                                         gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                            <  2:781714‑M1/118‑75 (MQ=255)
                                                                                                         gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                            >  1:781714‑M1/1‑44 (MQ=255)
                                                                                                          aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttttttttttcttgctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagc                                            >  2:228608‑M1/1‑43 (MQ=255)
                                                                                                            cAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtga                                                                          <  1:539704‑M1/117‑77 (MQ=255)
                                                                                                            cAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtga                                                                          >  2:539704‑M1/1‑41 (MQ=255)
                                                                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacgggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         >  1:399916‑M1/1‑40 (MQ=255)
                                                                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggctggtggatggagctggaagggga                                         >  2:228928‑M1/1‑40 (MQ=255)
                                                                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgagaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         >  1:911727‑M1/1‑40 (MQ=255)
                                                                                                             aaaaaaTCTACAGCATCAGGCCCGTATATGCGTGAAGGTAagttcaaatataatcgattctctttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         <  2:173519‑M1/149‑110 (MQ=255)
                                                                                                              aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcag                                        >  1:406479‑M1/1‑39 (MQ=255)
                                                                                                              aaaaaTTTACCGCACCTGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcag                                        <  2:62339‑M1/149‑111 (MQ=255)
                                                                                                                 aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         >  2:396706‑M1/1‑36 (MQ=255)
                                                                                                                 aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         <  1:396706‑M1/145‑110 (MQ=255)
                                                                                                                  aTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaat                                    >  1:925605‑M1/1‑35 (MQ=255)
                                                                                                                          cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                               >  2:240604‑M1/1‑27 (MQ=255)
                                                                                                                          cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                               <  1:240604‑M1/146‑120 (MQ=255)
                                                                                                                           aCTAGGCCCGTATATTCGTGAAGGTAagtgacaagaaaaccttttgtttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttta                           <  1:425828‑M1/149‑124 (MQ=255)
                                                                                                                                gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                            >  1:256819‑M1/1‑21 (MQ=255)
                                                                                                                                gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                            <  2:256819‑M1/95‑75 (MQ=255)
                                                                                                                                gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgactggctttatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac                      >  1:758073‑M1/1‑21 (MQ=255)
                                                                                                                                 cccGTATATCCGTGAAGGTAagtggaaagatactcgattttttttagattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattccggggctggtggatggagctggaagcgaaggaatcccgttttacctaca                     <  2:132703‑M1/149‑130 (MQ=255)
                                                                                                                                 cccGTATAGTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca                     <  2:379051‑M1/149‑130 (MQ=255)
                                                                                                                                  ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacag                    <  1:277189‑M1/149‑131 (MQ=255)
                                                                                                                                           cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgaggaaggtggatggtgctggaagcgcaggaatcccgttttacctacagttaccagtt           <  1:672310‑M1/149‑140 (MQ=255)
                                                                                                                                            gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                               >  2:827533‑M1/1‑9 (MQ=255)
                                                                                                                                            gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgt                               <  1:827533‑M1/128‑120 (MQ=255)
                                                                                                                                                 gggaagtgcaaagataatcgattctttatcgagtgtctggctgtatgcgtgaacgtgaagccggcaccggaagtgcgtgaattctggggctgggggatggagcgggaagcgaaggaatcccgttttacctacagttaccagtttggtct     <  2:818994‑M1/146‑146 (MQ=255)
                                                                                                                                                 ggTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtct     >  2:678423‑M1/1‑4 (MQ=255)
                                                                                                                                                  gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggcctg    >  2:892255‑M1/1‑3 (MQ=255)
                                                                                                                                                   tAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggt                                                            <  2:484301‑M1/92‑91 (MQ=255)
                                                                                                                                                   tAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggt                                                            >  1:484301‑M1/1‑2 (MQ=255)
                                                                                                                                                    aagtgcaaagataatctattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtt  <  2:463371‑M1/149‑149 (MQ=255)
                                                                                                                                                    aagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggtggttttggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtt  <  1:228608‑M1/149‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_hsa_tpiA/257759‑258055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: