Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,354 C→A 8.1% intergenic (+46/‑48) 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,3540CA 8.1% 156.0 / 2.7 66intergenic (+46/‑48)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base C (35/25);  new base A (4/2);  total (39/27)
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.88e-01

GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATG  >  NC_000913/464208‑464501
                                                                                                                                                  |                                                                                                                                                   
gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCa                                                                                                                                                   >  2:234857/1‑149 (MQ=255)
gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCa                                                                                                                                                   >  2:879419/1‑149 (MQ=255)
 cgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAg                                                                                                                                                  >  1:179067/1‑149 (MQ=255)
        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)
                                                                                                                                                  |                                                                                                                                                   
GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATG  >  NC_000913/464208‑464501

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: