Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,356 A→T 9.3% intergenic (+48/‑46) 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,3560AT 9.3% 154.9 / 7.6 65intergenic (+48/‑46)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base T (4/2);  total (38/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 = 3.23e-01

GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGG  >  NC_000913/464208‑464503
                                                                                                                                                    |                                                                                                                                                   
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                       <  2:31723/149‑1 (MQ=255)
                                                                                                                              gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtctc                       >  1:9943/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                   
GCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGGG  >  NC_000913/464208‑464503

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: