Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 464,354 C→A 7.5% 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 7.5% 410.1 / 18.0 120intergenic (+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 (60/51);  new base A (3/6);  total (63/57)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 122 total aligned reads displayed.

CGCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGG  >  NC_000913/464207‑464502
                                                                                                                                                   |                                                                                                                                                    
cgcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAc                                                                                                                                                      <  2:176790/148‑1 (MQ=255)
 gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCa                                                                                                                                                    <  3:127787/149‑1 (MQ=255)
   gAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGt                                                                                                                                                  >  4:225751/1‑149 (MQ=255)
     agagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGTCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                                >  4:487508/1‑149 (MQ=255)
     agagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                                <  3:225751/149‑1 (MQ=255)
     agagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                                <  4:634455/149‑1 (MQ=255)
      gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                                <  2:268599/148‑1 (MQ=255)
      gagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                                                                               >  4:123732/1‑149 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                         <  4:226197/149‑1 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                         <  3:8725/149‑1 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                         <  3:53030/149‑1 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                         <  3:211566/149‑1 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                          <  1:804277/148‑1 (MQ=255)
             gTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACctct                                                                                                                                        >  4:29495/1‑149 (MQ=255)
               ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTg                                                                                                                                       >  2:783279/1‑148 (MQ=255)
                   ttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                                                                  <  3:310842/149‑1 (MQ=255)
                   ttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                                                                   >  1:813831/1‑148 (MQ=255)
                   ttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                                                                   <  2:366596/148‑1 (MQ=255)
                   ttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                                                                   >  1:548655/1‑148 (MQ=255)
                    tAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                                                                  >  2:672364/1‑148 (MQ=255)
                       aaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                                                               <  1:877712/148‑1 (MQ=255)
                       aaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                                                                               <  2:810499/148‑1 (MQ=255)
                       aaCTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                              >  4:469461/1‑149 (MQ=255)
                          tgtgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                                                                            <  1:264551/148‑1 (MQ=255)
                           gtgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaa                                                                                                                           <  1:5437/148‑1 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                       >  3:110702/1‑149 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                                                                            >  2:443334/1‑144 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                                                                            <  1:443334/144‑1 (MQ=255)
                                  aTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                                                                   >  4:724766/1‑149 (MQ=255)
                                  aTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                                                                   <  3:29495/149‑1 (MQ=255)
                                   tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                                                                   <  2:699419/148‑1 (MQ=255)
                                   tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                  >  3:96756/1‑149 (MQ=255)
                                   tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                  >  4:355330/1‑149 (MQ=255)
                                   tCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                  >  4:567548/1‑149 (MQ=255)
                                    cTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                  >  2:549176/1‑148 (MQ=255)
                                    cTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                                                                  >  1:774669/1‑148 (MQ=255)
                                      aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGtt                                                                                                               <  3:156990/149‑1 (MQ=255)
                                      aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGt                                                                                                                >  2:445152/1‑148 (MQ=255)
                                              tGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCTAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCt                                                                                                       <  3:487508/149‑1 (MQ=255)
                                                ccGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                                                                      >  2:50139/1‑148 (MQ=255)
                                                ccGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                                                                     <  3:57147/149‑1 (MQ=255)
                                                  ggggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGc                                                                                                   >  3:101472/1‑149 (MQ=255)
                                                     gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                                <  4:317991/149‑1 (MQ=255)
                                                     gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                                 <  1:893085/148‑1 (MQ=255)
                                                     gATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                                 <  1:934466/148‑1 (MQ=255)
                                                           cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaa                                                                                           <  2:164246/148‑1 (MQ=255)
                                                            aaTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                                                                          >  2:868988/1‑148 (MQ=255)
                                                                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaaa                                                                                     <  3:922764/149‑1 (MQ=255)
                                                                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                      <  2:52300/148‑1 (MQ=255)
                                                                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                                                                    <  1:912908/148‑1 (MQ=255)
                                                                    ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                                                                  >  1:621310/1‑148 (MQ=255)
                                                                    ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                                                                 <  4:110702/149‑1 (MQ=255)
                                                                       ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                                                              <  4:291076/149‑1 (MQ=255)
                                                                          aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaac                                                                                                                       <  3:63702/105‑1 (MQ=255)
                                                                            gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCACACCAGTTATGACctct                                                                                                                                        >  4:677857/1‑86 (MQ=255)
                                                                            gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCACACCAGTTATGACctct                                                                                                                                        <  3:677857/86‑1 (MQ=255)
                                                                                tCTGGCGGTATCAACCCTGTAAGTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGACATACCAGTTATGAACTCCGTACTTAAAGCATCAACGTTAGGTCATAGCGCTATGCAAACACAAACCAAAGATCGTGGATaccat                                                                      <  2:348429/148‑1 (MQ=255)
                                                                                   ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                                                                   >  2:497967/1‑148 (MQ=255)
                                                                                      ggTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGCCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                                                               >  3:577104/1‑149 (MQ=255)
                                                                                       gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg                                                               >  2:247428/1‑148 (MQ=255)
                                                                                        tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGtt                                                             <  3:666345/149‑1 (MQ=255)
                                                                                        tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGtt                                                             >  3:423841/1‑149 (MQ=255)
                                                                                           tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTa                                                           <  2:248326/148‑1 (MQ=255)
                                                                                           tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTa                                                           >  2:802146/1‑148 (MQ=255)
                                                                                           tAACCATGTAATTAATTTGCATAGTGGAAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTATGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACGAATCAACGTTTGTCGATCTCATCTCGCCGttga                                                           <  2:938746/148‑3 (MQ=255)
                                                                                            aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc                                                          <  1:497967/148‑1 (MQ=255)
                                                                                            aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc                                                          >  2:422227/1‑148 (MQ=255)
                                                                                             aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc                                                         <  1:672364/148‑1 (MQ=255)
                                                                                                  gTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc                                                    >  2:255151/1‑148 (MQ=255)
                                                                                                    aattaattTGCATAGTGGCAATTTTTTGCCAGCCTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACAAAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGt                                                  >  1:933113/1‑148 (MQ=255)
                                                                                                       taattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaac                                              <  3:355330/149‑1 (MQ=255)
                                                                                                       taattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                               >  1:8690/1‑148 (MQ=255)
                                                                                                          tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaacaa                                            <  1:767657/148‑1 (MQ=255)
                                                                                                           ttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCaacaac                                           >  2:619457/1‑148 (MQ=255)
                                                                                                              cATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGcc                                        <  1:183016/148‑1 (MQ=255)
                                                                                                               aTAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaaa                                                                                     >  3:920592/1‑102 (MQ=255)
                                                                                                                  gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTa                                   >  3:885062/1‑149 (MQ=255)
                                                                                                                       aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCtt                               >  1:728198/1‑148 (MQ=255)
                                                                                                                       aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTg                              >  3:880192/1‑149 (MQ=255)
                                                                                                                                ccAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCg                     <  4:750855/149‑1 (MQ=255)
                                                                                                                                 cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCg                     >  2:60106/1‑148 (MQ=255)
                                                                                                                                  aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCga                    <  1:540086/148‑1 (MQ=255)
                                                                                                                                   gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc                                                                <  3:62002/103‑1 (MQ=255)
                                                                                                                                   gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc                                                                >  4:62002/1‑103 (MQ=255)
                                                                                                                                   gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaag                  >  3:851419/1‑149 (MQ=255)
                                                                                                                                   gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaa                   >  1:760416/1‑148 (MQ=255)
                                                                                                                                     cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaaga                 <  1:619457/148‑1 (MQ=255)
                                                                                                                                       gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGc              <  3:724766/149‑1 (MQ=255)
                                                                                                                                        aaGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGcc             >  4:249670/1‑149 (MQ=255)
                                                                                                                                         agagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGcc             <  1:579424/148‑1 (MQ=255)
                                                                                                                                          gagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGccc            >  2:582733/1‑148 (MQ=255)
                                                                                                                                           agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg           >  1:262130/1‑148 (MQ=255)
                                                                                                                                            ggTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGc          >  2:871316/1‑148 (MQ=255)
                                                                                                                                            ggTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGc          >  2:177549/1‑148 (MQ=255)
                                                                                                                                             gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTg        >  3:151413/1‑149 (MQ=255)
                                                                                                                                               cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAATCCCGCTgg       >  2:791576/1‑148 (MQ=255)
                                                                                                                                                aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGa     >  4:515474/1‑149 (MQ=255)
                                                                                                                                                  aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGAt    <  1:445152/148‑1 (MQ=255)
                                                                                                                                                  aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATg   >  4:38151/1‑149 (MQ=255)
                                                                                                                                                   ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATgg  <  4:851419/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGGGATGG  >  NC_000913/464207‑464502

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: