Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,481 T→G 6.7% V58G (GTC→GGC)  ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,4810TG 6.7% 184.5 / 2.0 71V58G (GTC→GGC) ybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base G (2/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATG  >  NC_000913/774336‑774628
                                                                                                                                                 |                                                                                                                                                   
gCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGGTTTGCCGCGAAGACCAGTGCACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtg                                                                                                                                                  >  1:365309/1‑149 (MQ=255)
  ggTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgtg                                                                                                                                                >  8:363337/1‑149 (MQ=255)
     aaTGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTg                                                                                                                                             >  6:67132/1‑149 (MQ=255)
      aTGGACTAGCGCCCCCTGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGTAGACCTGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTg                                                                                                                                             <  8:153707/148‑1 (MQ=255)
      aTGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTg                                                                                                                                             <  1:347627/148‑1 (MQ=255)
            aaGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGGTTTGGGACCCAGCTCGTGTTGCCGCGAAGACGAGTGAACTGGAAATCGGGCAGGGTTTAGTGCGGAGGTGGGCCGGCGGTGCTGGTGTGAt                                                                                                                                      >  3:243769/1‑149 (MQ=255)
             cgcgcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCGTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAtt                                                                                                                                     <  6:81218/148‑1 (MQ=255)
                 cccccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACg                                                                                                                                 <  5:67132/149‑1 (MQ=255)
                     cTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGt                                                                                                                             <  5:253875/149‑1 (MQ=255)
                        cGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTggg                                                                                                                          <  1:35318/149‑1 (MQ=255)
                        cGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTg                                                                                                                            <  1:15134/147‑1 (MQ=255)
                        cGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTg                                                                                                                            <  5:39688/147‑1 (MQ=255)
                                 tCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTc                                                                                                                  <  3:327642/148‑1 (MQ=255)
                                   cTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCg                                                                                                                >  7:332968/1‑148 (MQ=255)
                                    ttCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGGTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCg                                                                                                                >  2:36330/1‑147 (MQ=255)
                                       gTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGGTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGaa                                                                                                           <  1:135668/149‑1 (MQ=255)
                                        tGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAg                                                                                                          >  8:368710/1‑149 (MQ=255)
                                           tGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGtt                                                                                                       <  7:363337/149‑1 (MQ=255)
                                            ggCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGGATGGTTTATTGCTGATGTGGGGCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCGGCAGAAGGGTc                                                                                                      >  3:308600/1‑149 (MQ=255)
                                              cgcTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGGTCtt                                                                                                    >  5:234144/1‑149 (MQ=255)
                                               gcTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCttt                                                                                                   <  5:406932/149‑1 (MQ=255)
                                               gcTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTATTGCCGCGAATACAAGAGAAATTTAAATCTGGCTTGGTTTATTGGTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTGCGCAGTAGGttggtg                                                                                                   >  5:375671/1‑145 (MQ=255)
                                                     gTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAg                                                                                             >  1:390617/1‑149 (MQ=255)
                                                      ttAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCaa                                                                                              <  5:72437/147‑1 (MQ=255)
                                                         gCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGAt                                                                                         >  8:111404/1‑149 (MQ=255)
                                                             gATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGTAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGAttttt                                                                                     <  7:111404/149‑1 (MQ=255)
                                                             gATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGAttttt                                                                                     <  6:234144/149‑1 (MQ=255)
                                                               tgtatgtTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgt                                                                                                                                           >  2:7531/1‑93 (MQ=255)
                                                               tgtatgtTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgt                                                                                                                                           <  1:7531/93‑1 (MQ=255)
                                                                  atgtTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTgg                                                                                                                           <  8:320099/106‑1 (MQ=255)
                                                                  atgtTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTgg                                                                                                                           >  7:320099/1‑106 (MQ=255)
                                                                     ttttGGGCCCCATCGCGTTTTGCCGCGACGACCAGTCCCCTGGAAATCTGGCATGGTTTACTGCTGACGGGGGCCGGCGGTGCTGGCGCGCATCACGGCGGGGGCTTTCGTCCGCAGAACGTTCTTTGGCAAGGGAGTTTTTGCCCAtt                                                                             <  4:243769/149‑1 (MQ=255)
                                                                        tGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGGCTGTGCTGGTGTGATTCACGGCGTGGGCtt                                                                                                                       >  2:55758/1‑104 (MQ=255)
                                                                        tGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGGCTGTGCTGGTGTGATTCACGGCGTGGGCtt                                                                                                                       <  1:55758/104‑1 (MQ=255)
                                                                           gACCCATCTCGGTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCAttgcttgc                                                                       <  2:337416/149‑1 (MQ=255)
                                                                              ccATCTCGTCTTGCCGCCACGACCAGTGAACTCGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGGGGGATTCACGGCGTGGGCTTTCGTCCGCCGAAGGGGCTTTGGCAAGGGATTTTTAGCCCATTGCTTGCCGa                                                                    <  2:365309/149‑1 (MQ=255)
                                                                                 gctcGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTGTATTGCTGAAGTGGGCCGTCTGGGCTGGTGTGATTCACTGCGTCGGCTGTCGTCCGCCGCACGTTCTTTGGCAAGGGAGTTTTTGGCCATTGCTTGCCGatta                                                                 >  2:208270/2‑147 (MQ=255)
                                                                                  ctcGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATAtt                                                                >  7:435038/1‑149 (MQ=255)
                                                                                   tcGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTg                                                               >  1:236669/1‑149 (MQ=255)
                                                                                     gTTTCGCCGCGAAGACCAGTGAACTGGAAATCTGGCAGGGTTTCTTGCTGATGGGGGCCGTCTGTGCTGGTGTGAGTCACGGCGGGGGCTGTCGTCGGCAGACGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGatat                                                                 <  4:328494/145‑1 (MQ=255)
                                                                                          gCCGGGTAGACCAGTGAACTGGAATTCTGGCTTGGTTTATTGCTGATGTGGGCGGTCTGTTCTGGTGTGATTCTCGGCGTGGGCTTTCGTCCGCAGATGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCa                                                         <  7:394835/148‑1 (MQ=255)
                                                                                            cgcgAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTg                                                      >  3:313282/1‑149 (MQ=255)
                                                                                             gcgAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCa                                                         <  8:360639/145‑1 (MQ=255)
                                                                                             gcgAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCa                                                         >  7:360639/1‑145 (MQ=255)
                                                                                              cgAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTc                                                    <  6:375671/149‑1 (MQ=255)
                                                                                                       gTGAACTGGAAATCTGGCATGCTTTATTGCTGATGTGGGCCGTCTGTGCTTGTGTGATTCACCGCGTGGGCTATCGTCCGCAGAAGGTTCTTTGGAAAGGGATTTTTTGCCCATTGCTTGCAGATAGTGTGCTCATTGTCGGGCTGAtt                                           >  4:109456/1‑149 (MQ=255)
                                                                                                        tGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGAttt                                          >  5:163776/1‑149 (MQ=255)
                                                                                                          aaCTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCGGATTttc                                        <  8:332968/149‑1 (MQ=255)
                                                                                                          aaCTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGGTGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTttc                                        <  3:109456/149‑1 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGtt                                                                                                       <  5:176210/83‑1 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGtt                                                                                                       >  6:176210/1‑83 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCa                                                         <  1:382327/129‑1 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCa                                                         >  2:382327/1‑129 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTttcttc                                     >  4:138763/1‑149 (MQ=255)
                                                                                                             tGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTttcttc                                     >  4:62101/1‑149 (MQ=255)
                                                                                                                         aTGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAAt                          >  1:221149/1‑148 (MQ=255)
                                                                                                                             tttATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTg                                                                                                  >  1:162438/1‑72 (MQ=255)
                                                                                                                             tttATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTg                                                                                                  <  2:162438/72‑1 (MQ=255)
                                                                                                                                 ttGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCa                                                                                                              <  7:223857/56‑1 (MQ=255)
                                                                                                                                 ttGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCa                                                                                                              >  8:223857/1‑56 (MQ=255)
                                                                                                                                 ttGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATtc                  >  5:151944/1‑148 (MQ=255)
                                                                                                                                   gCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATtctc                <  7:368710/148‑1 (MQ=255)
                                                                                                                                    cTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTtc                                                            <  3:90575/103‑1 (MQ=255)
                                                                                                                                    cTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTtc                                                            >  4:90575/1‑103 (MQ=255)
                                                                                                                                    cTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAAt                          >  3:281417/1‑137 (MQ=255)
                                                                                                                                    cTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAAt                          <  4:281417/137‑1 (MQ=255)
                                                                                                                                     tGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCtt              >  3:91006/1‑148 (MQ=255)
                                                                                                                                       aTGGGGGCCGGCTGTGCTGGGGGGCGTCACGGCGTGGGCGTTCGGCCGCAGAAGGGTCTGTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTaa           <  2:390617/149‑1 (MQ=255)
                                                                                                                                          tgGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTaaaa         <  1:36330/148‑1 (MQ=255)
                                                                                                                                             gCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAAtt     <  6:163776/149‑1 (MQ=255)
                                                                                                                                                gTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATg  <  2:236669/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATG  >  NC_000913/774336‑774628

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