Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,826 +A intergenic (+1/+44) yjgL → / ← argI hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,8231.A100.0% 485.3 / NA 133*617K (TAA→AAA) yjgLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (94/39);  total (94/39)

Only 100 of 141 total aligned reads displayed.

AATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474733‑4474910
                                                                                           |                                                                                       
aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataa                                                                                       >  1:4202174/1‑93 (MQ=255)
aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataa                                                                                       >  2:4422084/1‑93 (MQ=255)
  tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaa                                                                                     >  2:1940986/1‑92 (MQ=255)
  tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaa                                                                                     <  1:1883875/93‑2 (MQ=255)
   gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaat                                                                                    >  1:3845969/1‑91 (MQ=255)
   gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaat                                                                                    >  2:1438667/1‑91 (MQ=255)
     aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaataa                                                                                  >  1:4660487/1‑89 (MQ=255)
        gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:3195250/1‑93 (MQ=255)
        gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:911726/1‑93 (MQ=255)
        gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  2:4408131/1‑93 (MQ=255)
        gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  2:3659226/1‑93 (MQ=255)
          gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:2077610/1‑93 (MQ=255)
          gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:3701567/1‑93 (MQ=255)
          gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:1239359/1‑93 (MQ=255)
          gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              <  1:19385/93‑1 (MQ=255)
          gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              <  1:1170769/93‑1 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGGCCTATCACTGCATAAAATAATGTGAt                                                                             >  2:3794904/1‑93 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:4235568/1‑93 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             <  2:3849345/93‑1 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:4516209/1‑93 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:1121119/1‑93 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             <  2:3116556/93‑1 (MQ=255)
           tttAAATAAACTAAATCCACCGGAAAATTTACGGATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:1070353/1‑93 (MQ=255)
                 taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCggg                                                                       >  1:2110197/1‑93 (MQ=255)
                 taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCggg                                                                       >  2:1771206/1‑93 (MQ=255)
                 taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCggg                                                                       <  1:2146860/93‑1 (MQ=255)
                  aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      >  1:2933855/1‑93 (MQ=255)
                  aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      <  2:3932604/93‑1 (MQ=255)
                   aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  2:4332184/1‑93 (MQ=255)
                           ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             <  2:2986233‑M1/93‑4 (MQ=255)
                           ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             <  2:3185685‑M1/93‑4 (MQ=255)
                            cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            >  1:4474425‑M1/1‑89 (MQ=255)
                            cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            >  1:894148‑M1/1‑89 (MQ=255)
                             aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                           >  2:2238406‑M1/1‑88 (MQ=255)
                              ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                          >  2:4482590‑M1/1‑87 (MQ=255)
                                ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:1425747‑M1/1‑85 (MQ=255)
                                ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  2:1323879‑M1/1‑85 (MQ=255)
                                 gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  2:4243506‑M1/1‑84 (MQ=255)
                                  aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  2:1540049‑M1/93‑11 (MQ=255)
                                     aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                   >  2:2499401‑M1/1‑80 (MQ=255)
                                     aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                   <  2:3982757‑M1/93‑14 (MQ=255)
                                      tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  2:3707926‑M1/1‑79 (MQ=255)
                                         aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  1:2219465‑M1/93‑18 (MQ=255)
                                         aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  1:2774652‑M1/93‑18 (MQ=255)
                                          cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  2:4185529‑M1/1‑75 (MQ=255)
                                          cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              <  2:153976‑M1/93‑19 (MQ=255)
                                            tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            >  1:2270879‑M1/1‑73 (MQ=255)
                                            tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            >  1:3182610‑M1/1‑73 (MQ=255)
                                            tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            >  2:1837729‑M1/1‑73 (MQ=255)
                                            tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            <  2:8879‑M1/93‑21 (MQ=255)
                                             ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  2:2991770‑M1/1‑72 (MQ=255)
                                             ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  2:1965393‑M1/1‑72 (MQ=255)
                                             ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  2:2066895‑M1/93‑22 (MQ=255)
                                              taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          >  1:2737823‑M1/1‑71 (MQ=255)
                                               aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                         >  1:1676271‑M1/1‑70 (MQ=255)
                                                gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  1:1134889‑M1/1‑69 (MQ=255)
                                                gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACCGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  2:1580292‑M1/1‑69 (MQ=255)
                                                 gaatAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  1:2254500‑M1/92‑26 (MQ=255)
                                                 cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       >  2:3658637‑M1/1‑68 (MQ=255)
                                                  aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      <  2:972726‑M1/93‑27 (MQ=255)
                                                    tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    <  2:826165‑M1/93‑29 (MQ=255)
                                                      gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  2:2483862‑M1/1‑63 (MQ=255)
                                                      gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  1:644844‑M1/1‑63 (MQ=255)
                                                      gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  1:180411‑M1/1‑63 (MQ=255)
                                                      gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  <  2:3267609‑M1/93‑31 (MQ=255)
                                                      gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGCGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  1:4331591‑M1/1‑63 (MQ=255)
                                                       aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 >  2:2780875‑M1/1‑62 (MQ=255)
                                                        acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  2:4057134‑M1/1‑61 (MQ=255)
                                                          acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              >  1:4228579‑M1/1‑59 (MQ=255)
                                                            tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTggatttcccggcatctttatagcgatagcaattatt                            >  1:1810042‑M1/1‑57 (MQ=255)
                                                              tGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                          <  1:2450993‑M1/93‑39 (MQ=255)
                                                              tGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGGTTgtatttcccggcatctttatagcgatagcaattattta                          >  2:2564435‑M1/1‑55 (MQ=255)
                                                                ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        >  2:3437744‑M1/1‑53 (MQ=255)
                                                                  tGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      <  2:235544‑M1/93‑43 (MQ=255)
                                                                   ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:4381195‑M1/1‑50 (MQ=255)
                                                                    gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    <  1:2220152‑M1/93‑45 (MQ=255)
                                                                    gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    >  1:109183‑M1/1‑49 (MQ=255)
                                                                      agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  2:3506969‑M1/1‑47 (MQ=255)
                                                                      agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  1:1735279‑M1/1‑47 (MQ=255)
                                                                      agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  <  2:3698694‑M1/93‑47 (MQ=255)
                                                                       gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  2:236800‑M1/93‑48 (MQ=255)
                                                                       gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  2:2510969‑M1/93‑48 (MQ=255)
                                                                       gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  2:2889571‑M1/1‑46 (MQ=255)
                                                                       gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  2:126941‑M1/1‑46 (MQ=255)
                                                                       gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  2:3291608‑M1/1‑46 (MQ=255)
                                                                         aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               >  1:2937785‑M1/1‑44 (MQ=255)
                                                                         aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               <  2:3470619‑M1/93‑50 (MQ=255)
                                                                         aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGGTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               >  2:2760704‑M1/1‑44 (MQ=255)
                                                                          tAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcg              >  2:1494802‑M1/1‑43 (MQ=255)
                                                                             aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacc           >  2:1878491‑M1/1‑40 (MQ=255)
                                                                               cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         >  1:4642013‑M1/1‑38 (MQ=255)
                                                                                tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  2:327715‑M1/1‑37 (MQ=255)
                                                                                tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        <  2:4607469‑M1/93‑57 (MQ=255)
                                                                                 aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatca       <  1:2453552‑M1/93‑58 (MQ=255)
                                                                                   cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     >  2:4380151‑M1/1‑34 (MQ=255)
                                                                                    aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg    >  1:682287‑M1/1‑33 (MQ=255)
                                                                                    aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg    <  2:2902582‑M1/93‑61 (MQ=255)
                                                                                     cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc   >  1:2602574‑M1/1‑32 (MQ=255)
                                                                                     cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc   <  1:1571763‑M1/93‑62 (MQ=255)
                                                                                      tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  1:3088476‑M1/1‑31 (MQ=255)
                                                                                           |                                                                                       
AATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474733‑4474910

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