Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,250,156 C→T 5.0% T201T (ACC→ACT malM → maltose regulon periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,250,1560CT 5.0% 257.1 / 2.3 79T201T (ACC→ACTmalMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base T (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTGAA  >  NC_000913/4250010‑4250302
                                                                                                                                                  |                                                                                                                                                  
gggCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGa                                                                                                                                                  <  3:168554/149‑1 (MQ=255)
 ggCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGAt                                                                                                                                                 <  1:139715/149‑1 (MQ=255)
 ggCAGCAAAAACTTTATGTACTCGTCATTACCAAGGAAAAAGATCTCCAGCACACGACCCAACTGCTCGACCCGCCTAAAGCCTATGACAAGGCCGTCGGTAACTCGCACCCGCAAATCCCCCATCCGTTTGCTCGTCATACAACCGAt                                                                                                                                                 >  2:166567/1‑149 (MQ=255)
   cagcaAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAACGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCTGTTGCTCGTCATACCACCGATgg                                                                                                                                               <  6:212212/149‑1 (MQ=255)
        aaaaCTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGCCCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTac                                                                                                                                          >  6:232029/1‑149 (MQ=255)
             tttATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaa                                                                                                                                     >  1:352728/1‑149 (MQ=255)
               tATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGGTATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaac                                                                                                                                    <  1:317449/148‑1 (MQ=255)
                  gTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTactgaaactga                                                                                                                                 >  6:206586/1‑148 (MQ=255)
                     cTGGTCTTTACCACGGAAAAAGATCTCCAGCAGCCGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGt                                                                                                                             >  1:80705/1‑149 (MQ=255)
                        ttCTTTACCACGGCAAACGATCTCCATCATAGGACCCCACCTCTCTACCCGTCTACATCCTAAGCCAACTGCGTCCCTCACTAGATCCCGGAAATCCCCGCTCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGa                                                                                                                           <  2:298397/147‑1 (MQ=255)
                         tCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGaaa                                                                                                                         <  1:238663/149‑1 (MQ=255)
                         tCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGaaa                                                                                                                         <  5:206586/149‑1 (MQ=255)
                               cctcGGAAAAAGATCTCGAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAAc                                                                                                                   <  1:166567/146‑1 (MQ=255)
                               ccACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAAc                                                                                                                   <  1:364265/149‑1 (MQ=255)
                                cACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGt                                                                                                                             <  6:189408/138‑1 (MQ=255)
                                cACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGt                                                                                                                             >  5:189408/1‑138 (MQ=255)
                                cACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGATAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCACGGATATCCCCGATACGGTTGATCGTCATACCACCGATGGCTAACTGAAACTGAAAGTGAAAACGAAct                                                                                                                  >  4:274201/1‑149 (MQ=255)
                                 aCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActc                                                                                                                 >  7:69828/1‑149 (MQ=255)
                                 aCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActc                                                                                                                 >  3:325310/1‑149 (MQ=255)
                                 aCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActc                                                                                                                 >  7:75133/1‑149 (MQ=255)
                                  cGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActc                                                                                                                 >  7:258376/1‑148 (MQ=255)
                                   ggAAAAGGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAACCTATGCCAAGGGCGGCCGTAACTCGCACCCGGGTATCCCCGCTCCGGTTGCTGGTCATTCCACCGCTGGCGTACGGAAACGGCAAGTGAAAACCAAcgccac                                                                                                              >  3:292501/1‑148 (MQ=255)
                                   ggAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActcca                                                                                                               <  7:370730/149‑1 (MQ=255)
                                         aGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctcca                                                                                                         <  4:325310/149‑1 (MQ=255)
                                                cagcagACGACGCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCgtgt                                                                                                   <  1:325277/148‑1 (MQ=255)
                                                cagcagACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGt                                                                                                                             >  2:177715/1‑122 (MQ=255)
                                                cagcagACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGt                                                                                                                             <  1:177715/122‑1 (MQ=255)
                                                  gcagACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACACCAGCTCCAGCGTGTTgg                                                                                                >  4:324351/1‑149 (MQ=255)
                                                      acgacCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAg                                                                                             >  3:290470/1‑148 (MQ=255)
                                                        gacCCAACTTCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGCAACTGAAAGTGAAATCGAACGCCAGCTCCCGCGTGTTGGGAGGAc                                                                                          >  2:216100/1‑149 (MQ=255)
                                                         acCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGAcc                                                                                         <  3:55135/149‑1 (MQ=255)
                                                         acCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGAcc                                                                                         >  7:199071/1‑149 (MQ=255)
                                                           ccAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCtt                                                                                       >  2:64524/1‑149 (MQ=255)
                                                           ccAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCtt                                                                                       <  4:333246/149‑1 (MQ=255)
                                                             aaCTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGa                                                                                                                           <  2:72038/111‑1 (MQ=255)
                                                             aaCTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGa                                                                                                                           >  1:72038/1‑111 (MQ=255)
                                                             aaCTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTa                                                                                     >  2:348479/1‑149 (MQ=255)
                                                               cTGCTCGACCCGGCTAAAGCCTATGCCATTGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTAtt                                                                                   <  2:308086/149‑1 (MQ=255)
                                                                 gcccGACCCGGCTAACCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTgg                                                                                <  4:292501/146‑1 (MQ=255)
                                                                 gCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTg                                                                                 >  3:317769/1‑149 (MQ=255)
                                                                 gCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTg                                                                                 >  7:99419/1‑149 (MQ=255)
                                                                   tCGACCCCGCTAAAGCCTATGCCAAGGGCGGCGGTTACTCGATCCCGGTAATCCCCGATCCGGGTGCTCGTCATACCACCCCTGGGTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCgt                                                                                                     <  1:269521/127‑1 (MQ=255)
                                                                     gACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACAGAAAGTGAAAACGCACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttc                                                                             >  6:202445/1‑149 (MQ=255)
                                                                      aCCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTGTATTTGGttct                                                                            >  1:235657/1‑149 (MQ=255)
                                                                       cccGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttctt                                                                           <  1:348479/149‑1 (MQ=255)
                                                                       cccGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttctt                                                                           <  8:258376/149‑1 (MQ=255)
                                                                        ccGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttcttc                                                                          <  7:341554/149‑1 (MQ=255)
                                                                        ccGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttcttc                                                                          >  8:341554/1‑149 (MQ=255)
                                                                           gCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGc                                                                       <  8:99419/149‑1 (MQ=255)
                                                                               aaGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGAGCCGGTTGATCGTCATACCCCCGATGGGTAAATGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACGTTTTTTTGGTTCTGCCGATCCa                                                                   >  4:283334/1‑149 (MQ=255)
                                                                                 gCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCCTACCACCGATGGCTTACTGAAACTGAAAGGGAAACCGAACTCCAGCTCCCGCGTGTTGGTAGGACCTGTAGTTGGTTCTTCCGCTCCAGc                                                                 >  4:223777/1‑149 (MQ=255)
                                                                                    tATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTcc                                                              >  1:31733/1‑149 (MQ=255)
                                                                                      tGCCAAGGGCGTCGGTACCGCGATCCCGGATATCCCCGATCCGGTTGCTCGGCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTcc                                                              <  3:223777/147‑1 (MQ=255)
                                                                                        ccAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGt                                                           >  7:253662/1‑148 (MQ=255)
                                                                                             ggCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGCTGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACggt                                                     <  1:216100/149‑1 (MQ=255)
                                                                                                gTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACggtagg                                                  <  4:275576/149‑1 (MQ=255)
                                                                                                gTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACggtagg                                                  >  3:275576/1‑149 (MQ=255)
                                                                                                     tAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGATAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAaca                                             <  7:152775/149‑1 (MQ=255)
                                                                                                             tCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCa                                     >  3:250423/1‑149 (MQ=255)
                                                                                                              cccGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCAc                                    <  2:80705/149‑1 (MQ=255)
                                                                                                               ccGGATATCCCCGATCCGGTTGCTCGTCATACCACCGAGGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTGGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCAc                                    <  5:232029/148‑1 (MQ=255)
                                                                                                                  gATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGc                                <  8:75133/149‑1 (MQ=255)
                                                                                                                  gATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGc                                <  8:69828/149‑1 (MQ=255)
                                                                                                                      tCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCtgtg                            <  4:250423/149‑1 (MQ=255)
                                                                                                                      tCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCtgt                             <  3:324351/148‑1 (MQ=255)
                                                                                                                        cccGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGTTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGc                          >  4:313100/1‑149 (MQ=255)
                                                                                                                        cccGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTgg                           <  8:199071/148‑1 (MQ=255)
                                                                                                                        cccGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGc                          >  4:343030/1‑149 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACTGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGt                                                           >  1:251994/1‑115 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACTGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGt                                                           >  1:246814/1‑115 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACTGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGt                                                           <  2:246814/115‑1 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACTGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGt                                                           <  2:251994/115‑1 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttctt                                                                           >  4:65139/1‑99 (MQ=255)
                                                                                                                         ccGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGttctt                                                                           <  3:65139/99‑1 (MQ=255)
                                                                                                                          cGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTg                        >  6:130149/1‑149 (MQ=255)
                                                                                                                             tCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCAc                     <  4:290470/149‑1 (MQ=255)
                                                                                                                             tCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAAATGTAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGATCCTGTTACGGTAGGTAACACGGCGGCACCAGCTGTTGCTGCAc                     <  3:274201/149‑1 (MQ=255)
                                                                                                                                 gTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGc                 >  3:54592/1‑149 (MQ=255)
                                                                                                                                      tCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGTGGCACCAGCTGTGGCTGCACCCGCTCCgg            >  4:206784/1‑149 (MQ=255)
                                                                                                                                              ccaccGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTg    <  5:130149/149‑1 (MQ=255)
                                                                                                                                               caccGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctcc                                                                                                          >  2:340884/1‑46 (MQ=255)
                                                                                                                                               caccGATGGCTTACTGAAACTGAAAGTGAAAACGAActccagctcc                                                                                                          <  1:340884/46‑1 (MQ=255)
                                                                                                                                               caccGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTga   <  8:253662/149‑1 (MQ=255)
                                                                                                                                               caccGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTg    >  7:276648/1‑148 (MQ=255)
                                                                                                                                                accGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTgaa  <  8:276648/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCCGATCCGGTTGCTCGTCATACCACCGATGGCTTACTGAAACTGAAAGTGAAAACGAACTCCAGCTCCAGCGTGTTGGTAGGACCTTTATTTGGTTCTTCCGCTCCAGCTCCGGTTACGGTAGGTAACACGGCGGCACCAGCTGTGGCTGCACCCGCTCCGGCACCGGTGAA  >  NC_000913/4250010‑4250302

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: