Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,358 T→G 11.8% intergenic (+50/‑44) ybaV → / → fadM putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913464,3580TG11.8% 136.9 / 8.0 65intergenic (+50/‑44)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base G (5/3);  total (37/28)
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.84e-01

TACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGTT  >  NC_000913/464216‑464505
                                                                                                                                              |                                                                                                                                                   
tACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTATCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAc                                                                                                                                               >  1:820411/1‑149 (MQ=255)
     tCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTg                                                                                                                                          <  2:898592/149‑1 (MQ=255)
     tCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGATATGACCTCTg                                                                                                                                          >  1:804198/1‑149 (MQ=255)
     tCCGTTTAAAACTGGGGAGGATCTAATGCAGGTGCCGGGGAGGGGCTATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGTCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctctt                                                                                                                                          <  2:1025144/149‑2 (MQ=255)
      ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                                                                         <  2:474129/149‑1 (MQ=255)
       cGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                        <  2:1010275/149‑4 (MQ=255)
       cGTTTAAAACTGGGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                                                                        <  2:815105/149‑1 (MQ=255)
        gTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGTCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                                                                       <  2:734591/149‑1 (MQ=255)
                  gtgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaac                                                                                                                             >  1:733960/1‑149 (MQ=255)
                   tgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                                                                            >  1:550585/1‑149 (MQ=255)
                     gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGTCCTCTGTACTTATaacaacaac                                                                                                                          >  1:204516/1‑149 (MQ=255)
                       ggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGCCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGt                                                                                                                        >  2:986431/1‑149 (MQ=255)
                          tCTAAAGCCGGTGCAGGGCATGGGCAATTCGCTGGTTGAACGTAATCTGGAGGTATTAACCCTGTAATTAATTTGCATACTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGAACTTATAACAACAACGTAAg                                                                                                                     <  2:650530/149‑1 (MQ=255)
                                gCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTAtt                                                                                                               <  2:114133/149‑1 (MQ=255)
                                      gCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                         <  1:732626/149‑1 (MQ=255)
                                        cGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATg                                                                                                       >  2:588839/1‑149 (MQ=255)
                                             attggaaaTTCGCTGGGGGAACGTAATCTGGCGGTATTTACCCTGTAATTAATTTGCATAGTGGCAATTTTTTTCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                                  <  1:105048/143‑1 (MQ=255)
                                             aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                                  <  1:306761/149‑1 (MQ=255)
                                                 gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                              >  1:724418/1‑149 (MQ=255)
                                                  cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGACAGACTGAATAGGTCAAACCACTTATGACCTCTGTACTTATAACAACAACCTAAGGTTATTGCGCTATGCAAACACaaa                                                                                             >  2:925960/1‑149 (MQ=255)
                                                    aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCCCTATGCAAACACAAATc                                                                                           >  1:613355/1‑149 (MQ=255)
                                                      tCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                         <  2:27654/149‑1 (MQ=255)
                                                        gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                                                                       >  2:364883/1‑149 (MQ=255)
                                                         cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAAATTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGt                                                                                      >  1:130530/1‑149 (MQ=255)
                                                              ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                                 <  1:234857/149‑1 (MQ=255)
                                                                aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGa                                                                               >  2:158263/1‑149 (MQ=255)
                                                                 aCGTAATCTGGCTGTTTTATCCCTGTAATTAATTTGCATAGTGTCAATTTTTTGCCAGAATGAAGATGTCATACCAGTTATTACATATGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTTGat                                                                              <  2:633498/149‑1 (MQ=255)
                                                                   gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGatat                                                                            >  2:606154/1‑149 (MQ=255)
                                                                     aaTCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                                                                          >  1:731268/1‑149 (MQ=255)
                                                                       tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcatc                                                                        <  1:879419/149‑1 (MQ=255)
                                                                                aTTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGttt                                                               >  1:963052/1‑149 (MQ=255)
                                                                                          tatttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCGGTACTTATACCAACCACGTAAGGTTATTGCGATATGCAACCACAACTCATCGTTCGTGGATAGCATCTCGACGTTTCCCGGCACga                                                     >  2:564469/4‑148 (MQ=255)
                                                                                          taattaatTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGATATGACGTCTGTACTTATTACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGt                                                     <  2:804198/149‑1 (MQ=255)
                                                                                           aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTAATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAATTTCGTGGATATCATCTCGACTTTTACCAGCACGTc                                                    >  2:957309/1‑149 (MQ=255)
                                                                                             ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                                  <  1:718038/149‑1 (MQ=255)
                                                                                               aattTGCATATTGGAAAATTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTATGTACTTCTCACAAGAACGTAAGGTTATTGCGCTATGCAAACACAAAGCACCGTTCGTGGATCTCATCTCGACGTTTACCAGCACGTcaaca                                                <  2:197549/149‑1 (MQ=255)
                                                                                                    gCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGcc                                           >  1:795020/1‑149 (MQ=255)
                                                                                                     cctAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAAAAACGTAAGGTTATTGCGCTATGCAAACAGAAATCAAAGTTCGGGGATATCATCTCGACGTTTACCAGCCCGTCAACCACGccc                                          >  2:798035/3‑149 (MQ=255)
                                                                                                      aTAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAACAAAGTTCGTGGATATCATCTCGACGTTTACAAGCAAGTCAACACCGCCCg                                         >  2:906398/1‑149 (MQ=255)
                                                                                                       tAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGc                                        <  1:158263/149‑1 (MQ=255)
                                                                                                              aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                                 >  2:543079/1‑149 (MQ=255)
                                                                                                              aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                                 >  2:639369/1‑149 (MQ=255)
                                                                                                                ttttttGCCACACTGAAGACGTCATACCAGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAGGCAAACACAAATCACAGTTCGTGGATATCATCTCGACGGTTACCAGCACGTCAACAACGCCCGCTACCTTGaa                               <  2:601785/149‑1 (MQ=255)
                                                                                                                 tttttGCCAGACTGAAGAGGGCATACCAGTTATGACCTCTGTACTTATAACAACAACGTACGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAt                              <  2:951423/149‑1 (MQ=255)
                                                                                                                 tttttGCCAGACTGAAGAGGGCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACACATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAAt                              <  2:948046/149‑1 (MQ=255)
                                                                                                                    ttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtc                           <  2:731268/149‑1 (MQ=255)
                                                                                                                      gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtctc                         >  1:9943/1‑149 (MQ=255)
                                                                                                                      gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtctc                         <  2:31723/149‑1 (MQ=255)
                                                                                                                        cAGACTGAAGAGGTCATACCAGTTTTGCCCTCTGTACTTATACCAACAACGTAAGGTTATTGCGCTATGCACCCACACATCAATGTTCGTGGAAATCATCTCGCCGTTTACCAGCACGTCCACCACGCCCGCTCCCTTGAATTTCTCga                       >  2:948196/1‑149 (MQ=255)
                                                                                                                        cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCga                       >  2:521710/1‑149 (MQ=255)
                                                                                                                         aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                      >  2:274854/1‑149 (MQ=255)
                                                                                                                           aCTGAAGAGGTCATACCAGGTATGCCCTCTGGACTTATAACAACAACGGAAGGGTATTGCGCTATCCAAACACAAATCAAAGTTCGTCGAGATCATCTCGCCGTTAACCAGCACGTCAACAACCCCCGCAACCTTCACATTCGCgaagc                    >  2:572621/1‑148 (MQ=255)
                                                                                                                            atGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaagaa                   <  1:5573/148‑1 (MQ=255)
                                                                                                                              acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGc                 >  1:731590/3‑149 (MQ=255)
                                                                                                                               cagagGTCATAACTGGTATGACCTCTGTACTTATAACCACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTCCCAGCCCGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGcc                >  2:459490/2‑149 (MQ=255)
                                                                                                                               cagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATGATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACGTTGAATTTCTCGAAGAAGcc                >  2:226367/2‑149 (MQ=255)
                                                                                                                               aaGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTAGAAGAAGcc                >  1:632193/1‑149 (MQ=255)
                                                                                                                                agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGccc               >  1:353173/1‑149 (MQ=255)
                                                                                                                                    gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAATCAGAAATCAAAGTTCGTGGATATCATCTCGACGTTTAGCAGCTCGTGAACAACGCCCGCTACGTTGAATTTCTCGAAGAAGCCCGCTg           <  1:925960/149‑1 (MQ=255)
                                                                                                                                    gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGATATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTg           >  2:335550/1‑149 (MQ=255)
                                                                                                                                      cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTggg         <  2:733960/149‑1 (MQ=255)
                                                                                                                                       aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGa        >  2:488585/1‑149 (MQ=255)
                                                                                                                                         aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATg      >  1:764174/1‑149 (MQ=255)
                                                                                                                                           cAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACAAGCACGTCAACAACGCCCGCTACCTTGAAGTTCTCGAAGAAGCCCGGTGGGATggg    >  2:867971/1‑149 (MQ=255)
                                                                                                                                             gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGtt  <  1:588839/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGGTT  >  NC_000913/464216‑464505

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: