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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_pae_pgi 257908 258683 776 67 [0] [0] 66 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 2579070 (0.000)57 (0.890) 48/262 0.4 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_pae_pgi 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_pae_pgi/257759‑258055
                                                                                                                                                    |                                                                                                                                                    
aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA                                                                                                                                                      >  1:1046223‑M1/1‑149 (MQ=255)
 cGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAa                                                                                                                                                     <  2:105505‑M1/149‑2 (MQ=255)
     tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTAACGCACTAGGCCCGTATATTCGTGAAGGTAagtgc                                                                                                                                                 >  2:631152‑M1/1‑144 (MQ=255)
         ttAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                                                                                                                             >  2:709678‑M1/1‑140 (MQ=255)
                 ttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                                                                                                                     <  1:185291‑M1/149‑18 (MQ=255)
                      aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                                <  2:1046223‑M1/149‑23 (MQ=255)
                       aaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttt                                                                                                                               >  2:611146‑M1/1‑126 (MQ=255)
                         tcttgcCAATTTGGTAAAACAGTTGCAGCACAACAGGAGATAGCAATGACGTTACCGAGGGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                                                                             <  1:368025‑M1/147‑26 (MQ=255)
                               cAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGTGAagtgcacagataatcgattctttttcgattg                                                                                                                       >  2:766115‑M1/1‑118 (MQ=255)
                                aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                                                                                                      >  2:587037‑M1/1‑117 (MQ=255)
                                aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                                                                                                      >  2:204882‑M1/1‑117 (MQ=255)
                                      gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggct                                                                                                                >  2:396096‑M1/1‑111 (MQ=255)
                                        aaaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGCTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattctttttcgattgtctggctgt                                                                                                              >  1:366171‑M1/1‑109 (MQ=255)
                                          aaCAGTTGCAACACAACAAGATATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTAACGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                            <  2:958626‑M1/149‑43 (MQ=255)
                                             aGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTAATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattggctggctgtatgcg                                                                                                         >  2:277708‑M1/1‑104 (MQ=255)
                                                 gCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                                                                     <  2:710790‑M1/149‑50 (MQ=255)
                                                      acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtga                                                                                                >  2:15180‑M1/1‑95 (MQ=255)
                                                        aacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCCCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgactctttttcgattgtctggctgtatgcgtcaacgtgaaa                                                                                              <  1:129645‑M1/149‑57 (MQ=255)
                                                             gagaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                                                                         >  2:167946‑M1/1‑88 (MQ=255)
                                                               gaTAGCAATGACGTTACCGAGTGGACACGCGAAGAGCAGATTGATCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcac                                                                                       <  2:220139‑M1/149‑64 (MQ=255)
                                                                  aGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgg                                                                                    >  2:711180‑M1/1‑83 (MQ=255)
                                                                    cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaa                                                                                  <  1:204882‑M1/149‑69 (MQ=255)
                                                                     aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATGTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccagaag                                                                                 >  1:604672‑M1/1‑80 (MQ=255)
                                                                       tGACGTAACCGACTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtg                                                                               >  2:146370‑M1/1‑78 (MQ=255)
                                                                                 gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattct                                                                     >  1:448738‑M1/1‑68 (MQ=255)
                                                                                   gTGGACACCCGAAGAGAAAAATAAAAAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgggaactctgg                                                                   <  2:1053018‑M1/149‑84 (MQ=255)
                                                                                         acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgg                                                             <  1:711180‑M1/149‑90 (MQ=255)
                                                                                         acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgg                                                             >  2:640071‑M1/1‑60 (MQ=255)
                                                                                         acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatgctttttcgattgtatcgctgtatgtgtcaacgtgaaaccggcaccggaagttcgtgaattctggggctgg                                                             >  2:1022940‑M1/1‑60 (MQ=255)
                                                                                          cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggt                                                            <  1:685412‑M1/149‑91 (MQ=255)
                                                                                           ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                                                                          <  1:771146‑M1/102‑45 (MQ=255)
                                                                                           ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                                                                          >  2:771146‑M1/1‑58 (MQ=255)
                                                                                           ccGAAGAGCAGATTGATCAAAAAATTTACCGCACGAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtgaacgtgaaaccggcacgggacgtgcgtgaattctggggctgttg                                                           >  2:959943‑M1/1‑58 (MQ=255)
                                                                                            cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgg                                                          >  1:259445‑M1/1‑57 (MQ=255)
                                                                                            cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg                                                                                                                   <  2:809596‑M1/92‑36 (MQ=255)
                                                                                            cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg                                                                                                                   >  1:809596‑M1/1‑57 (MQ=255)
                                                                                              aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggg                                                                  <  2:501334‑M1/139‑85 (MQ=255)
                                                                                              aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggg                                                                  >  1:501334‑M1/1‑55 (MQ=255)
                                                                                                     agttaatccaaAAATTTACCGCACCAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctg                                                 <  1:254157‑M1/144‑102 (MQ=255)
                                                                                                      aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattatttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgagttctggggctggtggatggagctgg                                                >  2:739040‑M1/1‑47 (MQ=255)
                                                                                                       ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgga                                               >  1:761693‑M1/1‑46 (MQ=255)
                                                                                                          aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagc                                            >  2:545992‑M1/1‑43 (MQ=255)
                                                                                                          aTCAAAAAATTTACCGCACTAGGCCCGAATCCTCGTGAAGGTAagtgcaaagataatcgaatctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggtgctgctggatggagccagaagc                                            >  2:500126‑M1/1‑43 (MQ=255)
                                                                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         >  2:1153575‑M1/1‑40 (MQ=255)
                                                                                                             aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatggtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                         >  1:921823‑M1/1‑40 (MQ=255)
                                                                                                              aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcag                                        <  1:709678‑M1/149‑111 (MQ=255)
                                                                                                               aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcagg                                       <  1:184585‑M1/149‑112 (MQ=255)
                                                                                                               aaaaTTTACCGCACTAGGCCCGTATAGTCGTGAAGGTAagtgcaaagataatcgattcttttgcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcagg                                       <  2:1174301‑M1/149‑112 (MQ=255)
                                                                                                                          cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggcagcgcaggaatcccgtttt                            >  2:967091‑M1/1‑27 (MQ=255)
                                                                                                                            cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttac                          <  1:396096‑M1/149‑125 (MQ=255)
                                                                                                                            cTAGGCCCGTATATTCGTGAAGGTAagggcaaagataatcgattctttttcgattgtccggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgcgaattctggggctggtggatggagctggaagcgcaggaatcccgttttac                          <  1:936925‑M1/149‑125 (MQ=255)
                                                                                                                                 cccGTATATTCGTGAAGGTAagtgcaaagataaacgaatcttaatagagtttttggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca                     <  1:277708‑M1/149‑130 (MQ=255)
                                                                                                                                  ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacag                    >  1:1168453‑M1/1‑19 (MQ=255)
                                                                                                                                    gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagtt                  <  1:745169‑M1/149‑133 (MQ=255)
                                                                                                                                     tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagtta                 <  1:640071‑M1/149‑134 (MQ=255)
                                                                                                                                        atTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctccagttacca              >  2:458797‑M1/1‑13 (MQ=255)
                                                                                                                                        atTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacca              >  1:4327‑M1/1‑13 (MQ=255)
                                                                                                                                           cGTGAAGGTAagtgcaaagataatcgatgctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtaaattctggggctggtggatggagctgtaagcgcaggaatcccgttttacctacagttaccagtt           <  2:514398‑M1/149‑140 (MQ=255)
                                                                                                                                             ttaaggTAatggcaaagataatcgattctttttcgaatgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttg         <  1:766115‑M1/147‑142 (MQ=255)
                                                                                                                                             tGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttg         <  1:611146‑M1/149‑142 (MQ=255)
                                                                                                                                                aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc      <  1:1153575‑M1/149‑145 (MQ=255)
                                                                                                                                                 ggTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtct     <  1:808491‑M1/149‑146 (MQ=255)
                                                                                                                                                  gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggcagtgcgtgaattctggggctggtcgaagtcatgttaatcgtaggactcccgttttacctacagtgaccagtttggtctg    >  1:277603‑M1/1‑3 (MQ=255)
                                                                                                                                                   tAagtgcaaagataatcgattcttttttttttttttttctgtatgcgtaaacgttaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttctcagttttgtcttt   >  2:911662‑M1/1‑2 (MQ=255)
                                                                                                                                                    aagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggt       <  2:814870‑M1/144‑144 (MQ=255)
                                                                                                                                                    aagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggt       >  1:814870‑M1/1‑1 (MQ=255)
                                                                                                                                                    aagtgcaaagataatcgattctttttcgagtgtctggctggatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtt  <  2:838504‑M1/149‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_pae_pgi/257759‑258055

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: