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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908 258683 776 60 [0] [0] 61 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579070 (0.000)57 (0.710) 42/280 1.5 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 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

CTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTC  >  NC_000913/258537‑258832
                                                                                                                                                   |                                                                                                                                                    
gcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                                                                                                                                     <  1:1063999‑M2/2‑1 (MQ=255)
gcttgtgttaaaaatttgctaaattttgccaagttggtaaaacagttgcattacaacaggggatagcaatgacgttaccgagtggacactcgaagtgaagattaaaaaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                                                                                                                                     <  2:1385617‑M2/2‑1 (MQ=255)
  ttgtgctaaaaattttctaaattttgtcaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                                                                                                   <  2:1255554‑M2/4‑1 (MQ=255)
        taaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                                                                                                             <  2:884462‑M2/10‑1 (MQ=255)
         gaaaatgtgctagattgtgccaatttggtaaaacagttgcatcacaacaggggatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGAt                                                                                                                                            <  2:267531‑M2/11‑1 (MQ=255)
            aatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                                                                                         <  1:798489‑M2/14‑1 (MQ=255)
            aagttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                                                                                         <  2:316552‑M2/14‑1 (MQ=255)
             atttgccaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgctaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                                                                                                                        <  1:1154820‑M2/15‑1 (MQ=255)
                 gctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAtt                                                                                                                                    <  2:858231‑M2/19‑1 (MQ=255)
                  ctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                                                                                                                                   >  2:928362‑M2/130‑149 (MQ=255)
                   taaattttgccaatttggtaaaacagttgcatcacaacaggagatagcactgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGCTAATCGAATCt                                                                                                                                  >  2:349418‑M2/129‑149 (MQ=255)
                   taaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCt                                                                                                                                  >  1:925838‑M2/129‑149 (MQ=255)
                   taaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATAgattat                                                                                                                                  >  2:1390321‑M2/129‑147 (MQ=255)
                     aattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATCCttt                                                                                                                                >  2:607161‑M2/127‑149 (MQ=255)
                       ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgccgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              <  1:1224337‑M2/25‑1 (MQ=255)
                       ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              >  2:393143‑M2/125‑149 (MQ=255)
                       ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              >  1:109666‑M2/125‑149 (MQ=255)
                       tttcgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              <  1:519505‑M2/25‑1 (MQ=255)
                          tgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGa                                                                                                                           >  2:658696‑M2/122‑149 (MQ=255)
                               atttggtaaaacagttgcatcacaacaggagatagcaattacgttaccgagtggacgcccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                      >  2:301900‑M2/117‑149 (MQ=255)
                               atttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaagtaAATGCACAGATAATCGATTCTTTTTCTCTTgtt                                                                                                                      >  1:1357223‑M2/117‑148 (MQ=255)
                                          cagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaggagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc                                                                                                           <  1:1309897‑M2/44‑1 (MQ=255)
                                            ggtgcatcacaacaagagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGt                                                                                                         <  1:1055889‑M2/46‑1 (MQ=255)
                                             ttgcatcacaacaggagatagcaatgacgttaccgggtggacacccgcagaggagagtgatcaaaaaatttacggcactaggcccgtatattcgtgacggtaAGTGCAAAGATCATCGATTCTTTTTCGATTGTCTGGCTGTAGGCGTc                                                                                                        >  2:6852‑M2/103‑149 (MQ=255)
                                                catcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                                                                                                     <  1:393959‑M2/50‑1 (MQ=255)
                                                 atcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactgggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTGTTTTTCGATTGTCTGGCTGTATGCGTCAACg                                                                                                    >  1:192610‑M2/99‑149 (MQ=255)
                                                 atcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg                                                                                                    <  1:674040‑M2/51‑1 (MQ=255)
                                                  tcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                                                   >  1:571550‑M2/98‑149 (MQ=255)
                                                     caacaggagatagcaatgtcgttaccgagtggacacccgacgagcagatataccaaaaaatttacggcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTCTTTCGATTGTCTGGCGGTATGCGTCAACGTGaa                                                                                                >  1:1023770‑M2/95‑149 (MQ=255)
                                                                tgacacgacgtgaccgagtggacaccctaagagcagattgatcaaaaaatttaccgcactgggcccgtatattcgtgaaggtaAGTGCAAAGATAAGCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCgg                                                                                     <  2:493658‑M2/66‑1 (MQ=255)
                                                                   aatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGGTTGTCTGGCTGTATGCCGCAACGGGAAACGGGCACCGGAAg                                                                                  >  2:1274244‑M2/81‑149 (MQ=255)
                                                                    atgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGt                                                                                 <  1:658696‑M2/70‑1 (MQ=255)
                                                                       acgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCg                                                                              <  1:928362‑M2/73‑1 (MQ=255)
                                                                         gttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgccctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGCAACCGGCACCGGAAGTGCGTg                                                                            >  2:530445‑M2/75‑149 (MQ=255)
                                                                          ttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGa                                                                           <  2:1249706‑M2/76‑1 (MQ=255)
                                                                                 gtggacacccgaagagcagattgatcaaacaatttaccgcactaggcccgtatatacgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTgg                                                                    <  2:815164‑M2/83‑1 (MQ=255)
                                                                                 gtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTgg                                                                    <  1:791424‑M2/83‑1 (MQ=255)
                                                                                 gtggacacccgaagagcagattgagcaaaaaatttaccgcactcggcccgtatattcgtgaaggtaAGTGCAAAGATTATCGATTCTTTTTCTATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAGTTCTgg                                                                    <  2:902888‑M2/83‑1 (MQ=255)
                                                                                       acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtgg                                                              <  1:1368663‑M2/89‑1 (MQ=255)
                                                                                       acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGGGAGttctgggggcgg                                                              >  2:30658‑M2/61‑145 (MQ=255)
                                                                                       acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatagtcgtgaaggtaAGTGCAACGATAATCGATTCTGTTTCGATTGTCTGGCTGTATGCGGCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtgg                                                              >  2:352560‑M2/61‑149 (MQ=255)
                                                                                         ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggtg                                                            >  2:800298‑M2/59‑149 (MQ=255)
                                                                                            aagagcagaatgaacaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGTAGTGCGTGAATTCTGGGGCTGTtggat                                                         <  2:1166750‑M2/94‑1 (MQ=255)
                                                                                                gcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAg                                                     >  1:918276‑M2/52‑149 (MQ=255)
                                                                                                 cagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGc                                                    >  1:484045‑M2/51‑149 (MQ=255)
                                                                                                 cagattgatcaaaaaatttaacgaactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATGGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGCAGTGCGTGAATTCTGGGGCCGGGGGATGGAGc                                                    >  2:328882‑M2/51‑149 (MQ=255)
                                                                                                   gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCAACACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTg                                                  >  2:639021‑M2/49‑149 (MQ=255)
                                                                                                           aaaaaatttaccgcactaggccggtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTGCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCa                                          >  1:247226‑M2/41‑149 (MQ=255)
                                                                                                           aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTGAACGTGAAACCGGCACCGGAAGTGGGTGAATTCTGGGGCTGGTGGATGGAGCAGGAagggcc                                          >  1:249504‑M2/41‑148 (MQ=255)
                                                                                                                atttacctcactcggcccgtatattcgtgaaggtaCGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAt                                     <  2:1023770‑M2/114‑1 (MQ=255)
                                                                                                                atttaccccaccaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAt                                     <  1:1276722‑M2/114‑1 (MQ=255)
                                                                                                                          ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTGCTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                           <  1:357623‑M2/124‑1 (MQ=255)
                                                                                                                               cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTGCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACa                      >  1:680778‑M2/21‑149 (MQ=255)
                                                                                                                                ccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCtacgg                     >  2:586926‑M2/20‑147 (MQ=255)
                                                                                                                                    atatccgggaaggtaAGTGCAAAGATATTCGATTCTTTTTCGAGTGTCGGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc                 <  1:62688‑M2/134‑1 (MQ=255)
                                                                                                                                          gtgaatgtaAGTGCAAAGATAATCGATTCGTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt           <  2:210983‑M2/140‑1 (MQ=255)
                                                                                                                                          gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAGTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt           <  1:146811‑M2/140‑1 (MQ=255)
                                                                                                                                             aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGt        <  2:571550‑M2/143‑1 (MQ=255)
                                                                                                                                             aaggtaAGTGCAAAGATAATCGATTCTTTTTCCATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGt        <  1:71923‑M2/143‑1 (MQ=255)
                                                                                                                                                gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTg     >  1:296850‑M2/4‑149 (MQ=255)
                                                                                                                                                   aGTGGAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTc  <  2:1226001/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTC  >  NC_000913/258537‑258832

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